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1. 


Foreword 


This  final  report  describes  results  obtained  with  support  from  ARO  grant  DAAH04-96-I- 
0096  between  05/0 1  /96  and  04/30/00.  This  grant  has  made  considerable  impact  in  that  it  has  enabled 
us  to  expand  and  diversify  our  research  program.  The  principal  thrust  of  work  supported  by  this  grant 
focused  on  molecular  mechanisms  of  chemoreception  in  the  mammalian  vomeronasal  organ.  We 
characterized  the  transduction  pathway  for  the  recognition  of  pheromones  in  the  vomeronasal  organ 
and  also  characterized  subpopulations  of  olfactory  neurons  expressing  different  axonal  G  proteins 
in  the  main  and  accessory  olfactory  projections.  These  observations  bear  directly  on  the  formation 
of  chemotopic  sensory  maps  in  the  brain.  These  studies  were  recently  extended  with  a 
characterization  of  genetically  modified  mice  deficient  in  the  CC  subunit  of  the  G-protein,  Go,  which 
show  profound  anosmia.  These  studies  resulted  in  publications  in  Endocrinology  and  Brain  Research 
and  have  also  been  reported  in  abstract  form  at  several  national  and  international  meetings.  This 
work  was  performed  with  the  involvement  of  postdoctoral  research  associate,  Dr.  Kennedy  S. 
Wekesa,  currently  an  Assistant  Professor  at  Alabama  State  University,  Montgomery,  AL.  This  work 
is  currently  being  continued  with  support  from  the  W.  M.  Keck  Foundation. 

In  addition  to  studies  on  the  mammalian  olfactory  system,  grant  DAAH04-96-I-0096  enabled 
us  to  expand  work  on  the  genomic  architecture  of  olfactory  behavior  in  Drosophila  melanogaster, 
which  is  currently  continued  with  support  from  a  grant  from  the  National  Institute  of  General 
Medical  Sciences.  This  work  led  to  the  discovery  of  several  new  gene  products  that  are  essential  for 
processing  olfactory  information,  including  a  sodium  channel,  a  novel  dual  specificity  tyrosine 
phosphorylation  regulated  kinase  ( Dyrk2 ),  a  postsynaptic  density  protein  (Scribble),  an  olfactory 
receptor,  and  an  odorant  binding  protein,  the  latter  two  interacting  with  the  repellent  odorant, 
benzaldehyde.  These  studies  have  been  published  in  Genetics  and  documented  in  abstract  form. 

Finally,  we  were  able  to  continue  work  on  olfactomedin.  Olfactomedin  was  originally 
identified  as  the  major  mucus  component  of  the  olfactory  neuroepithelium.  Subsequently,  other 
investigators  identified  an  olfactomedin  homologue  that  is  expressed  as  different  splice  variants 
throughout  the  rat  brain.  These  proteins  became  especially  important  when  an  olfactomedin-related 
protein  was  discovered  in  the  anterior  segment  of  the  eye  and  found  to  be  closely  associated  with 
the  pathogenesis  of  glaucoma  (the  trabecular  meshwork  inducible  glucocorticoid  response  protein, 
TIGR).  Our  laboratory  has  shown  that  olfactomedin-related  proteins  occur  ubiquitously  from  C. 
elegans  to  humans  and  that  they  have  conserved  C-terminal  domains  with  characteristic  sequence 
motifs.  Moreover,  we  showed  that  a  family  of  olfactomedin-related  proteins  is  encoded  in  the 
human  genome  with,  thus  far,  at  least  five  members.  These  olfactomedin-related  proteins  appear  to 
be  members  of  a  diverse  family  of  tissue-specific  extracellular  matrix  components.  The  link  between 
TIGR  and  ocular  hypertension  and  the  expression  of  several  of  these  proteins  in  mucus-lined  tissues 
suggest  that  they  play  an  important  role  in  regulating  physical  properties  of  the  extracellular 
environment.  Our  recent  discovery  of  an  olfactomedin  homologue  in  Drosophila  melanogaster  will 
facilitate  studies  on  the  functions  of  olfactomedin-related  proteins  and  their  interactions.  W ork  on 
olfactomedin  was,  in  part,  performed  by  postdoctoral  research  associate.  Dr.  Christa  Karavanich, 
currently  a  faculty  member  at  Richland  College,  Dallas,  TX.  This  work  has  been  published  in 
Molecular  Biology  and  Evolution,  in  Genetical  Research  and  as  several  abstracts,  and  is  currently 
supported  by  the  Glaucoma  Research  Foundation. 
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4.  Body  of  Report 


A.  Statement  of  the  Problems  Studied 

(i)  Molecular  and  cellular  mechanisms  of  chemoreception  in  the  mammalian 
olfactory  system 

Odor  recognition  is  essential  for  the  survival  and  procreation  of  most  animals.  Chemical 
information  transferred  via  olfactory  neurons  to  the  olfactory  bulb  is  transformed  into  a  chemotopic 
map  (Joerges  et  al. ,  1 997 ;  Friedrich  and  Korsching,  1 997  and  1998;  Bozza  and  Kauer,  1 998 ;  Rubin 
and  Katz,  1999).  This  map  is  represented  by  modular  neural  arrays,  glomeruli,  which  in  the  mouse 
number  approximately  1,800  (Pomeroy  et  al.,  1990)  and  each  of  which  represents  the  convergent 
projection  of  neurons  of  similar  chemosensory  specificity  (reviewed  by  Hildebrand  and  Shepherd, 
1997). 

The  olfactory  system  contains  two  interacting  components:  the  main  olfactory  system  is 
dedicated  to  general  odorant  discrimination,  whereas  the  accessory  olfactory  system  primarily 
processes  chemical  cues  that  guide  social  behaviors  (reviewed  by  Meredith,  1998).  In  the  main 
olfactory  system  odor  recognition  is  mediated  by  receptors,  which  belong  to  the  superfamily  of  G- 
protein  coupled  receptors  and  are  encoded  by  as  many  as  1,000  different  genes  (Buck  and  Axel, 
1991;  Levy  et  al,  1991;  Strotmann  et  al.,  1994;  Axel,  1995;  Buck,  1996;  Sullivan  et  al,  1995  and 
1996).  Two  distinct  families  of  G  protein  coupled  receptors  have  also  been  identified  as  putative 
pheromone  receptors  in  chemosensory  neurons  of  the  vomeronasal  organ  (VNO),  the  chemosensory 
organ  of  the  accessory  olfactory  system  (Dulac  and  Axel,  1 995 ;  Matsunami  and  Buck,  1 997;  Herrada 
and  Dulac,  1997;  Ryba  and  Tirindelli,  1997). 

In  the  main  olfactory  system  each  olfactory  neuron  expresses  a  single  odorant  receptor  from 
the  large  repertoire  of  receptors  encoded  in  the  genome  (Chess  et  al.,  1994).  Thus,  binding  of  an 
odorant  to  its  receptor  translates  directly  into  a  distinct  pattern  of  neural  activity  that  encodes  its 
structure.  Receptor  activation  results  in  stimulation  of  adenylate  cyclase  (Pace  et  al.,  1985;  Sklar  et 
al. ,  1 986;  reviewed  by  Anholt,  1993),  followed  by  the  opening  of  cyclic  nucleotide-gated  channels 
that  carry  the  generator  current  (Nakamura  and  Gold,  1 987;  Firestein  et  al.  1991).  Calcium,  entering 
through  these  channels  (Frings  et  al.,  1995),  amplifies  the  response  by  opening  calcium-gated 
chloride  channels  (Kleene  and  Gesteland,  1991 ;  Kurahashi  and  Yau,  1993;  Reuter  et  al.,  1998)  and, 
after  binding  to  calmodulin,  limits  the  response  by  lowering  the  affinity  of  the  cyclic  nucleotide- 
activated  channel  for  cyclic  AMP  (Chen  and  Yau,  1994;  Kurahashi  and  Menini,  1997).  Although 
inositol- 1,4, 5, -triphosphate  also  has  been  implicated  as  second  messenger  in  olfaction  (Boekhoff  et 
al.  ,1990;  Ronnett  et  al. ,  1 993),  knock-out  mice  deficient  in  the  a  subunit  of  the  cyclic  nucleotide¬ 
gated  channel  were  found  to  be  generally  anosmic  (Brunet  et  al. ,  1 996)  as  were  mice  deficient  in  Golf 
(Belluscio  et  al. ,  1 998),  the  G  protein  that  links  odorant  receptor  activation  to  stimulation  of  adenylyl 
cyclase  (Jones  and  Reed,  1989). 

Whereas  cyclic  AMP  appears  to  be  the  primary  second  messenger  in  olfactory  transduction, 
inositol- 1,4, 5, -triphosphate has  been  implicated  in  signal  transduction  in  the  VNO  (Jiang  et  al. ,  1990; 
Luo  etal.,  1994:  Taniguchi  etal,  1995;  Wekesa  and  Anholt,  1997),  although  a  role  for  cyclic  AMP 
has  also  been  suggested  in  this  system  (Jiang  et  al.,  1990;  Luo  et  al.,  1994;  Okamoto  et  al.,  1996). 


5 


Our  laboratory  played  a  major  role  in  elucidating  the  second  messenger  pathway  that  operates  in  the 
VNO  by  obtaining  direct  measurements  of  second  messengers  in  membrane  preparations  from 
porcine  VNO  and  demonstrating  dose-dependent,  tissue-specific  and  sex-dependent  generation  of 
inositol- 1 ,4, 5, -triphosphate  in  VNO  membranes  from  prepubertal  female  pigs  in  response  to  seminal 
fluid  and  boar  urine  (Wekesa  and  Anholt,  1997).  This  generation  of  inositol- 1,4, 5, -triphosphate  is 
GTP-dependent  and  is  thought  to  be  mediated,  at  least  in  part,  through  a  Gq/ll  -related  G-protein 
(Wekesa  and  Anholt,  1997). 

The  dendritic  and  axonal  compartments  of  olfactory  neurons  fulfill  distinct  functions  for  the 
acquisition  of  chemosensory  information.  Whereas  dendritic  specializations  mediate  odorant 
recognition  and  chemosensory  transduction,  the  axonal  compartment  regulates  signal  propagation, 
axon  sorting  and  target  innervation.  The  roles  of  G  proteins  in  chemosensory  transduction  at  the 
dendritic  compartments  of  chemosensory  neurons  have  been  studied  extensively,  but  the  functions 
of  axonal  G  proteins  have  not  been  investigated  in  detail. 

An  important  role  for  axonal  G  proteins  in  mammalian  chemoreception  is  implicated  by  the 
differential  distribution  of  Gi2  and  G0  in  the  accessory  olfactory  system.  The  apical  layer  of  the  VNO 
expresses  a  distinct  family  of  putative  pheromone  receptors  (VN 1  receptors;  Dulac  and  Axel,  1 995), 
whereas  VNO  neurons  in  the  basal  layer  express  a  different  family  of  receptors  (VN2  receptors)  that 
resemble  metabotropic  glutamate  receptors  (Matsunami  and  Buck,  1997;  Herrada  and  Dulac,  1997; 
Ryba  and  Tirindelli,  1997).  In  situ  hybridization  studies  in  the  murine  VNO  (Berghard  and  Buck, 
1 996)  and  immunohistochemical  studies  in  opossum  (Halpern  et  al. ,  1 995)  showed  that  neurons  in 
the  apical  layer  of  the  VNO  express  Gi2  and  project  to  the  rostral  region  of  the  accessory  olfactory 
bulb  (AOB),  whereas  neurons  in  the  basal  layer  express  G0  and  project  to  the  caudal  region  of  the 
AOB.  Our  laboratory  has  demonstrated  that  differential  expression  patterns  of  Gi2  and  G0  are  not 
unique  to  the  AOB,  but  also  occur  in  the  main  olfactory  bulb  (MOB;  Wekesa  and  Anholt,  1999). 


(ii)  Functional  genomics  of  odor-guided  behavior 

Behavior  supports  the  stage  on  which  the  interplay  between  the  environment  and  the  genome 
guides  evolution.  Behavioral  adaptations  to  environmental  cues  determine  survival  and  reproductive 
success.  Among  those  environmental  cues  none  are  more  important  than  chemical  signals.  Indeed, 
odor-guided  behavior  is  absolutely  essential  for  most  organisms  for  food  localization,  avoidance  of 
environmental  toxins  or  predators,  oviposition  site  selection,  kin  recognition,  species  recognition  and 
mate  selection,  and  reproduction.  The  importance  of  molecular  recognition  of  environmental 
chemicals  is  underscored  by  the  discovery  that  an  unusually  large  percentage  of  the  genome  is 
dedicated  to  olfaction:  in  the  mammalian  genome,  odor-recognition  alone  is  mediated  by  a  multigene 
family  of  about  1,000  odorant  receptors  (Buck  and  Axel,  1991),  that  has  undergone  rapid  evolution 
through  gene  duplication  and  diversification  (Hughes  and  Hughes,  1 993 ;  Ben-Arie  et  al. ,  1 993 ;  Issel- 
Tarver  and  Rine,  1997;  Rouquier  et  al,  1998).  Initial  chemosensory  recognition  events,  however, 
are  only  the  first  step  in  a  series  of  complex  processes,  that  involves  processing  of  chemosensory 
information  in  the  central  nervous  system,  evaluating  the  nature  of  the  chemosensory  perception,  and 
initiating  and  executing  an  appropriate  behavioral  response.  It  is  intuitively  evident  that  this  process 
involves  a  vast  ensemble  of  genes  and  that  variation  in  the  expression  of  any  of  these  genes  can 
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generate  individual  variation  in  olfactory  responsiveness  within  a  natural  population.  Thus,  odor- 
guided  behavior  is  a  quantitative  trait  and  understanding  its  complex  genetic  architecture  requires 
quantitative  genetic,  statistical,  and  genomic  analyses. 

We  have  chosen  olfactory  avoidance  behavior  of  Drosophila  melanogaster  as  a  model  system 
for  studies  on  the  genetic  architecture  of  odor-guided  behavior.  D.  melanogaster  is  a  genetic  model 
system  par  excellence ,  because  its  generation  time  is  short,  there  is  no  recombination  in  males,  and 
the  availability  of  balancer  chromosomes  allows  mutations  to  be  stably  propagated  and  enables  the 
manipulation  of  entire  chromosomes  to  construct  “designer  genotypes”.  Furthermore,  highly  inbred 
lines  can  be  readily  generated,  which  eliminates  genetic  variance  and  greatly  facilitates  the  analysis 
of  complex  traits.  In  addition,  extensive  resources  are  available  for  genomic  studies  on  Drosophila, 
including  deficiency  lines  and  lines  that  carry  transposable  element  insertions  at  different  defined 
locations  in  the  genome;  these  resources  are  available  through  large  national  stock  centers. 
Sequencing  of  the  entire  genome  has  recently  been  completed  (Adams  et  al.,  2000)  and  a  large 
genomic  data  base  is  available  (http://flybase.bio.indiana.edu/). 

There  are  additional  reasons  that  make  Drosophila  eminently  suited  for  studies  on  olfactory 
behavior,  as  described  in  more  detail  below.  The  functional  organization  of  its  olfactory  system  is 
similar  to  that  of  the  vertebrate  olfactory  system,  suggesting  similar  principles  of  odor  coding,  but 
it  is  more  tractable  since  it  contains  about  10,000-100,000  times  fewer  neurons  and  30  times  fewer 
glomeruli  (Laissue  et  al.,  1999).  Furthermore,  D.  melanogaster  is  the  only  insect  species  to  date  in 
which  odorant  receptors  have  been  identified  (Clyne  et  al.,  1999a;  Vosshall  et  al.,  1999).  Finally  - 
and  most  importantly  -  the  onset  of  avoidance  responses  to  repellent  odorants  can  be  precisely 
controlled  and  a  reproducible  statistical  behavioral  assay  has  been  devised  to  accurately  measure 
small  olfactory  impairments  quantitatively  (Anholt  et  al,  1996). 

(Hi)  Characterization  of  olfactomedin-related proteins 

Olfactomedin  was  originally  identified  as  the  major  mucus  component  of  the  olfactory 
neuroepithelium  (Snyder  et  al.,  1991;  Yokoe  and  Anholt,  1993).  Subsequently,  other  investigators 
identified  an  olfactomedin  homologue  that  is  expressed  as  different  splice  variants  throughout  the 
rat  brain  (Danielson  etal. ,  1 994).  These  proteins  became  especially  important  when  an  olfactomedin- 
related  protein  was  discovered  in  the  anterior  segment  of  the  eye  and  found  to  be  closely  associated 
with  the  pathogenesis  of  glaucoma  (the  trabecular  meshwork  inducible  glucocorticoid  response 
protein,  TIGR;  Stone  et  al.,  1997;  Polanski  et  al.,  1997;  Kubota  et  al.,  1997).  We  characterized  the 
family  of  olfactomedin-related  proteins  both  in  the  human  genome  and  in  Drosophila  to  set  the  stage 
for  studies  aimed  at  determining  their  functions.  These  studies  involved  a  collaboration  with  Dr. 
William  Atchley  from  the  Department  of  Genetics  and  with  Dr.  Teresa  Borras  from  the  Department 
of  Ophthalmology  at  Duke  University  Medical  Center. 


B.  Summary  of  the  Most  Important  Results 

(i)  Molecular  and  cellular  mechanisms  of  chemoreception  in  the  mammalian 
olfactory  system 

After  many  years  of  previous  research  on  signal  transduction  in  the  main  olfactory  system, 
we  turned  our  attention  to  pheromonal  signal  transduction  in  the  vomeronasal  organ  (VNO),  a 
chemosensory  organ  specialized  for  the  perception  of  social  chemical  signals.  We  identified  the  pig 
as  a  good  model  system  for  biochemical  studies  because  of  the  large  size  of  its  YNOs.  We  developed 
a  membrane  preparation  enriched  in  microvillar  membranes  from  chemosensory  vomeronasal 
neurons.  These  preparations  allowed  us  to  measure  for  the  first  time  directly  pheromone-regulated 
increases  in  the  second  messenger,  inositol  trisphosphate,  in  the  mammalian  VNO.  We  demonstrated 
in  membranes  from  prepubertal  female  animals  sex-specific,  dose-dependent  and  G-protein  mediated 
increases  in  inositol  trisphosphate  in  response  to  boar  urine  and  seminal  fluid.  Furthermore,  our 
studies  implicated  Gq/U  as  one  of  the  G-proteins  likely  to  link  pheromone  receptor  activation  to 
stimulation  of  phospholipase  C.  These  experiments  provided  for  the  first  time  a  biochemical  assay 
for  the  activity  of  mammalian  pheromones  and  consolidated  through  direct  measurements  the  notion 
that  the  main  and  accessory  olfactory  systems  use  different  signal  transduction  mechanisms  to  evoke 
neural  activation  (Wekesa  and  Anholt,  1997).  Thus,  in  contrast  to  the  main  olfactory  system  where 
cyclic  AMP  is  the  principal  neurotransmitter,  inositol  trisphosphate  is  the  major  neurotransmitter  in 
the  VNO,  where  it  most  likely  mediates  influx  of  calcium  via  the  opening  of  a  recently  identified 
VNO-specific  TRP2  ion  channel. 

Whereas  odorant  recognition  and  chemosensory  transduction  occur  at  the  dendrites  of 
olfactory  neurons,  signal  propagation,  axon  sorting  and  target  innervation  are  functions  of  their 
axons.  Previous  studies  by  other  investigators  had  shown  that  in  the  accessory  olfactory  system 
vomeronasal  neurons,  whose  cell  bodies  are  located  in  the  superficial  layers  of  the  vomeronasal 
organ,  send  their  axons  to  the  anterior  region  of  the  accessory  olfactory  bulb  and  express  Gj2, 
whereas  neurons  in  the  basal  layer  of  the  vomeronasal  neuroepithelium  send  their  axons  to  the 
posterior  accessory  olfactory  bulb  and  express  G0  (Halpern  et  al.,  1995;  Berghard  and  Buck,  1996) 
We  showed  that  differential  expression  of  these  G  proteins  is  not  unique  to  the  accessory  olfactory 
system,  but  also  occurs  in  the  main  olfactory  projection.  Whereas  in  the  main  olfactory  projection 
all  neurons  express  G0,  a  subpopulation  expresses  Gj2.  We  showed  that  the  projection  of  Gi2- 
expressing  neurons  originates  from  the  dorsomedial  region  of  the  nasal  cavity  and  fans  out  across 
the  medial  aspect  of  the  main  olfactory  bulb,  sparing  the  penetration  corridor  of  the  vomeronasal 
nerve  and  extending  sparsely  to  scattered  glomeruli  in  the  lateral  olfactory  bulb.  These  results  are 
important  in  that  olfactory  neurons  expressing  the  same  odorant  receptor  have  been  shown  to 
converge  on  two  glomeruli  (spherical  synaptic  neural  output  modules)  in  the  olfactory  bulb,  one 
medial  and  one  lateral.  Thus,  the  olfactory  projection  forms  two  chemotopic  representations  in  each 
olfactory  bulb,  one  lateral  map  and  one  medial  map.  Our  results  suggest  that  Gi2  may  play  a  role  in 
signal  processing  or  axonal  targeting  to  the  medial  chemotopic  map.  In  collaboration  with  Dr.  John 
Vandenbergh,  we  are  using  homologous  recombinant  mice  that  lack  the  CC  subunit  of  G0  (V alenzuela 
et  al. ,  1997)  to  further  explore  the  functions  of  these  axonal  G-proteins  in  the  olfactory  projection. 
The  original  knock-out  mice  were  generously  provided  by  Dr.  Eva  Neer  from  Harvard  University 
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as  C57B1/129  chimeras.  We  found  that  intercrossing  heterozygotes  of  these  mice  generated  few 
homozygous  G0-/-  offspring,  presumably  due  to  inbreeding  depression.  Therefore,  we  backcrossed 
these  mice  into  the  CD-I  outbred  genetic  background.  This  increased  the  number  of  homozygous 
knock-out  offspring  and  their  viability.  Homozygous  knock-outs  are  slightly  smaller  in  body  weight 
compared  to  heterozygous  littermates  at  birth,  but  appear  otherwise  normal,  as  reported  by  Neer  and 
colleagues  (Valenzuela  et  al.,  1997).  Whereas  virtually  all  of  the  heterozygotes  and  wild-type 
progeny  were  able  to  locate  a  buried  food  pellet  after  an  overnight  starvation  period,  the  G „-/-  mice 
were  less  successful  in  localizing  the  pellet  within  the  4  min  assay  period.  Olfactory  ability  of  the 
G0-/-  mice  was  further  assessed  by  an  olfactory  habituation-dishabituation  test.  When  control 
heterozygous  mice  are  presented  with  an  unfamiliar  odor  on  a  cotton  wool  swab  protruding  from  the 
cage  lid,  they  rear  up  to  investigate  the  odor.  The  number  of  rearings  and  total  rearing  time  abates 
as  they  habituate  to  the  odor,  but  can  be  elicited  de  novo  by  introducing  a  different  odor.  The  number 
of  rearings  in  response  to  sequential  odor  exposures  was  markedly  lower  or  virtually  absent  in  G0-/- 
mice,  although  their  mobility  was  not  impaired.  Preliminary  studies  further  indicate  that  absence  of 
G0  may  be  accompanied  by  overexpression  of  Gi2.  Although  upon  superficial  examination  olfactory 
bulbs  in  the  knock-out  animals  appear  intact,  closer  examination  indicates  that  the  packing  and  shape 
of  glomeruli  in  these  animals  may  be  altered.  Dr.  Peter  Mombaerts  from  Rockefeller  University  has 
agreed  to  provide  us  with  his  P2-IRES-tau-lacZ  knock-in  mice,  in  which  the  convergent  projection 
of  olfactory  neurons  expressing  the  P2  odorant  receptor  can  be  readily  visualized  by  staining  for  P- 
galactosidase  (Mombaerts  et  al.,  1996).  Introducing  the  P2-IRES-tau-lacZ  knock-in  construct  into 
our  G0-/-  mice  will  enable  us  to  assess  directly  whether  the  pattern  of  olfactory  projections  in  these 
mice  has  been  altered.  In  the  future  we  will  continue  to  investigate  the  impact  of  deletion  of  G0  on 
the  formation  of  chemotopic  projections  and/or  compensatory  expression  of  other  G-proteins. 

(i i)  Functional  genomics  of  odor-guided  behavior 

1.  Behavioral  assay 

Despite  considerable  advances  in  our  understanding  of  the  functional  organization  of  the 
olfactory  system  of  Drosophila,  the  link  between  odor  perception  at  the  molecular  and  cellular  level 
and  odor-guided  behavior  at  the  level  of  the  organism  remains  to  be  explored.  When  exposed  to  a 
strong  repellent  odorant,  flies  will  immediately  migrate  away  from  the  odor  source.  This  ability  to 
avoid  odor  sources  associated  with  toxic  compounds  in  the  environment  obviously  has  important 
survival  value.  To  investigate  the  genetic  basis  for  variation  in  odor-guided  behavior,  we  made  use 
of  this  natural  avoidance  response  and  developed  a  behavioral  paradigm  that  is  rapid,  simple  and 
readily  quantifiable  (Anholt  et  al.,  1996).  Moreover,  the  use  of  isogenic  strains  of  flies  eliminates 
within-strain  genetic  variance.  All  observed  variation  in  measurements  between  individuals  of  the 
same  strain  is  attributable  to  environmental  variation.  Controlling  the  background  genotype  precisely 
increases  the  statistical  power  to  detect  differences  in  behavior  between  different  genotypes. 

We  have  used  benzaldehyde,  a  standard  repellent  odorant,  to  measure  olfactory  responses  in 
flies.  We  have  also  used  structurally  unrelated  odorants,  such  as  2-isobutylthiazole  and  2-n- 
propylpyrazine,  to  test  whether  olfactory  defects  of  lines  with  aberrant  olfactory  behavior  are  specific 
to  benzaldehyde.  Prior  to  measuring  behavioral  responses  to  odorants,  single  sex  groups  of  five 
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individuals  are  placed  in  test  vials  without  food  for  2  hours.  The  test  vials  are  divided  into  two 
compartments  by  placing  a  mark  on  the  wall  3  cm  from  the  bottom  of  the  vial.  The  animals  are  then 
exposed  to  an  aqueous  solution  containing  the  desired  concentration  of  benzaldehyde.  The  odorant 
is  introduced  into  the  vial  on  a  cotton  swab  and  the  number  of  flies  migrating  to  a  compartment 
remote  from  the  odor  source  is  measured  at  5  second  intervals,  from  15  seconds  to  60  seconds  after 
introduction  of  the  odor  source  (Fig.  1).  Distilled  water  is  used  as  a  control.  The  “avoidance  score” 
of  the  replicate  is  the  average  of  these  10  counts,  giving  a  possible  range  of  avoidance  scores 
between  0  (maximal  attraction  to  the  odor  source)  and  5  (all  flies  are  in  the  compartment  away  from 
the  odor  source  for  the  entire  assay  period,  i.e.  a  maximal  repellent  response).  Many  replicate  assays 
are  done  for  each  line,  each  consisting  of  5  males  and  5  females.  The  elimination  of  genetic  variance 
through  the  use  of  an  isogenic  genetic  background  together  with  our  ability  to  rapidly  accumulate 
large  data  sets  for  each  line  through  repeated  measurements  provides  us  with  the  statistical  power 
to  reproducibly  resolve  not  only  large  phenotypic  effects,  but  also  small  smell  impairments  (Anholt 
et  al.,  1996). 

2.  Responding  to  chemosensory  information:  Odor-guided  behavior  as  a  quantitative  trait 

Can  variation  in  odor-guided  behavior  in  a  natural  population  be  resolved  and  quantified? 
If  so,  how  much  of  this  variation  is  due  to  genetic  variance  and  how  much  is  accounted  for  by 
environmental  variance?  How  is  genetic  variation  for  olfactory  behavior  maintained?  What  is  the 
relationship  between  variation  in  odor-guided  behavior  and  variation  in  fitness?  To  address  these 
questions,  chromosomes  were  extracted  from  a  natural  population  and  substituted  into  a  common 
inbred  genetic  background.  The  avoidance  responses  to  benzaldehyde  were  quantified  for  43  X 
chromosome  substitution  lines  and  35  third  chromosome  substitution  lines  (Mackay  et  al.,  1996). 
There  was  significant  genetic  variation  in  avoidance  scores  in  this  sample  of  chromosomes. 
Estimates  of  quantitative  genetic  parameters  showed  that  heritabilities  of  olfactory  avoidance 
behavior  were  low  (averaged  for  both  sexes  h2  =  0.084  for  the  X chromosome  and  h2  =  0.1 34  for  the 
third  chromosome),  whereas  coefficients  of  genetic  (CVG)  and  environmental  variance  (CVE)  were 
large,  as  is  characteristic  for  fitness  traits. 

Competitive  fitness  estimates  were  made  for  the  Chromosome  3  substitution  lines,  estimating 
viability  and  fertility,  against  a  competitor  strain  using  the  balancer  equilibrium  technique  described 
by  Sved  (1 975).  There  was  significant  variation  in  fitness  among  these  lines.  However,  no  significant 
correlations  were  observed  between  odor-guided  behavior  and  fitness  estimates.  This  suggests  that 
the  number  of  loci  causing  variation  in  olfactory  avoidance  behavior  may  be  small  relative  to  the 
number  of  loci  affecting  variation  in  fitness.  Alternatively,  the  controlled  laboratory  environment, 
which  does  not  require  flies  to  make  food  selection  or  oviposition  site  choices  or  to  avoid 
environmental  toxins,  may  not  be  suitable  for  evaluating  the  relationship  between  variation  in 
responsiveness  to  repellent  odorants  and  fitness. 

Intriguingly,  the  genetic  correlations  between  the  sexes  for  olfactory  avoidance  behavior  were 
extremely  low  (Fig.  2),  suggesting  that  different  genes  contribute  to  variation  in  avoidance  scores 
in  males  and  females  (Mackay  et  al. ,  1 996).  This  is  intuitively  easy  to  understand,  since  odor-guided 
behavior  subserves  different  functions  in  males  and  females,  e.g.  females  only  must  select  or  reject 
oviposition  sites.  The  lack  of  genetic  correlation  of  avoidance  scores  between  the  sexes  has, 
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however,  profound  evolutionary  implications.  Because  the  genetic  architecture  for  odor-guided 
behavior  appears  sexually  dimorphic,  it  is  clear  that  the  olfactory  subgenome  in  males  and  females 
may  evolve  along  different  evolutionary  trajectories,  resulting  in  sexual  dimorphism  for  olfactory 
behavior.  However,  the  low  genetic  correlation  between  the  sexes  (the  sex  environment  x  genotype 
interaction)  may  facilitate  the  maintenance  of  genetic  variation  for  olfactory  behavior.  Maintaining 
variation  in  the  trait  ensures  survival  of  at  least  some  members  of  the  species  under  diverse 
environmental  conditions  encountered  over  evolutionary  time.  At  the  same  time  selection  pressure 
is  maintained  at  each  generation,  since  no  single  genotype  can  satisfy  optimal  fitness  requirements 
for  both  sexes. 

3.  Transposon  tagging  and  the  identification  of  smell  impaired  genes 

Early  studies  on  olfaction  in  D.  melanogaster  led  to  the  identification  of  several  mutants, 
primarily  located  on  the  Jf-chromosome  (Rodrigues  and  Siddiqi,  1978;  Aceves-Pina  and  Quinn, 
1979;  Lilly  and  Carlson,  1989;  Helfand  and  Carlson,  1989;  McKenna  et  al.,  1989;  Ayer  and  Carlson, 
1992;  Woodard  et  al.,  1992;  Lilly  et  al.,  1994).  These  genes  were  discovered  using  either  larval 
assays,  in  which  mutant  larvae  were  identified  based  on  their  inability  to  locate  an  attractant  odor 
source  (Aceves-Pina  and  Quinn,  1979;  Rodrigues,  1980;  Monte  et  al.,  1989),  or  using  an  olfactory 
jump  assay,  in  which  flies  of  the  Canton-S  strain  failed  to  jump  in  response  to  benzaldehyde 
(Helfand  and  Carlson,  1989;  McKenna  et  al.,  1989).  One  of  these  acj  (abnormal  chemosensory 
jump)  mutations,  acj6,  encodes  a  POU-domain  transcription  factor,  which  controls  the  expression 
of  a  subpopulation  of  putative  odorant  receptors  (Clyne  et  al.,  1999b).  The  identification  of  acj6 
(Clyne  et  al.,  1999b),  several  genes  encoding  transduction  proteins  involved  in  chemosensory 
transduction  (Woodard  et  al,  1 992;  Riesgo-Escovar  et  al.,  1995;  Talluri  et  al. ,  1 995 ;  Stortkuhl  et  al. , 
1 999),  genes  encoding  odorant  binding  proteins  (Pikielny  et  al. ,  1 994;  McKenna  et  al. ,  1 994;  Kim 
et  al. ,  1 998),  the  smellblind  (sbl)  locus,  an  allele  of  paralytic  (para)  that  encodes  a  voltage  gated 
sodium  channel  (Lilly  and  Carlson,  1989;  Lilly  et  al. ,  1 994),  an  array  of  genes  that  mediate  olfactory 
learning  (reviewed  by  DeZazzo  and  Tully,  1995),  and,  finally,  members  of  at  least  one  family  of 
putative  odorant  receptors  (Clyne  et  al.,  1999a;  Vosshall  et  al.,  1999),  all  revealed  important 
components  that  participate  in  chemosensory  pathways  in  D.  melanogaster. 

A  comprehensive  understanding  of  the  genetic  architecture  of  odor-guided  behavior  will, 
ultimately,  require  the  identification  of  all  the  genes  that  contribute  to  this  trait  and  characterization 
of  their  interactions.  One  strategy  that  can,  in  principle,  accomplish  this  daunting  task  is  the  use  of 
P-element  insertional  mutagenesis,  which  enables  phenotypic  effects  to  be  linked  directly  to  gene 
expression  (Cooley  et  al. ,  1988;  Bellen  et  al. ,  1 989).  Introduction  of  a  transposon  in  the  genome  of 
D.  melanogaster  can  result  in  gene  disruption  at  or  near  the  site  of  insertion  of  the  transposable 
element.  Introduction  of  a  reporter  gene  in  genetically  engineered  transposable  element  constructs 
(e.g.  PjlArB ]),  which  can  be  driven  by  promoter/enhancer  elements  near  the  insertion  site  can  reveal 
expression  patterns  of  the  affected  gene  (enhancer  trap).  Furthermore,  the  use  of  a  cloning  vector, 
such  as  pBluescript,  in  the  construct  can  facilitate  cloning  of  flanking  sequences  adjacent  to  the  site 
of  L-element  insertion.  An  enhancer  trap  study  that  surveyed  6,400  lines  containing  different  P- 
element  insertion  sites  showed  that  about  45%  of  these  lines  displayed  expression  of  the  lacZ 
reporter  gene  in  the  third  antennal  segment  or  maxillary  palps,  and  several  of  these  lines  showed 
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specific,  non-uniform  staining  patterns  (Riesgo-Escovar  et  al.,  1992).  This  study,  however,  did  not 
report  whether  any  of  the  /’-element  insertions  affected  odor-guided  behavior. 

To  address  this  question,  we  introduced  the  P[lArB]  construct  at  random  locations  into  the 
autosomal  genome  of  a  highly  inbred  line  (Anholt  et  al. ,  1996).  The  co-isogenic  background  of  the 
resulting  P[lArB]  insertion  lines  together  with  the  statistical  behavioral  assay,  described  above, 
provided  sufficient  resolution  to  identify  not  only  mutations  of  large  effect,  but  also  smell 
impairments  of  smaller  effect  that  might  have  been  missed  with  conventional  mutant  screens.  We 
discovered  1 4  novel  olfactory  genes  by  screening  a  relatively  small  number  (379)  of  /-element  insert 
lines  (Anholt  et  al.,  1996).  The  frequency  with  which  smell  impaired  (smi)  genes  were  discovered 
was  about  4%,  suggesting  that  at  least  4%  of  the  genome  is  essential  in  mediating  the  simple 
avoidance  response  to  benzaldehyde.  /’-element  insertion  sites  can  be  readily  determined  through 
in  situ  hybridization  to  larval  salivary  gland  chromosomes  or  by  cloning  of  flanking  sequences  and 
comparing  these  sequences  to  the  genomic  data  base.  Candidate  genes  can  then  be  identified  for  each 
of  the  /-element  insertion  sites  (Table  1).  Conclusive  evidence  that  a  candidate  smi  gene  is  indeed 
the  locus  affected  by  the  transposon  and  that  reduction  in  the  expression  of  its  gene  product  is  causal 
to  observed  aberrations  in  the  olfactory  avoidance  response  requires  extensive  further 
experimentation.  First,  the  ./-element  can  be  mobilized  and  phenotypic  revertants  generated  to 
demonstrate  that  the  P [l ArBJ  insertion,  rather  than  an  unrelated  mutation,  is  responsible  for  the 
smell-impaired  phenotype.  Second,  olfactory  behavior  in  flies  carrying  different  alleles  in  the  same 
candidate  gene  must  be  evaluated  and  complementation  tests  between  the  original  PflArB ]  insertion 
line  and  flies  that  contain  deficiencies  or  other  ./-element  insertions  in  the  region  of  interest  must  be 
performed  to  provide  further  genetic  evidence  that  the  transposon-tagged  candidate  gene  is 
associated  with  the  impaired  phenotype.  Third,  the  expression  of  the  encoded  gene  product  must  be 
quantified,  either  by  Northern  blots,  quantitative  PCR,  or  Western  blots  (if  antibodies  are  available) 
to  show  that  the  smi  gene  message  or  SMI  protein  is  reduced  in  mutant  flies  as  compared  to  wild- 
type  or  phenotypic  revertants.  Fourth,  in  situ  hybridization  can  be  used  to  show  that  expression 
patterns  in  wild-type  flies  resemble  reporter  gene  expression  patterns  in  the  smi  mutant.  Finally, 
introduction  of  the  wild-type  gene  into  the  mutant  background  is  expected  to  rescue  the  mutant 
phenotype  if  the  candidate  gene  indeed  is  responsible  for  the  observed  deficiency  in  odor-guided 
behavior. 

4.  Epistasis  among  smi  genes 

The  availability  of  smell  impaired  lines  in  a  co-isogenic  genetic  background  provides  a 
unique  opportunity  to  investigate  epistasis  among  loci  that  all  affect  the  same  phenotype.  Small 
epistatic  effects  would  not  be  resolvable  in  the  presence  of  large  genetic  variance  in  diverse  genetic 
backgrounds  where  such  effects  would  be  confounded  by  segregation  of  unknown  modifiers  of  the 
phenotype.  However,  the  absence  of  genetic  variation  in  our  inbred  /-element  host  strain  and  the 
sensitivity  of  our  behavioral  assay  provided  enough  statistical  power  to  identify  enhancer  and 
suppressor  effects  among  smi  loci.  Double  heterozygous  hybrids  were  constructed  among  12 
independent  smi  mutations  in  a  classic  diallel  cross  design  (Griffing,  1 956)  and  avoidance  responses 
were  scored  for  these  transheterozygotes  (Fedorowicz  et  al.,  1998).  Heterozygous  effects  and 
epistasis  could  be  separated  by  determining  the  General  Combining  Ability  (GCA)  and  Specific 
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Combining  Ability  (SCA)  for  each  smi  mutant.  The  GCA  is  defined  as  the  average  avoidance  score 
of  a  smi  mutant  as  a  transheterozygote  in  combination  with  all  other  smi  mutations.  The  SCA  is 
defined  for  each  transheterozygous  genotype  as  the  difference  between  the  observed  avoidance  score 
of  each  specific  genotype  from  that  expected  from  the  sum  of  the  GCAs  of  each  smi  parent. 
Statistical  analyses  of  the  data  revealed  significant  epistatic  interactions  for  nine  transheterozygous 
genotypes  involving  10  of  the  12  smi  loci.  An  ensemble  of  eight  of  these  loci  could  be  represented 
as  an  epistatic  interaction  diagram  (Fig.  3;  Fedorowicz  et  al.,  1998). 

The  detection  of  epistasis  in  this  unique  set  of  smi  mutants  which  share  a  common  genetic 
background  vividly  illustrates  the  power  of  quantitative  genetic  analyses  to  detect  subtle  phenotypic 
effects.  These  observations  also  point  to  an  extensive  network  of  epistatic  interactions  among  genes 
in  the  olfactory  subgenome,  of  which  the  set  of  genes  analyzed  by  Fedorowicz  et  al.  (1998) 
represents  only  a  small  fraction.  Since  the  genes  that  constitute  the  interaction  diagram  shown  in  Fig. 
3  are  all  tagged  with  the  PflArBJ  transposon  they  can  be  characterized  at  the  molecular  level,  and 
epistasis  at  the  level  of  phenotype  can  then  be  evaluated  within  the  context  of  the  functions  of  their 
gene  products. 

5.  Understanding  the  genetic  architecture  of  behavior:  Linking  genetic  variation  to  phenotype 

One  of  the  best  understood  complex  traits  in  terms  of  its  underlying  genetic  architecture  is 
mechanosensory  bristle  number  in  Drosophila  (reviewed  by  Mackay,  1 996).  Although  a  vast  number 
of  genes  are  involved  in  determining  the  number  of  sternopleural  or  abdominal  bristles,  a  smaller 
number  accounts  for  most  of  the  variation  in  bristle  number  (Mackay,  1996).  These  genes  include 
many  neurogenic  genes,  such  as  Delta,  daughterless,  extramacrochaetae,  Hairless,  Enhancer-of- 
split,  hairy,  scabrous  and  the  achaete-scute  complex.  V ariation  in  bristle  number  could  be  correlated 
with  polymorphisms  in  Delta  (Long  et  al.,  1998;  Lyman  and  Mackay,  1998),  scabrous  (Lai  et  al., 
1 994;  Lyman  et  al. ,  1 999)  and  achaete-scute  (Mackay  and  Langley,  1 990;  Long  et  al. ,  2000)  and  the 
extent  to  which  alleles  of  each  of  these  loci  contribute  to  variation  in  bristle  number  in  natural 
populations  could  be  estimated. 

Similar  to  studies  on  the  contributions  of  bristle  genes  to  variation  in  bristle  number,  the 
extent  to  which  variation  in  smi  loci  contributes  to  naturally  occurring  variation  in  olfactory  behavior 
can  be  assessed  by  linkage  disequilibrium  mapping  (Lander  and  Schork,  1994;  Neimann-Sorensen 
and  Robertson,  1961;  Risch  and  Merikangas,  1996).  Linkage  disequilibrium  between  polymorphic 
loci  is  expected  to  diminish  over  time  as  recombination  occurs  between  them.  Over  a  long  period 
of  time,  only  very  closely  linked  polymorphisms  will  remain  in  linkage  disequilibrium.  Therefore, 
if  a  polymorphic  site  at  a  smi  gene  is  associated  with  differences  in  odor-guided  behavior,  that  site 
is,  or  is  closely  linked  to,  the  causal  site.  The  use  of  linkage  disequilibrium  mapping  requires  (1) 
identification  of  a  candidate  smi  gene  that  is  associated  with  naturally  occurring  variation  in  odor- 
guided  behavior;  (2)  identification  of  molecular  polymorphisms  in  that  smi  gene  that  are  segregating 
in  a  natural  population  (either  insertion/deletion  or  single  nucleotide  polymorphisms);  and,  (3) 
assessment  whether  alternative  molecular  variants  are  associated  with  variation  in  odor-guided 
behavior. 

The  genetic  architecture  of  behavior  is  likely  to  increase  in  complexity  when  compared  to 
that  of  morphological  traits.  Thus,  the  number  of  genes  that  contribute  prominently  to  variation  in 
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olfactory  avoidance  behavior  is  likely  to  be  far  greater  than  the  ensemble  that  contributes  to  variation 
in  bristle  number. 

The  first  layer  of  complexity  is  at  the  developmental  level  and  may  involve  several  of  the 
same  proneural  and  neurogenic  genes  that  regulate  bristle  number.  For  example,  lozenge  (Iz)  mutants 
fail  to  generate  basiconic  sensilla  (Stocker  and  Gendre,  1988).  Variation  in  the  number  of 
chemosensory  sensilla  among  individuals  in  natural  populations  has  not  been  investigated  and  could 
contribute  to  variation  in  olfactory  discrimination  and  responsiveness  to  odorants. 

The  next  level  of  complexity  is  the  perception  of  odorants,  and  involves  odorant  binding 
proteins  (Pikielny  et  al. ,  1 994;  McKenna  et  al.  ,1994;  Kim  et  al. ,  1 998)  and  odorant  receptors  (Clyne 
et  al.,  1999a;  Vosshall  et  al.,  1999).  Mutations  in  the  LUSH  odorant  binding  protein  resulted  in  a 
specific  olfactory  impairment,  responsiveness  to  ethanol,  propanol  and  butanol  (Kim  et  al.,  1998). 
Deficits  in  odorant  receptors  are  likely  to  be  odorant-specific  and  the  magnitude  of  smell  impairment 
will  depend  on  the  redundancy  of  recognition  of  the  odorant  by  other  receptors.  The  diversity  and 
relatively  small  size  of  the  Drosophila  odorant  receptor  family  suggests  that  this  system  has  less 
redundancy  in  odorant  recognition  than  its  mammalian  counterpart.  Nonetheless,  variations  in  the 
expression  of  single  odorant  receptors  may  contribute  less  phenotypic  variation  in  natural 
populations  than  polymorphisms  in  genes  that  encode  proteins  common  to  the  perception  of  all 
odorants.  The  best  documented  example  of  a  specific  olfactory  deficit  due  to  absence  of  an  odorant 
receptor  has  been  described  for  Caenorhabditis  elegans,  where  a  null  mutation  in  the  odr-10  gene, 
which  encodes  an  odorant  receptor  for  diacetyl,  results  in  impaired  chemotaxis  of  mutant  nematodes 
to  diacetyl  (Sengupta  et  al. ,  1 996).  To  what  extent  variation  in  the  expression  of  the  odorant  receptor 
repertoire  contributes  to  variation  in  olfactory  behavior  in  natural  populations,  however,  remains 
to  be  investigated. 

Genes  encoding  odorant  binding  proteins  and  odorant  receptors  confer  specificity  to  the 
behavior.  To  perform  the  behavior,  however,  several  additional  processes  must  occur:  signal 
transduction,  transfer  of  the  signal  from  the  periphery  to  the  central  nervous  system,  signal 
integration  and  the  generation  of  a  behavioral  response.  These  processes  are  mediated  via  gene 
products  that  are  not  unique  to  olfaction  and  phenotypes  resulting  from  defects  in  these  genes  are 
likely  to  be  pleiotropic.  For  example,  mutations  in  the  retinal  degeneration  B  ( rdgB )  gene,  which 
encodes  a  phosphatidyl  inositol  transfer  protein,  and  in  the  norpA  gene,  which  encodes  a 
phospholipase  C,  result  in  both  visual  and  olfactory  impairments  (Smith  et  al.,  1991;  Vihtelic  et  al. , 
1993;  Woodard  et  al,  1992;  Riesgo-Escovar  et  al.,  1994  and  1995). 

E-element  insertional  mutagenesis  can,  in  principle,  identify  genes  from  all  these  categories. 
Candidate  genes  for  previously  identified  smi  lines  are  listed  in  Table  1 .  Preliminary  characterization 
of  several  of  these  transposon  tagged  genes  implicate  genes  that  mediate  signal  propagation, 
including  a  voltage-gated  sodium  channel  (Kulkami,  Mackay  and  Anholt,  unpublished  observations) 
and  a  protein  containing  multiple  leucine  rich  repeats  and  PDZ  domains  likely  to  be  involved  in 
postsynaptic  organization  in  the  olfactory  pathway  (I.  Ganguly,  personal  communication).  Several 
novel  genes  of  unknown  function  have  also  been  implicated,  including  a  novel  tyrosine  regulated 
protein  kinase  (G.  Fedorowicz,  personal  communication).  Of  greatest  interest  is  our  recent 
realization  that  the  smi21F  insertion  tags  a  gene  encoding  a  novel  putative  pheromone  binding 
protein.  Previous  reporter  gene  expression  showed  differential  expression  in  the  proximal  dorsal 
region  of  the  antenna  as  well  as  specificity  in  smell  impairment  to  some,  but  not  all  odors  (Anholt 
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et  al.,  1996).  In  addition,  the  smi45E  insertion  appears  to  disrupt  the  function  of  a  nearby  gene 
encoding  an  olfactory  receptor.  Specificity  of  the  smell  impairment  suggests  that  this  receptor  has 
affinity  for  benzaldehyde.  The  characterization  of  these  candidate  genes  is  in  progress  and  will  be 
completed  in  the  near  future.  The  discovery  and  characterization  of  novel  gene  products  that  have 
major,  hitherto  unappreciated  effects  on  olfactory  behavior  will  provide  new  insights  in  the 
generation  and  regulation  of  odor-guided  behavior.  A  complete  understanding  of  the  genetic 
architecture  of  odor-guided  behavior,  ultimately,  requires  identification  of  all  the  genes  involved, 
characterization  of  their  interactions  in  shaping  the  phenotype  at  the  genetic  level,  and  identification 
of  polymorphisms  that  generate  phenotypic  variation  at  the  population  level.  Although  this  seems 
a  monumental  enterprise,  ongoing  improvements  in  genomic  technologies  are  bringing  the 
realization  of  this  goal  within  reach. 

(Hi)  Characterization  of  olfactomedin-related proteins 

Studies  on  the  molecular  evolution  of  olfactomedin-related  proteins  showed  that 
olfactomedin  homology  domains  are  highly  conserved  across  species  (Karavanich  and  Anholt, 
1998).  Olfactomedin  homologues  were  found  in  many  tissues,  including  the  brain  (Danielson  et  al. , 
1994),  and  in  species  ranging  from  Caenorhabditis  elegans  to  Homo  sapiens  (Karavanich  and 
Anholt,  1998).  These  observations  led  to  the  hypothesis  that  olfactomedin  homology  domains  may 
mediate  homophilic  protein-protein  interactions.  We  sought  to  provide  evidence  for  this  hypothesis 
in  the  eye,  where  the  glaucoma-associated  TIGR/myocilin  protein,  an  olfactomedin-related  protein, 
may  interact  with  other  olfactomedin-related  proteins.  Elucidating  such  interactions  would  also 
contribute  to  enhancing  our  understanding  of  the  pathogenesis  of  glaucoma. 

1.  Identification  and  characterization  of a  new  family  of olfactomedin-related proteins  in  the  human 
genome. 

The  prevalence  of  olfactomedin-related  proteins  among  species  and  their  identification  in 
different  tissues  prompted  us  to  investigate  whether  a  gene  family  exists  within  a  species,  specifically 
Homo  sapiens.  A  GenBank  search  indeed  revealed  an  entire  human  gene  family  of  olfactomedin- 
related  proteins  with  at  least  five  members,  designated  hOlfA  through  hOlfD,  and  the  TIGR/myocilin 
protein  (Kulkarni  et  al.,  2000).  hOlfA  corresponds  to  the  previously  described  rat  neuronal  AMZ 
protein  (Danielson  et  al.,  1994).  Phylogenetic  analyses  of  18  olfactomedin-related  sequences 
resolved  four  distinct  subfamilies  (Fig.  4).  Among  the  human  proteins,  hOlfA  and  hOlfC,  both 
expressed  in  brain,  are  most  closely  related.  Northern  blot  analyses  of  16  human  tissues 
demonstrated  highly  specific  expression  patterns:  hOlfA  is  expressed  in  brain,  hOlfB  in  pancreas  and 
prostate,  hOlfC  in  cerebellum,  hOlfD  in  colon,  small  intestine  and  prostate,  and  TIGR/myocilin  in 
heart  and  skeletal  muscle  (Fig.  5).  The  link  between  TIGR/myocilin  and  ocular  hypertension  and  the 
expression  of  several  of  these  proteins  in  mucus-lined  tissues  suggest  that  they  play  an  important  role 
in  regulating  physical  properties  of  the  extracellular  environment.  These  observations  led  to  the 
obvious  question  of  whether  any  of  these  olfactomedin-related  proteins  are  co-expressed  with 
TIGR/myocilin. 
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2.  Expression  of  multiple  olfactomedin-related proteins  by  trabecular  meshwork  cells. 


The  identification  of  additional  olfactomedin-related  proteins  enabled  us  to  design  primers 
for  PCR  amplification  and  to  assess  whether  messages  for  these  proteins  are  produced  in  trabecular 
meshwork  cells.  Indeed,  messages  that  encode  hOlfA  and  hOlfB,  but  not  hOlfC,  could  be  readily 
amplified  from  a  cDNA  library  from  trabecular  meshwork  cells.  Sequencing  of  the  amplification 
products  confirmed  that  they  indeed  represent  hOlfA  and  hOlfB.  Since  the  expression  of 
TIGR/myocilin  is  induced  by  glucocorticoids,  we  investigated  whether  glucocorticoids  would  also 
upregulate  expression  of  hOlfA  and  hOlfB.  Northern  blots  of  control  eyes  and  eyes  exposed  to 
dexamethasone  showed  dramatic  upregulation  of  TIGR/myocilin,  but  did  not  show  differences  in 
the  constitutive  expression  levels  of  hOlfA  and  hOlfB.  These  observations  suggest  a  possible  novel 
mode  of  action  for  TIGR/myocilin  in  the  pathogenesis  of  glaucoma.  In  this  model  hOlfA  and  hOlfB 
form  an  extracellular  matrix  in  the  intraocular  environment  similar  to  the  way  in  which  olfactomedin 
polymers  constitute  the  scaffold  of  the  lower  mucus  layer  that  covers  the  olfactory  neuroepithelium. 
Upregulation  of  TIGR/myocilin  would  disrupt  or  expand  this  matrix  through  interactions  between 
the  olfactomedin  domains  of  TIGR/myocilin  and  those  of  hOlfA  and  hOlfB.  It  is  this  alteration  in 
molecular  partners  between  different  intraocular  olfactomedins,  that  would  alter  the  flow 
characteristics  of  the  intraocular  fluid  and,  hence,  increase  the  intraocular  pressure.  The  first  step 
toward  verifying  that  such  molecular  interactions  actually  occur  is  to  assess  by  double 
immunofluorescence  whether  or  not  hOlfA,  hOlfB  and  TIGR/myocilin  co-localize  in  the  anterior 
segment  of  the  eye.  The  production  of  antibodies  against  these  three  proteins  is  in  progress. 

3.  Identification  of  olfactomedin  in  Drosophila  melanogaster 

We  used  a  signature  motif  of  the  olfactomedin  homology  domain  to  search  the  Drosophila 
genome  data  base  for  homologues.  In  contrast  to  Caenorhabditis  elegans,  which  has  two 
olfactomedin-related  proteins  (Karavanich  and  Anholt,  1998),  Drosophila  melanogaster  contains 
a  single  gene  that  expresses  a  protein  containing  an  olfactomedin  homology  domain.  This  is  the 
CG6867  gene  product  located  on  the  X-chromosome  at  cytological  location  16F1.  It  encodes  a 
protein  of  935  amino  acids  that  contains  a  coiled  coil  region,  a  collagen-like  domain,  two 
immunoglobulin-like  c2  domains,  and  a  large,  well-defined  olfactomedin  homology  domain  at  its 
C-terminus.  One  of  the  closest  mammalian  homologues  of  the  CG6867  gene  product  is  the 
TIGR/myocilin  protein.  The  CG6867  gene  product  is  clearly  an  extracellular  matrix  protein  designed 
to  engage  in  multiple  protein-protein  interactions,  a  feature  it  appears  to  share  with  all  other  members 
of  the  family  of  olfactomedin-related  proteins.  Identification  of  olfactomedin  in  Drosophila  is  likely 
to  facilitate  studies  on  its  function. 
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ABSTRACT 

Social  behaviors  of  most  mammals  are  profoundly  affected  by 
chemical  signals,  pheromones,  exchanged  between  conspeeifics.  Pher¬ 
omones  interact  with  dendritic  microvilli  of  bipolar  neurons  in  the 
vomeronasal  organ  (VNO).  To  investigate  vomeronasal  signal  trans¬ 
duction  pathways,  microvillar  membranes  from  porcine  VNO  were 
prepared.  Incubation  of  such  membranes  from  prepubertal  females 
with  boar  seminal  fluid  or  urine  results  in  an  increase  in  production 
of  inositol-(l,  4,  5)-trisphosphate  (IP3).  The  dose  response  for  IP3 
production  is  biphasic  with  a  GTP-dependent  component  at  low  stim¬ 
ulus  concentrations  and  a  nonspecific  increase  in  IP3  at  higher  stim¬ 
ulus  concentrations.  The  GTP-dependent  stimulation  is  mimicked  by 


GTPyS  and  blocked  by  GDPj3S.  Furthermore,  the  GTP-dependent 
component  of  the  stimulation  of  IP3  production  is  sex  specific  and 
tissue  dependent.  Studies  with  monospecific  antibodies  reveal  a 
Gttq/i! -related  protein  in  vomeronasal  neurons,  concentrated  at  their 
microvilli.  Our  observations  indicate  that  pheromones  in  boar  secre¬ 
tions  act  on  vomeronasal  neurons  in  the  female  VNO  via  a  receptor 
mediated,  G  protein-dependent  increase  in  IP3.  These  observations 
set  the  stage  for  further  investigations  on  the  regulation  of  stimulus- 
excitation  coupling  in  vomeronasal  neurons.  The  pheromone-induced 
IP3  response  also  provides  an  assay  for  future  purification  of  mam¬ 
malian  reproductive  pheromones.  ( Endocrinology  138:  3497-3504, 
1997) 


MOST  MAMMALS  use  pheromones  to  coordinate  re¬ 
production  (1).  These  chemical  signals  can  be  clas¬ 
sified  into  two  categories:  those  with  short-term  effects  on  the 
behavior  of  the  recipient  (signaling  pheromones),  and  those 
with  long  term  effects  on  the  physiology  of  the  recipient 
(priming  pheromones)  (1).  For  example,  signaling  phero¬ 
mones  in  urine  or  glandular  secretions  play  a  role  in  the 
initiation  of  copulatory  behavior  (2),  whereas  priming  pher¬ 
omones  are  responsible  for  puberty  acceleration  (3-8)  and 
reproductive  activation  (9,  10).  Although  the  majority  of 
studies  on  reproductive  pheromones  have  been  done  on 
rodents,  pheromone-dependent  effects  on  reproduction  have 
been  documented  also  for  sheep  (11),  cattle  (12),  and  pigs  (5, 
13,  14). 

One  of  the  most  extensively  studied  pheromonal  effects  is 
the  acceleration  of  puberty  in  the  female  house  mouse  (3,  4, 
6,  7).  The  presence  of  an  adult  male,  or  his  urine,  accelerates 
the  onset  of  puberty  in  prepubertal  female  mice  as  evident 
from  a  rapid  and  dramatic  increase  in  uterine  weight  (3,  4, 
6,  7).  The  induction  of  the  puberty  accelerating  pheromone 
is  androgen  dependent  because  urine  from  prepubertal 
males,  castrated  adult  males,  or  adult  females  fails  to  accel¬ 
erate  puberty  (4).  Pheromone-dependent  puberty  accelera¬ 
tion  is  not  unique  to  rodents  but  occurs  also  in  pigs  (5, 13), 
sheep  (11),  and  cows  (12). 

Whereas  the  perception  of  signaling  pheromones  may  be 
mediated  by  the  main  olfactory  system,  the  physiological 
effects  of  most  priming  pheromones  are  initiated  in  the 
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vomeronasal  organ  (VNO;  1,  15-17).  The  VNOs  are  paired, 
cartilage-encased  elongated  organs  associated  with  the 
vomer  bone  in  the  rostral  nasal  cavity.  The  VNO  contains  a 
lumen  that  communicates  via  a  duct  with  the  oral  (most 
mammals,  including  pigs)  or  nasal  (e.g.  horses)  cavity  (15, 
16).  Chemical  stimuli  in  urine  and  glandular  secretions  of 
conspeeifics  act  upon  the  dendritic  microvilli  of  bipolar  che- 
mosensory  neurons  in  the  VNO.  The  VNO  is  the  chemore- 
ceptive  organ  of  the  accessory  olfactory  system,  which  is 
functionally  and  anatomically  distinct  from  the  main  olfac¬ 
tory  system  (16, 18, 19).  The  main  olfactory  bulb  sends  pro¬ 
jections  to  the  primary  olfactory  cortex,  the  nucleus  of  the 
lateral  olfactory  tract,  the  olfactory  tubercle,  and  the  peri- 
amygdaloid  region,  while  afferent  neurons  from  the  VNO 
project  to  the  accessory  olfactory  bulb,  from  which  secondary 
neurons  extend  into  the  bed  nucleus  of  the  stria  terminalis, 
the  medial  amygdala,  and  the  hypothalamus,  enabling  pher¬ 
omones  to  influence  reproductive  physiology  and  behavior 
(16, 18).  Removal  of  the  VNO  or  the  accessory  olfactory  bulb 
impairs  reproductive  behavior,  whereas  lesions  of  the  main 
olfactory  epithelium  or  main  olfactory  bulb  that  result  in 
anosmia  do  not  impair  pheromonal  effects  (17, 19, 20).  Thus, 
whereas  the  main  olfactory  bulb  analyses  odors  in  the  en¬ 
vironment,  the  accessory  olfactory  bulb  is  specialized  to  de¬ 
tect  conspecific  chemical  signals  that  activate  stereotypic  in¬ 
stinctive  behaviors  via  the  neuroendocrine  system. 

In  the  main  olfactory  system,  odorants  bind  to  heptahe- 
lical,  G  protein-coupled  receptors  on  the  ciliated  dendrites  of 
olfactory  neurons  (21).  These  odorant  receptors  are  encoded 
by  a  diverse  array  of  about  1,000  genes  (21,  22).  Binding  of 
an  odorant  to  its  receptor  activates  a  heterotrimeric  G  protein 
(Gaolf;  23),  which  leads  to  activation  of  adenylyl  cyclase  (22, 
24).  The  resulting  increase  in  cAMP  elicits  the  generator 
potential  by  directly  opening  cyclic  nucleotide-gated  chan¬ 
nels  in  the  ciliary  plasma  membrane  (22,  24-26). 
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The  VNO  and  olfactory  epithelium  both  derive  from  the 
olfactory  placode.  Both  the  VNO  and  main  olfactory  ep¬ 
ithelium  possess  bipolar  neurons  that  are  functionally  re¬ 
placed  from  neurogenic  precursor  cells  throughout  life 
(27).  In  addition,  in  both  systems  primary  chemosensory 
neurons  form  convergent  projections  onto  the  main  olfac¬ 
tory  bulb  or  accessory  olfactory  bulb.  Similarities  in  em¬ 
bryonic  development,  anatomical  organization,  and  che¬ 
mosensory  function  initially  suggested  that  transduction 
mechanisms  in  the  VNO  would  resemble  those  in  the 
olfactory  epithelium.  However,  screens  of  complementary 
DNA  (cDNA)  libraries  from  murine  VNO  failed  to  yield 
cDNAs  encoding  Gaolf,  adenylate  cyclase  type  III,  and  the 
a-subunit  of  the  olfactory  cyclic  nucleotide-gated  channel 
(28,  29).  Studies  in  rat  also  indicate  that  the  olfactory 
epithelium  and  VNO  differ  in  signal  transduction  path¬ 
ways  (30).  Electrophysiological  characterization  of  che¬ 
mosensory  neurons  from  the  murine  VNO  failed  to  detect 
functional  cyclic  nucleotide-gated  channels  in  vomerona¬ 
sal  neurons  (31,  32).  Finally,  a  family  of  putative  phero¬ 
mone  receptors  has  been  identified  in  the  VNO  (33).  Al¬ 
though  these  belong  to  the  superfamily  of  heptahelical 
receptors,  they  do  not  share  motifs  characteristic  of  the 
family  of  odorant  receptors  (21).  These  studies  all  support 
the  notion  that  chemosensory  neurons  of  the  main  olfac¬ 
tory  system  and  the  accessory  olfactory  system  use  dif¬ 
ferent  signal  recognition  and  transduction  pathways  (32). 

Functional  biochemical  studies  on  the  VNO  that  would 
complement  molecular  biological  approaches  have  been 
hampered  by  the  small  size  of  the  VNO  in  most  common 
laboratory  animals.  To  eliminate  this  problem,  we  selected 
the  domestic  pig  ( Sus  scrofa)  as  our  model  system.  Pigs 
have  large  VNOs  that  are  well  separated  from  the  main 
olfactory  system,  and  their  reproductive  physiology  and 
behavior,  like  that  of  rodents,  is  regulated  by  pheromones 
(5,  13,  14).  We  have  developed  a  procedure  for  the  prep¬ 
aration  of  a  VNO  membrane  fraction  enriched  in  dendritic 
microvillar  membranes  that  allows  routine  measurements 
of  second  messengers.  Here  we  report  that  pheromones 
contained  in  seminal  fluid  and  boar  urine  stimulate  the 
production  of  inositol-(l,  4,  5)-trisphosphate  (IP3)  when 
applied  to  microvillar  VNO  membranes  from  female  ju¬ 
venile  pigs  (gilts).  This  stimulation  is  dose  dependent,  GTP 
dependent,  sex  dependent,  and  tissue  specific.  We  show 
further  that  this  microvillar  membrane  preparation  con¬ 
tains  a  G  protein  of  the  Gaq  class,  which  commonly  me¬ 
diates  receptor-activated  increases  in  IP3  (34).  Finally,  im- 
munohistochemical  studies  show  prominent  Gnq/11- 
immunoreactivity  concentrated  at  the  microvilli  of 
vomeronasal  neurons,  supporting  the  notion  that  phero¬ 
mones  in  seminal  fluid  and  boar  urine  stimulate  vomer¬ 
onasal  neurons  of  female  pigs  via  a  receptor  mediated 
Gttq/n  coupled  IP3  pathway. 

Materials  and  Methods 

Membrane  preparations 

Freshly  collected  boar  urine,  gilt  urine,  and  seminal  fluid  were  pro¬ 
vided  by  Dr.  W.  L.  Flowers  from  the  Animal  Science  Department  at 
North  Carolina  State  University.  The  seminal  fluid  was  centrifuged  to 
remove  cells.  The  seminal  fluid,  boar  urine  and  gilt  urine  were  stored  as 


aliquots  under  argon  at  -80  C  until  used.  Pigs  were  made  available 
immediately  following  euthanasia  by  Drs.  R.  A.  Argenzio,  N.  A.  Mon- 
teiro-Riviere,  and  R.  A  Abdullahi,  from  the  College  of  Veterinary  Med¬ 
icine  at  North  Carolina  State  University.  VNOs  from  gilts,  up  to  six 
months  old,  were  dissected  from  their  crevices  in  the  nasal  cavity, 
removed  from  the  cartilaginous  capsule,  and  frozen  on  dry  ice.  The 
tissues  were  then  minced  and  crushed  with  a  razor  blade  and  subjected 
to  sonication  for  2-5  min  in  ice-cold  PBS  in  a  Bransonic  bath  sonicator. 
The  resulting  suspension  was  layered  on  a  45%  (wt/ wt)  sucrose  cushion 
and  centrifuged  at  4  C  for  30  min  at  40,000  rpm  in  a  Beckman  SW55Ti 
rotor.  The  membrane  fraction  on  top  of  the  sucrose  was  collected  and 
centrifuged  as  before  for  15  min  to  pellet  the  membranes.  The  mem¬ 
branes  were  resuspended  in  100  pi  of  ice-cold  PBS.  Protein  was  then 
determined  according  to  the  method  of  Lowry  et  ah  (35),  using  BSA  as 
standard.  Membranes  from  olfactory  tissue  were  prepared  according  to 
the  same  procedure.  Membranes  from  liver,  brain,  kidney,  and  lung 
were  prepared  by  homogenizing  the  tissue  in  PBS  with  a  Teflon  ho- 
mogenizer.  Membranes  were  collected  by  centrifugation,  washed  once, 
and  resuspended  in  PBS. 

Second  messenger  assays 

Forskolin  and  nucleotides  were  purchased  from  Boehringer  Mann¬ 
heim  (Indianapolis,  IN),  [3H]  cAMP  and  [32P]-a-ATP  were  from  Am- 
ersham  Radiochemical  Corporation  (Arlington  Heights,  IL).  Adenylate 
cyclase  activity  was  measured  according  to  the  method  of  Salomon  et  ah 
(36),  in  the  presence  of  10  pM  forskolin,  boar  urine  or  10  pM  guanosine 
5'-0-(3-thiotriphosphate  (GTP7S).  For  IP3  assays,  reactions  were  incu¬ 
bated  for  1  min  at  37  C  in  25  mM  Tris-acetate  buffer,  pH  7.2,  5  mM 
Mg-acetate,  1  mM  dithiothreitol,  0.5  mM  ATP,  0.1  mM  CaCl2, 0.1  mg /ml 
BSA,  10  (am  GTP,  and  20  pg  VNO  membrane  protein.  Reactions  were 
terminated  by  adding  1  m  trichloroacetic  acid.  IP3  was  measured  with 
a  kit  from  New  England  Nuclear,  Inc.  (Boston,  MA)  according  to  the 
manufacturer's  instructions  and  is  based  on  displacement  of  [3H]  IP3 
from  a  specific  IP3  binding  protein.  Differences  between  experimental 
and  control  animals  were  analyzed  using  the  Student's  t  test. 

Western  Blotting 

VNO  membrane  samples  were  subjected  to  electrophoresis  on  a  10% 
SDS-polyacrylamide  gel,  followed  by  electrophoretic  transfer  onto  a 
nitrocellulose  membrane.  Strips  of  the  membrane,  containing  approxi¬ 
mately  10  pg  protein  were  probed  with  a  1,000-fold  dilution  of  normal 
rabbit  serum  or  1,000-fold  dilutions  of  rabbit  antisera  against  specific  G 
protein  subunits  (Calbiochem,  La  Jolla,  CA).  Bound  antibody  was  vi¬ 
sualized  via  a  biotinylated  goat-antirabbit  secondary  antibody  com- 
plexed  with  avidin  and  biotinylated  horseradish  peroxidase,  using  Am- 
ersham's  chemiluminescent  ECL  detection  system.  Migration  distances 
were  calibrated  with  biotinylated  low  range  molecular  weight  markers 
(Bio-Rad,  Richmond,  CA). 

Immunohistochemistry 

Formalin-fixed  and  paraffin-embedded,  5-pm  thick,  coronal  sections 
through  the  VNO  were  deparaffinized  in  xylene  and  rehydrated  through 
graded  alcohols.  Sections  were  blocked  with  0.5%  casein  for  1  h  at  room 
temperature.  They  were  then  incubated  with  a  200-fold  dilution  of  either 
normal  rabbit  serum  or  anti-Gaq/11  antiserum  in  PBS,  0.05%  Triton 
X-100.  This  was  followed  by  several  washes  and  incubation  for  1  h  at 
room  temperature  with  biotinylated  goat  antirabbit  IgG  in  PBS  supple¬ 
mented  with  0.05%  Triton  X-100  and  1%  normal  rabbit  serum.  At  this 
stage,  endogenous  peroxidase  activity  was  abolished  by  exposing  the 
sections  for  10  min  to  0.3%  hydrogen  peroxide.  After  addition  of  strepta- 
vidin-biotinylated  horseradish  peroxidase  complex  (Zymed  Laborato¬ 
ries,  Inc.,  San  Francisco,  CA),  bound  antibody  was  visualized  as  brown 
immunoprecipitates  using  3,3  'diaminobenzidine  as  chromogenic  sub¬ 
strate.  Sections  were  then  counterstained  with  hematoxylin,  and  viewed 
and  photographed  under  a  Zeiss  Axioplan  microscope. 
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Results 

Dose-dependent  and  GTP-dependent  increases  in  IP3  levels 
induced  by  boar  seminal  fluid  and  urine  in  VNO 
membranes  from  gilts 

To  study  transduction  pathways  activated  by  pheromonal 
stimuli  from  the  male,  we  developed  a  preparation  enriched 
in  microvillar  membranes  from  VNOs  of  gilts.  The  tissue  is 
subjected  to  sonication  to  detach  microvilli  from  the  vome¬ 
ronasal  neuroepithelial  surface,  and  the  resulting  mem¬ 
branes  are  then  collected  by  centrifugation  on  a  sucrose  cush¬ 
ion.  The  yield  of  these  membranes  is  141  ±  9  jug  protein/ 
VNO  (n  =  11)  and  this  membrane  fraction  is  approximately 
3-fold,  enriched  in  both  the  specific  activities  of  adenylate 
cyclase  and  phospholipase  C.  The  basal  activity  of  adenylate 
cyclase  (68  ±  12  pmol/min-mg  n  =  5)  can  be  readily  activated 
by  the  nonhydrolyzable  GTP  analogue,  GTPyS,  and  forskolin 
(117  ±  13  and  228  ±  38  pmol/min-mg,  respectively,  n  =  5), 
but  we  could  not  detect  stimulation  with  pheromonal  stimuli 
from  either  seminal  fluid,  boar  urine,  or  the  known  boar 
pheromone,  5a-androst-16-en-3-one  (Sigma  Chemical  Com¬ 
pany,  St.  Louis,  MO).  Vomeronasal  adenylate  cyclase  activity 
is  approximately  50-fold  lower  than  activities  observed  in 
olfactory  cilia  preparations  (37,  38)  and  is  not  affected  by 
calcium/calmodulin  (38).  This  is  in  line  with  observations  by 
Berghard  and  Buck  which  indicate  that  the  vomeronasal 
adenylate  cyclase  is  the  calmodulin-insensitive  type  II  iso¬ 
form  rather  than  the  adenylate  cyclase  type  III  found  in 
olfactory  receptor  cells  (29). 

In  contrast  to  the  vomeronasal  adenylate  cyclase,  incuba¬ 
tion  of  microvillar  membranes  from  gilts  with  seminal  fluid 
results  in  a  robust,  dose-dependent  increase  in  IP3  (Fig.  1). 
The  dose-response  curve  is  biphasic.  At  lower  stimulus  con- 
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Fig.  1.  Dose  dependence  of  the  production  of  IP3  by  seminal  fluid  in 
female  VNO  membranes.  The  dose-response  curve  for  activation  of 
IP3  is  biphasic  with  a  specific  GTP-dependent  component  at  lower 
concentrations  of  seminal  fluid  and  a  nonspecific  GTP-independent 
component  at  higher  concentrations.  Each  data  point  represents  the 
mean  of  four  to  six  independent  experiments,  each  consisting  of  du¬ 
plicate  measurements.  SEs  are  within  18%  of  the  mean  for  all  mea¬ 
surements. 


centrations  a  component  of  the  response  saturates  at  the  level 
of  stimulation  observed  with  GTPyS.  At  higher  stimulus, 
concentrations  an  additional  nonsaturable,  nonspecific  in¬ 
crease  in  IP3  is  observed.  The  response  seen  at  low  concen¬ 
trations  of  seminal  fluid  (up  to  1 .5%  vol/ vol)  is  mimicked  by 
GTPyS  and  blocked  by  GDP/3S  (Fig.  2).  Similar  responses  are 
elicited  with  urine  from  the  boar,  but  here  higher  stimulus 
concentrations  (up  to  3%  vol /vol)  are  required  to  resolve  the 
GTP-dependent  component  of  the  IP3  response.  No  increases 
in  IP3  were  observed  in  response  to  5a-androst-16-en-3-one 
up  to  concentrations  of  100  mM.  This  is  consistent  with  pre¬ 
vious  studies  that  suggest  that  androstenone  is  a  signaling 
pheromone  which  mediates  its  effects  via  the  main  olfactory 
rather  than  the  accessory  olfactory  system  (39).  We  conclude 
from  these  observations  that  female  VNO  membranes  re¬ 
spond  to  stimuli  in  boar  seminal  fluid  and  urine  with  an 
increase  in  IP3  via  a  G  protein  coupled  pathway. 

Sex  specificity  and  tissue  specificity  of  the  IP3  response 

Relatively  little  information  is  available  about  the  chemical 
nature  of  mammalian  pheromones.  Because  the  active  com¬ 
ponents  in  boar  seminal  fluid  and  urine  have  not  been  iden¬ 
tified,  we  deemed  it  important  to  establish  whether  the  ob¬ 
served  IP3  responses  were  physiologically  relevant.  Previous 
studies  on  rodents  have  shown  that  production  of  repro¬ 
ductive  pheromones  by  the  male  is  androgen  dependent  (4). 
We  reasoned,  therefore,  that  only  boar  urine,  but  not  urine 
from  prepubertal  females,  should  elicit  an  increase  in  IP3  in 
VNO  membranes  from  gilts.  The  results  presented  in  Fig.  3 
demonstrate  that  this  prediction  is  correct.  Whereas  boar 
urine  elicits  a  robust  increase  in  IP3,  experiments  using  urine 
samples  from  gilts  failed  to  activate  IP3  production  above  the 
basal  level.  To  further  document  specificity  of  the  observed 
IP3  response,  we  investigated  the  tissue  specificity  of  seminal 
fluid-induced  increases  in  IP3  by  testing  the  effect  of  1.5% 
seminal  fluid  side-by-side  on  microvillar  VNO  membranes 
and  membranes  obtained  from  olfactory  tissue,  brain,  lung, 
liver,  and  kidney.  Basal  activity  of  phospholipase  C  was 
detected  in  all  samples  and  was  particularly  high  in  mem¬ 
branes  from  olfactory  tissue  (Fig.  4).  However,  whereas  sem¬ 
inal  fluid  caused  a  2-fold  increase  in  IP3  production  in  VNO 
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Fig.  2.  GTP  dependence  of  the  production  of  IP3  by  seminal  fluid  in 
female  VNO  membranes.  Reactions  were  performed  without  stimu¬ 
lus,  in  the  presence  of  1.5%  seminal  fluid  (vol/vol),  10  pM  GTPyS,  or 
1.5%  seminal  fluid  together  with  100  pM  GDP/3S.  Significant  stim¬ 
ulation  compared  with  basal  activity  is  observed  in  the  presence  of 
seminal  fluid  and  GTPyS  (*,  P  <  0.05)  and  not  in  the  presence  of 
seminal  fluid  together  with  GDP/3S. 
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Fig.  3.  Sex  dependence  of  the  produc¬ 
tion  of  IP3  in  female  VNO  membranes. 
Reactions  were  performed  without 
stimulus,  in  the  presence  of  3%  (vol/vol) 
gilt  urine,  3%  (vol/vol)  boar  urine,  1.5% 
(vol/vol)  seminal  fluid,  and  10  /am 
GTPyS.  All  assays  were  done  in  the 
presence  of  10  fxM  GTP.  Significant 
stimulation  compared  with  basal  activ¬ 
ity  is  observed  in  the  presence  of  boar 
urine,  seminal  fluid,  and  GTPyS  (*,  P  < 
0.05).  Levels  of  IP3  production  by  boar 
urine,  seminal  fluid,  and  GTPyS  are  not 
statistically  different.  There  is  also  no 
significant  difference  between  the  base¬ 
line  control  and  stimulation  by  gilt 
urine. 
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Fig.  4.  Tissue  specificity  of  the  produc¬ 
tion  of  IP3  by  seminal  fluid  in  female 
VNO  membranes.  Reactions  were  per¬ 
formed  using  20  fxg  of  membrane  pro¬ 
tein  and  1.5%  seminal  fluid  (vol/vol)  as 
the  stimulus.  There  were  no  significant 
differences  between  the  basal  level  of 
IP3  ( open  bars)  and  its  level  in  the  pres¬ 
ence  of  seminal  fluid  ( closed  bars )  with 
membranes  from  olfactory  tissue 
(OLF),  kidney,  liver,  or  lung.  Signifi¬ 
cant  differences  between  control  and 
stimulated  IP3  levels  are  observed  in 
VNO  membranes  and  brain  mem¬ 
branes  (*,  P  <  0.05). 
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membranes,  no  significant  stimulation  above  the  basal  level 
was  observed  in  membranes  from  olfactory  tissue,  liver, 
lung,  and  kidney.  Besides  the  VNO,  significant  increases  in 
the  level  of  IP3  upon  exposure  to  seminal  fluid  were  observed 
only  in  membranes  from  brain.  We  attribute  this  stimulation 
to  the  presence  of  neuroactive  substances  in  seminal  fluid. 

The  experiments  described  above  demonstrate  a  GTP-de- 
pendent  and  dose-dependent  increase  in  the  formation  of  IP3, 
but  not  cAMP,  in  female  VNO  membranes  upon  exposure  to 
boar  seminal  fluid  or  urine.  This  stimulation  shows  tissue 


specificity  and  sex  dependence.  Together,  these  observations 
indicate  that  pheromones  from  the  male  activate  the  pro¬ 
duction  of  IP3  in  the  female  VNO  via  specific  G  protein- 
coupled  receptors. 

Identification  of  a  Gaq/11  related  G  protein  on  the 
microvillar  surface  of  the  VNO 

Previously,  immunohistochemical  studies  identified 
a-subunits  of  Gi2  and  G0  in  distinct  subpopulations  of  vome- 
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ronasal  neurons.  Although  we  cannot  exclude  that  either  of 
these  two  G  proteins  could  mediate  the  observed  stimulation 
of  phospholipase  C,  we  decided  to  investigate  whether  our 
VNO  membrane  preparation  contains  a  G  protein  of  the  Gaq 
class,  known  to  cause  activation  of  at  least  the  J31 -isoform  of 
phospholipase  C  (34).  Using  subunit-specific  rabbit  antisera 
against  unique  peptide  sequences  of  Gas,  Gan/Gaa,  Gan/ 
Gai3/Ga0,  Ga  i3,  Gaq/11/  and  Gp,  we  confirmed  in  our  mi- 
crovillar  membrane  preparation  the  presence  of  Gan  and 
Ga0,  as  reported  previously  (29,  40;  Fig.  5).  In  addition,  we 
observe  Gas,  and  a  single  /3-subunit  at  35  kDa.  A  Gai3/GaQ 
specific  antiserum  sometimes  reveals  a  doublet  of  immuno- 
reactive  bands,  suggesting  the  presence  of  both  Gai3  and  Gac. 
Of  greatest  interest,  however,  is  the  observation  of  a  prom¬ 
inent,  previously  unreported,  immunoreactivity  revealed  by 
an  antiserum  against  the  Gaq/11  protein  (Fig.  5). 

We  investigated  whether  the  Gaq/11  protein  is  localized  to 
the  microvillar  surface  of  the  VNO.  The  apical  region  of 
porcine  vomeronasal  neurons  contains  a  dense  group  of  mi¬ 
crovilli  in  contrast  to  supporting  cells  that  carry  smaller 
groups  of  microvilli  (41).  Immunohistochemical  staining  of 
coronal  sections  through  the  VNO  with  the  Gaq/11  antiserum 
reveals  intense  staining  of  microvillar  tufts  at  the  surface  of 
the  vomeronasal  lumen  (Fig.  6A).  Cell  bodies  of  vomeronasal 
neurons  and  dendritic  processes  are  also  stained,  but  the 
staining  here  is  lighter,  and  it  is  obvious  that  immunoreac¬ 
tivity  is  concentrated  primarily  at  the  microvillar  surface.  To 
verify  the  specificity  of  this  staining,  adjacent  sections  were 
incubated  either  without  the  primary  antibody  or  with  nor¬ 
mal  rabbit  serum  at  the  same  concentration.  Under  these 
conditions,  no  staining  was  detected  (Fig.  6B).  Thus,  Gaq/11, 
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antisera  against  G-protein  subunits 


Fig.  5.  Identification  of  G  protein  subunits  in  vomeronasal  mem¬ 
branes.  Each  strip  contained  10  pg  of  VNO  membrane  protein  and 
was  probed  either  with  1,000-fold  dilution  of  normal  rabbit  serum 
(NRS)  or  monospecific  rabbit  antisera  against  subunits  of  G  proteins, 
as  indicated.  The  arrow  indicates  the  45-kDa  Gas  subunit,  identified 
by  the  Ga8-antiserum.  A  single  /3-subunit  is  detected  which  migrates 
at  35  kDa.  The  immunoreactive  bands  identified  by  antisera  against 
Gai3  and  Gui3/Gai0  can  sometimes  be  resolved  as  doublets,  suggesting 
the  presence  of  both  Gai3  and  Ga0.  Note  also  the  presence  of  Gai2  and 
the  prominent  presence  of  a  polypeptide  immunoreactive  with  anti- 
Gaq/n. 


Fig.  6.  Immunohistochemical  localization  ofGaq/11  to  the  microvillar 
surface  of  the  VNO.  A,  Section  stained  with  a  200-fold  dilution  of 
antiserum  against  Gaq/11.  Note  the  deposition  of  brown  reaction  prod¬ 
uct  on  microvillar  tufts  along  the  microvillar  surface  of  the  neuroep¬ 
ithelium  ( arrows )  and  in  the  cell  bodies  of  the  vomeronasal  neurons 
0 arrowheads ).  B,  Adjacent  section  incubated  with  normal  rabbit  se¬ 
rum.  Scale  bar,  100  pm. 

which  links  receptor  mediated  responses  to  activation  of 
phospholipase  C,  is  expressed  by  vomeronasal  neurons  and 
concentrated  at  the  microvillar  surface,  where  the  phero¬ 
mone  first  encounters  the  chemosensory  neurons. 

Discussion 

Our  results  suggest  that  in  the  porcine  VNO  IP3  is  the 
principal  second  messenger  that  mediates  pheromonal  sig¬ 
nal  transduction.  This  is  in  contrast  to  the  main  olfactory 
system,  where  cAMP  appears  to  be  the  dominant  second 
messenger  that  regulates  olfactory  transduction  (22,  24-26). 
Although  we  cannot  exclude  a  role  for  cAMP  in  vomeronasal 
signal  transduction,  we  were  unable  to  detect  effects  of  either 
boar  urine  or  boar  seminal  fluid  on  vomeronasal  adenylate 
cyclase.  Failure  to  detect  message  for  Gaolf,  adenylate  cyclase 
type  III  and  the  a-subunit  of  the  olfactory  cyclic  nucleotide 
gated  channel  in  cDNA  libraries  from  VNO  (29,  30)  further 
accentuates  the  differences  between  chemosensory  transduc¬ 
tion  in  olfactory  and  vomeronasal  neurons. 

Involvement  of  IP3  has  also  been  implicated  in  signal  trans¬ 
duction  in  the  reptilian  VNO.  In  garter  snakes,  a  polypeptide 
purified  from  secretions  of  earthworms,  the  snake's  prey. 
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induces  a  G  protein-dependent  increase  in  IP3  and  a  decrease 
in  c AMP,  suggesting  interactions  between  both  signal  trans¬ 
duction  pathways  (42, 43).  Furthermore,  patch-clamp  studies 
on  vomeronasal  neurons  from  the  turtle  have  shown  that 
intracellular  injection  of  IP3  elicits  a  membrane  conductance 
(44),  although  effects  of  cAMP  have  also  been  reported  in  this 
system  (45).  Thus,  the  role  of  IP3  in  vomeronasal  signal  trans¬ 
duction  may  be  universal  among  vertebrates,  whereas  the 
influence  of  cAMP  in  vomeronasal  signaling  remains  to  be 
evaluated  further. 

Despite  behavioral,  electrophysiological,  and  molecular 
biological  studies  on  the  mammalian  VNO,  no  direct  mea¬ 
surements  of  dose-dependent,  GTP-dependent  increases 
in  second  messengers  upon  exposure  to  pheromones  have 
been  reported  until  now.  Evaluation  of  the  dose  depen¬ 
dence  of  vomeronasal  stimulation  is  especially  important 
in  light  of  the  large  nonspecific  component  we  observed. 
In  light  of  the  previously  observed  diversity  of  putative 
pheromone  receptors  (33),  it  may  seem  puzzling  that  in 
our  case  specific  G  protein-dependent  production  of  IP3  in 
response  to  boar  seminal  fluid  approaches  the  level  of  IP3 
production  observed  with  GTPyS.  However,  this  obser¬ 
vation  can  be  explained,  if  we  assume  that  boar  seminal 
fluid  contains  a  cocktail  of  pheromones  and  that  the  high 
stimulus  concentrations  used  in  our  study  activate  the  vast 
majority  of  pheromone  receptors,  coupled  to  the  IP3  path¬ 
way.  Although  the  chemical  nature  of  the  pheromonal 
stimuli  in  boar  urine  and  boar  seminal  fluid  are  not  known, 
we  believe  that  the  IP3  responses  we  report  reflect  the 
physiological  response  of  the  vomeronasal  neuron  to  pher¬ 
omones  because:  1)  the  response  shows  a  saturable  dose- 
dependent  component;  2)  this  component  is  GTP  depen¬ 
dent,  indicating  the  involvement  of  a  G  protein,  and  by 
inference,  G  protein-coupled  receptors  (33);  3)  the  re¬ 
sponse  is  sex  dependent,  as  expected  from  pheromones 
produced  only  by  the  male  and  intended  to  affect  the 
female;  4)  the  response  shows  tissue-specificity,  in  agree¬ 
ment  with  the  notion  that  it  is  mediated  via  receptors 
expressed  selectively  in  the  VNO;  and  5)  a  Gaq/11  type  G 
protein,  classically  involved  with  receptor-mediated  acti¬ 
vation  of  phospholipase  C  (34),  is  expressed  in  vomero¬ 
nasal  neurons  and  concentrated  on  their  chemosensory 
microvilli. 

The  procedure  used  for  the  preparation  of  microvillar 
membranes  is  modeled  after  well  established  methods  for 
harvesting  olfactory  cilia  from  olfactory  neuroepithelium 
(46,  47).  Sonication  of  olfactory  membranes  results  not  only 
in  the  detachment  of  olfactory  cilia  but  also  in  the  detachment 
of  microvilli  from  sustentacular  cells  and  plasma  membrane 
fragments  from  other  components  of  the  neuroepithelium 
(46,  47).  Electron  microscopic  examination  of  these  prepara¬ 
tions  revealed  membrane  vesicles,  axonemal  structures  de¬ 
void  of  a  plasma  membrane,  and  axonemal  structures  asso¬ 
ciated  with  membrane  fragments  (47).  The  membrane 
preparation  we  refer  to  as  "microvillar  membranes"  is,  there¬ 
fore,  likely  to  contain  contaminants  derived  from  other  com¬ 
ponents  of  the  VNO,  including  microvillar  membranes  from 
supporting  cells,  and  it  is  difficult  to  estimate  the  purity  of 
this  preparation  precisely.  However,  our  preparation  ap¬ 
pears  to  be  sufficiently  enriched  in  chemosensory  mem¬ 


branes  for  the  purpose  of  our  studies.  This  assessment  is 
based  on  the  fact  that  the  observed  pheromone  induced 
responses  are  tissue  specific,  sex  dependent,  and  G  protein 
mediated.  It  is  further  supported  by  the  notion  that  the  prom¬ 
inent  expression  of  Gaq/n  immunoreactivity  at  the  microvil¬ 
lar  surface  of  the  neuroepithelium  (Fig.  6)  mirrors  the  prom¬ 
inent  visualization  of  Gaq/11  immunoreactivity  observed  in 
the  microvillar  membrane  preparation  on  Western  blots 
(Fig.  5). 

Immunohistochemical  studies  on  opossum  (40)  and  in 
situ  hybridization  studies  in  mouse  (29)  revealed  that  a 
population  of  neurons  in  the  apical  layer  of  the  vomero¬ 
nasal  neuroepithelium  expresses  Gai2.  These  neurons 
project  to  the  anterior  region  of  the  accessory  olfactory 
bulb  (40).  In  contrast,  a  population  of  neurons  located  in 
the  base  of  the  vomeronasal  neuroepithelium  expresses 
Ga0  (29)  and  projects  to  the  posterior  region  of  the  acces¬ 
sory  olfactory  bulb  (40).  Both  a-subunits  of  these  G  pro¬ 
teins  are  also  detected  in  the  VNO  of  the  pig.  It  remains 
to  be  determined  whether  these  G  proteins  play  roles 
directly  in  pheromonal  transduction  or  in  signaling  pro¬ 
cesses  that  coordinate  the  growth  and  differentiation  of 
vomeronasal  neurons.  In  an  extensive  screen  of  a  rat  VNO 
cDNA  library,  Berghard  and  Buck  (29)  detected  several 
cDNAs  encoding  Ga^.  It  seems,  therefore,  reasonable  to 
presume  that  the  Gaq/11  immunoreactivity  we  detect  most 
likely  represents  Gan.  Because  the  frequency  of  clones 
encoding  Galx  in  the  library  screened  by  Berghard  and 
Buck  (29)  was  low  relative  to  cDNAs  encoding  Gai2  and 
Ga0/  it  appears  that  low  levels  of  message  are  produced  for 
the  Ga1:l  protein,  which  may  reflect  a  slower  turnover  than 
Gai2  and  Ga0.  Although  a  direct  link  between  pheromone 
detection  and  activation  of  Gan  must  still  be  documented, 
the  uniform  presence  of  this  G  protein  on  all  microvillar 
tufts  suggests  a  role  for  this  G  protein  in  pheromonal 
signaling  in  all  mature  vomeronasal  neurons.  Localization 
of  Gaq/11  to  the  neuronal  compartment  of  the  VNO  is 
supported  by  the  observation  that  neuronal  cell  bodies  of 
the  vomeronasal  epithelium  also  stain  and  is  in  agreement 
with  immunohistochemical  observations  at  the  electron 
microscopic  level  by  Menco  et  al.  (48),  who  reported  pres¬ 
ence  of  Gaq  immunoreactivity  on  axons  of  vomeronasal 
neurons. 

Pheromone-induced  increases  in  IP3  imply  a  role  for  cal¬ 
cium  in  vomeronasal  signal  transduction  (49).  Although  the 
mechanisms  that  underlie  transduction-excitation  coupling 
in  the  VNO  remain  to  be  elucidated,  several  lines  of  circum¬ 
stantial  evidence  support  a  role  for  calcium  in  this  process. 
High  levels  of  three  calcium  binding  proteins,  calretinin, 
calbindin-D28k,  and  parvalbumin  are  found  in  vomeronasal 
neurons  (50,  51).  In  addition,  patch-clamp  studies  identified 
both  an  L-type  and  T-type  calcium  current  in  rat  vomero¬ 
nasal  neurons,  indicating  a  role  for  calcium  in  neuronal  ex¬ 
citation  (31).  Electron  microscopic  studies  show  that,  in  all 
species  thus  far  examined,  including  pigs,  the  apical  den¬ 
dritic  domes  of  vomeronasal  neurons  are  densely  populated 
with  intracellular  vesicles  (41).  It  is  tempting  to  speculate  that 
these  vesicles  may  serve  as  calcium  storage  depots,  from 
where  calcium  can  be  released  by  IP3.  It  will  be  of  interest  to 
determine  in  future  studies  whether  IP3  receptors  are  located 
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on  these  vesicular  membranes.  In  addition,  we  have  ob¬ 
served  that  incubation  of  microvillar  membranes  with  phor- 
bol  esters  results  in  protein  phosphorylation  (data  not 
shown).  Thus  diacylglycerol  formed  together  with  IP3  may 
also  play  a  role  in  pheromonal  signaling  through  activation 
of  protein  kinase  C. 

Although  behavioral  and  physiological  effects  of  phero¬ 
mones  have  been  well  documented,  mostly  in  rodents  (1-4, 
6-10),  but  also  in  pigs  (5, 13, 14),  unambiguous  identification 
of  pheromones  from  urine  has  been  problematic.  Several 
laboratories  have  reported  the  identification  of  putative 
pheromones  (6,  7),  but  independent  confirmation  of  these 
claims  has  been  notoriously  lacking.  Hitherto,  identification 
of  pheromones  has  depended  on  laborious  and  time-con¬ 
suming  bioassays  that  involve  many  animals  maintained 
under  controlled  environmental  conditions  to  limit  individ¬ 
ual  variation.  This  renders  systematic  isolation  of  phero¬ 
mones  extremely  difficult  and  virtually  impossible  if  the 
biological  response  depends  on  a  blend  of  pheromones,  of 
which  individual  components  may  separate  during  fraction¬ 
ation.  In  addition  to  setting  the  stage  for  further  biochemical 
studies  on  regulation  of  pheromonal  signal  transduction 
pathways  and  transduction-excitation  coupling  in  the  VNO, 
our  experiments  provide  a  biochemical  assay,  i.e.  a  robust 
GTP-dependent  increase  in  IP3,  for  the  future  identification 
of  mammalian  pheromones. 
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ABSTRACT 

Odor-guided  behavior  is  a  polygenic  trait  determined  by  the  concerted  expression  of  multiple  loci. 
Previously,  P-element  mutagenesis  was  used  to  identify  single  P[lArB]  insertions,  in  a  common  isogenic 
background,  with  homozygous  effects  on  olfactory  behavior.  Here,  we  have  crossed  12  lines  with  these 
smell  impaired  (smi)  mutations  in  a  half-diallel  design  (excluding  homozygous  parental  genotypes  and 
reciprocal  crosses)  to  produce  all  possible  66  doubly  heterozygous  hybrids  with  PflArB]  insertions  at  two 
distinct  locations.  The  olfactory  behavior  of  the  transheterozygous  progeny  was  measured  using  an  assay 
that  quantified  the  avoidance  response  to  the  repellent  odorant  benzaldehyde.  There  was  significant 
variation  in  general  combining  abilities  of  avoidance  scores  among  the  smi  mutants,  indicating  variation 
in  heterozygous  effects.  Further,  there  was  significant  variation  among  specific  combining  abilities  of  each 
cross,  indicating  dependencies  of  heterozygous  effects  on  the  smi  locus  genotypes,  ie.,  epistasis.  Significant 
epistatic  interactions  were  identified  for  nine  transheterozygote  genotypes,  involving  10  of  the  12  smi\ oci. 

Eight  of  these  loci  form  an  interacting  ensemble  of  genes  that  modulate  expression  of  the  behavioral 
phenotype.  These  observations  illustrate  the  power  of  quantitative  genetic  analyses  to  detect  subtle  pheno¬ 
typic  effects  and  point  to  an  extensive  network  of  epistatic  interactions  among  genes  in  the  olfactory 
subgenome. 


THE  fundamental  goal  of  quantitative  genetics  is  to 
understand  how  complex  traits  are  shaped  through 
the  interactions  of  multiple  genes  in  different  genetic 
backgrounds  and  under  varying  environmental  condi¬ 
tions.  Perhaps  the  most  complex  category  of  polygenic 
traits  is  represented  by  various  forms  of  animal  behavior. 
Drosophila  melanogaster  presents  an  ideal  model  system  to 
study  the  genetic  basis  of  behavioral  quantitative  traits, 
because  mutations  in  highly  inbred  strains  can  be  easily 
generated,  allowing  control  over  the  segregation  of 
many  individual  loci  that  contribute  to  the  trait  and 
enabling  the  effect  of  each  locus  to  be  studied  indepen¬ 
dently.  We  have  used  odor-guided  behavior  in  D.  melano- 
gasterdiS  a  model  system  to  study  the  quantitative  genetics 
of  behavior. 

Odor-guided  behavior  is  of  special  interest,  because 
the  ability  of  an  organism  to  respond  to  chemical  signals 
from  its  environment  is  essential  for  its  survival  and, 
often,  its  procreation.  Thus,  olfactory  behavior  contrib¬ 
utes  to  individual  fitness  (Mackay  et  al  1996).  In  recent 
years,  considerable  progress  has  been  made  in  elucidat¬ 
ing  the  molecular  mechanisms  that  underlie  odor  recog¬ 
nition,  olfactory  transduction,  and  neural  coding  of  ol¬ 
factory  information  both  in  vertebrates  (reviewed  by 
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Anholt  1993;  Axel  1995;  Buck  1996)  and  in  invertebrate 
model  systems,  such  as  Caenorhabditis  elegans  (Troemel 
et  al  1996;  Sengupta  et  al  1996)  and  lobster  (Fadool 
and  Ache  1992).  However,  the  genetic  basis  of  variation 
in  olfactory  responsiveness  and  the  genetic  mechanisms 
that  shape  behavioral  responses  to  odorants  are  still 
poorly  understood. 

Chemical  mutagenesis  has  been  used  to  induce  muta¬ 
tions  affecting  olfactory  behavior  in  D.  melanogaster, 
mostly  on  the  X  chromosome  (Rodrigues  and  Siddiqi 
1978;  Aceves-Pina  and  Quinn  1979;  Helfand  and 
Carlson  1989;  Lilly  and  Carlson  1989;  McKenna  et 
al  1989;  Ayer  and  Carlson  1992;  Woodard  et  al  1992; 
Lilly  et  al  1994a, b).  This  resulted  in  the  characteriza¬ 
tion  of  a  number  of  genes  that  encode  proteins  likely 
to  participate  in  olfactory  signal  transduction  in  Dro¬ 
sophila,  such  as  smellblind  (an  allele  of  paralytic),  which 
encodes  a  voltage-gated  sodium  channel  (Rodrigues 
and  Siddiqi  1978;  Aceves-Pina  and  Quinn  1979;  Lilly 
and  Carlson  1989;  Lilly  et  al  1994a, b),  norpA,  which 
encodes  a  phospholipase  C  (Woodard  etal  1992;  Riesgo- 
Escovar  et  al  1995),  and  rdgB,  which  encodes  a  phos¬ 
phatidyl  inositol  transfer  protein  (Vipitelic  et  al  1993). 
Although  mutations  in  any  of  these  genes  cause  exten¬ 
sive  impairment  of  olfactory  behavior,  it  is  not  clear  how 
these  genes  contribute  quantitatively  to  variation  in  olfac¬ 
tory  responsiveness  and  how  they  function  in  the  context 
of  the  genetic  background,  ie.,  the  entire  olfactory  subgen¬ 
ome.  Recently,  we  have  identified  14  loci  that  contribute 
to  olfactory  behavior  by  P-element  insertional  mutagenesis 
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in  an  isogenic  strain  (Anholt  et  al  1996).  Identification 
of  these  loci,  designated  smell  impaired  ( smi ),  was  achieved 
using  statistical  and  quantitative  genetic  analysis  of  mea¬ 
surements  of  olfactory  behavior.  These  analyses  are  capa¬ 
ble  of  detecting  small  phenotypic  effects  with  a  resolution 
limited  only  by  sample  size. 

As  the  smi  loci  have  similar  phenotypes,  they  are  likely 
to  be  functionally  related  and  participate  in  common  phys¬ 
iological  and/or  developmental  pathways  that  shape  olfac¬ 
tory  responsiveness.  One  genetic  method  for  identifying 
and  ordering  genes  in  functionally  interacting  groups 
is  to  screen  for  mutations  at  unlinked  loci  that  enhance 
or  suppress  the  mutant  effects  of  a  known  member  of 
the  pathway  ( Garcia- Bellido  1981).  Epistatic  interac¬ 
tions  between  such  genes  can  be  deduced  by  examining 
the  phenotypes  of  the  one-  and  two-locus  genotypes. 
For  independent  loci,  the  phenotypes  of  the  two-locus 
genotypes  are  the  sum  of  the  single-locus  phenotypes; 
i.e.,  the  loci  act  additively.  Departures  from  strict  additivity 
indicate  epistatic,  or  interacting,  loci.  A  simple  test  for 
interaction  that  can  be  used  for  recessive  mutations 
with  large,  qualitative  effects  that  have  similar  loss-of- 
function  phenotypes  and  that  therefore  affect  a  com¬ 
mon  process,  is  to  examine  the  phenotypes  of  double 
mutant  heterozygotes.  Epistasis  is  evident  when  the  dou¬ 
ble  heterozygote  has  the  same  loss-of-function  pheno¬ 
type  as  the  single  homozygous  mutations  as  a  result  of 
combined  haploinsufficiency  of  function.  This  approach 
and  variants  of  it  have  been  used  to  identify  epistatic 
interactions  and  to  identify  new  loci  that  modify  mutant 
phenotypes  of  other  loci  (Botas  et  al  1982;  Belote  et  al 
1985;  Kennison  and  Russell  1987;  Homyk  and  Emerson 
1988;  Tricoire  1988;  Dambly-Chaudiere  et  al  1988). 

Detecting  interactions  between  mutations  with  quanti¬ 
tative  effects  is  more  difficult,  because  the  mutations 
are  not  usually  completely  recessive  (Mackay  etal  1992; 
Lyman  et  al  1996).  Further,  the  background  genotype 
needs  to  be  controlled  to  enable  small  phenotypic  ef¬ 
fects  to  be  perceived  and  to  ensure  any  interactions  are 
due  to  epistasis  between  the  mutations  of  interest,  and 
are  not  confounding  nonadditive  interactions  among 
alleles  segregating  between  the  background  genotypes 
in  which  the  mutations  were  induced.  The  smi  mutations 
are  in  a  common  isogenic  background  and  therefore 
can  be  used  to  detect  epistasis.  We  have  generated  all 
possible  double  heterozygous  hybrids  among  12  inde¬ 
pendent  smi  mutations  that  appear  amenable  to  molecu¬ 
lar  characterization  in  a  diallel  cross  design  (Griffing 
1956),  which  is  the  quantitative  genetic  analogue  of  the 
transheterozygote  test  for  epistasis.  This  approach  is 
based  on  the  assumption  that  reduced  expression  of 
two  independent  P[lArB]-tagged  smi  genes  in  double 
heterozygotic  offspring  may  result  in  quantitative  failure 
to  complement  (Mackay  and  Fry  1996;  Long  et  al 
1996)  if  these  genes  interact. 

Significant  epistatic  interactions  were  identified  for 
nine  transheterozygote  genotypes,  involving  10  of  the 


12  smi  loci.  Interactions  between  eight  of  these  loci  show 
evidence  of  a  web  of  mutually  interactive  genes,  the 
coordinated  expression  of  which  modulates  the  behav¬ 
ioral  phenotype.  These  findings  illustrate  the  power  of 
quantitative  genetic  analyses  to  detect  subtle  phenotypic 
effects  and  indicate  that  phenotypic  determination  of 
odor-guided  behavior  in  D.  melanogaster  depends  quanti¬ 
tatively  on  an  extensive  network  of  genetic  interactions. 

MATERIALS  AND  METHODS 

Generation  of  transheterozygous  PflArB ]  insert  lines:  The 

parental  lines  used  to  generate  double  mutant  heterozygotes 
were  12  homozygous  mi  lines  obtained  by  P-element  mutagen¬ 
esis  of  the  isogenic  Samarkand if06  strain:  smi21F,  smi26D,  smi27E, 
smi28E,  smi35A,  smi45E ,  smi51A ,  smi60E,  smi61A,  smi65A,  smi97B, 
and  smi98B  (Anholt  et  al  1996).  The  mutations  were  named 
according  to  the  cytological  insertion  sites  of  the  P  elements. 
The  12  swzz  lines  were  crossed  in  a  half-diallel  design  (excluding 
homozygous  parental  lines  and  reciprocal  crosses)  to  produce 
all  66  possible  combinations  of  F:  transheterozygous  offspring 
with  two  P  elements  at  different  loci.  Crosses  were  initiated 
at  a  density  of  five  females  of  smi  line  i  and  five  males  of  smi 
line  /  (z  ^  j)  in  plastic  culture  vials.  All  animals  were  reared 
at  25°  on  agar-yeast-molasses  medium. 

Behavioral  assay:  To  quantify  odor-guided  behavior  we  used 
the  simple,  rapid,  and  highly  reproducible  “dipstick”  assay, 
described  previously  (Anholt  et  al  1996).  This  assay  was  cho¬ 
sen  because  it  has  several  advantages  over  other  commonly 
used  assays.  The  Y-maze  assay,  developed  by  Rodrigues  and 
Siddiqi  (1978) ,  is  better  suited  for  measurements  of  attraction 
and  odor  discrimination  than  repulsion  and  is  laborious  for 
large  behavioral  screens.  The  olfactory  jump  assay,  described 
by  McKenna  et  al  (1989) ,  was  in  our  hands  unreliable,  because 
in  contrast  to  Canton-S  flies  for  which  this  assay  was  developed, 
flies  of  both  our  inbred  Samarkand  strain  and  of  substitution 
lines  containing  chromosomes  from  natural  populations 
(Mackay  et  al  1996)  seldom  jumped  in  response  to  repellent 
odorants.  The  “dipstick  assay”  used  in  this  study  and  in  previ¬ 
ous  studies  (Anholt  etal  1996;  Mackay  etal  1996)  is  a  simple, 
rapid,  and  highly  reproducible  statistical  sampling  assay  that 
quantifies  odor-guided  behavior  with  a  resolution  limited  only 
by  sample  size  and,  hence,  can  detect  subtle  olfactory  impair¬ 
ments  (Anholt  et  al  1996).  Previously,  this  assay  has  led  to 
the  identification  of  14  novel  smi  loci.  In  addition,  known 
mutants,  such  as  smellblind  (Rodrigues  and  Siddiqi  1978; 
Aceves-Pina  and  Quinn  1979;  Lilly  and  Carlson  1989; 
Lilly  et  al  1994a, b),  are  immediately  apparent  and  readily 
quantifiable  in  this  assay  (Anholt  et  al  1996). 

After  2-4  hr  of  starvation,  2-10-day  post-eclosion  transhet¬ 
erozygous  progeny  were  tested  for  responsiveness  to  benzalde- 
hyde,  a  repellent  odorant,  exactly  as  described  by  Anholt  et 
al.  (1996).  Briefly,  one  replicate  assay  consisted  of  a  single¬ 
sex  group  of  five  individuals  in  a  test  vial.  The  animals  were 
exposed  to  1%  benzaldehyde  (v/v)  introduced  on  a  cotton 
wool  swab,  and  the  number  of  flies  migrating  to  a  compart¬ 
ment  remote  from  the  odor  source  was  measured  at  5-sec 
intervals,  from  15  to  60  sec  after  introduction  of  the  odor 
source.  The  “avoidance  score”  of  the  replicate  is  the  average 
of  these  10  counts,  giving  a  possible  range  of  avoidance  scores 
between  0  (all  flies  in  the  compartment  near  the  odor  source 
for  the  entire  assay  period)  and  5  (all  flies  in  the  compartment 
away  from  the  odor  source  for  the  entire  assay  period).  For 
each  of  the  66  crosses,  10  replicate  avoidance  score  estimates 
were  obtained  for  each  sex,  for  a  total  of  20  replicates  (100 
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individual  flies)  per  double  heterozygote  genotype,  and  a  total 
sample  size  of  6600  animals. 

Statistical  analyses:  The  avoidance  scores  of  transheterozy- 
gous  genotypes  were  analyzed  by  two-way  analysis  of  variance 
(ANOVA),  with  Genotype  and  Sex  the  fixed  cross-classified 
main  effects.  Sums  of  squares  were  partitioned  into  sources 
(degrees  of  freedom)  attributable  to  Genotype  (65),  Sex  (1), 
Genotype  X  Sex  interaction  (65),  and  Error  (1188).  As  this 
is  a  fixed  effects  model,  the  error  mean  square  was  used  as 
the  denominator  for  all  ETatio  tests  of  significance.  To  analyze 
epistatic  effects  between  smi  loci,  we  could  not  simply  compare 
the  responses  of  double  heterozygotes  with  the  single  heterozy¬ 
gotes  of  smi  lines  with  Sam ,  because  the  effect  of  P[lArB]  insert 
copy  number  (2  vs.  1)  could  be  confounding.  Rather,  the 
correct  control  for  this  analysis  is  measurement  of  the  devia¬ 
tion  from  the  average  of  all  other  transheterozygotes  with  the 
two  single  inserts  being  compared.  Thus,  the  general  combin¬ 
ing  ability  ( GCA)  of  a  mutation  is  its  average  avoidance  score 
as  a  transheterozygote  with  all  other  mutations,  expressed  as 
the  deviation  from  the  overall  mean  (Sprague  and  Tatum 
1942),  and  is  an  estimate  of  the  average  heterozygous  effect 
of  the  mutation  relative  to  the  heterozygous  effects  of  the 
other  mutations.  The  specific  combining  ability  (SCA)  of  a 
transheterozygous  genotype  is  the  difference  between  the  ob¬ 
served  avoidance  score  of  the  genotype,  x (where  i  and  j 
denote  two  different  smi  mutations),  and  the  score  expected 
from  the  sum  of  the  corresponding  GCAs  of  mutants  /  and/ 
The  sums  of  squares  due  to  Genotype  and  Genotype  X  Sex 
were  further  partitioned  into  sources  of  variation  (degrees  of 
freedom)  attributable  to  GCA  (11),  SC  A  (54),  GCA  X  Sex 
(11),  and  SCA  X  Sex  (54).  SCA  effects  are  due  to  variation 
in  heterozygous  effects  that  depend  on  the  genetic  back¬ 
ground  with  respect  to  other  smi  mutations  and  can  only  be 
caused  by  epistatic  interactions. 

This  fixed  effects  half-diallel  corresponds  to  Method  4, 
Model  I  of  Griffing  (1956).  Consequently,  the  GCA  for  each 
smi  mutant  was  estimated  as 

GCA,  =  TJ (n  -  2)  -  ^T/n(n  -  2)  (1) 

where  T{  is  the  sum  of  mean  avoidance  score  values  (averaged 
over  all  replicates)  of  heterozygotes  with  the  tth  mutation, 
£7"  is  twice  the  sum  of  mean  avoidance  score  values  of  all 
heterozygotes,  and  n  is  the  number  of  mutant  lines  (see  also 
Falconer  and  Mackay  1996).  The  SCA  effects  were  com¬ 
puted  using  the  method  of  Griffing  (1956)  for  each  heterozy¬ 
gous  genotype  as 

SCAtj  =  Xij~  (Tt  +  Tj)/(n  -  2)  +  ^T/(n  ~  1  )(n  -  2).  (2) 

The  significance  of  the  overall  GCA,  SCA,  GCA  X  Sex,  and 
SGA  X  Sex  effects  was  tested  using  an  F  variance  ratio  test 
statistic  with  the  error  mean  square  as  the  denominator.  Stan¬ 
dard  errors  of  individual  GCA  and  SC  A  effects  were  computed 
according  to  the  formulae  given  by  Griffing  (1956).  Analyses 
of  variance  and  tests  of  significance  were  calculated  using  SAS 
procedures  (SAS  Institute,  Inc.  1988),  and  GCA  and  SCA 
sums  of  squares  were  computed  using  the  diallel  cross  analysis 
program  of  Schaffer  and  Usanis  (1969). 


RESULTS  AND  DISCUSSION 

^  The  effects  of  1 2  P-element  insertional  mutations  with 

homozygous  effects  on  olfactory  behavior  were  evalu¬ 
ated  in  all  possible  double  heterozygote  combinations, 
in  a  half-diallel  design.  The  mean  avoidance  responses 
to  benzaldehyde,  averaged  over  sexes,  are  shown  for 
each  of  the  66  transheterozygote  genotypes  in  Table  1. 
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TABLE  2 


Analysis  of  variance  of  avoidance  responses  to 
benzaldehyde  of  transheterozygous  smi  lines 

Source 

d.f. 

SS 

F 

P 

Genotype 

65 

110.917 

3.82 

0.0001 

Sex 

1 

4.983 

11.15 

0.0009 

Genotype  X  Sex 

65 

32.493 

1.12 

0.2457 

Error 

1188 

530.695 

The  analysis  of  variance  of  these  data  is  given  in  Table 
2.  The  differences  in  mean  avoidance  responses  among 
the  heterozygous  genotypes  were  highly  significant  (P  = 
0.0001).  There  was  also  significant  sexual  dimorphism 
in  avoidance  response  to  benzaldehyde,  averaged  over 
all  genotypes  (P  =  0.0009) ,  with  a  mean  male  avoidance 
score  of  4.1  and  a  mean  female  score  of  3.9.  Sexual 
dimorphism  for  olfactory  avoidance  response  has  been 
observed  previously  for  homozygous  P-element  inser- 
tional  mutations  (Anholt  et  al  1996)  and  among  a 
sample  of  isogenic  Xand  third  chromosomes  extracted 
from  a  natural  population  and  substituted  into  the  same 
inbred  strain  used  for  P-element  mutagenesis  (Mackay 
et  al  1996).  Interestingly,  both  the  homozygous  P-ele¬ 
ment  insertions  and  the  naturally  occurring  alleles  af¬ 
fecting  olfactory  behavior  had  very  large  genotype  X  sex 
interaction  effects,  indicating  that  there  was  variation  in 
the  magnitude  of  the  sex  dimorphism  of  effects  among 
the  homozygous  genotypes.  However,  the  genotype  X 
sex  interaction  was  not  significant  for  the  double  hetero¬ 
zygote  genotypes;  therefore,  the  sex-specific  effects  ob¬ 
served  previously  are  on  average  recessive. 

Variation  among  the  transheterozygous  genotypes 
can  arise  from  two  sources:  variation  in  mean  heterozy¬ 
gous  effects  of  the  different  mutations,  and  variation 
from  epistatic  interactions.  Because  all  P-element  inser¬ 
tions  are  in  the  same  inbred  strain,  all  genetic  variation 
among  the  genotypes  is  attributable  to  one  of  these  two 
sources,  with  no  confounding  effects  contributed  by 
the  background  genotype.  Classical  diallel  cross  analysis 
enables  us  to  separate  heterozygous  from  epistatic  ef¬ 
fects  by  partitioning  the  variation  among  double  hetero¬ 
zygous  genotypes  into  their  general  ( GCA)  and  specific 
(SCA)  combining  abilities.  As  mentioned  above,  the 
GCA  of  a  mutation  is  an  estimate  of  its  mean  heterozy¬ 
gous  effect  in  the  background  of  each  of  the  other 
mutations.  Estimates  of  the  GCA  of  each  smi  mutation, 
expressed  as  deviations  from  the  overall  mean  of  the 
population  of  heterozygous  genotypes,  are  given  in  Ta¬ 
ble  1.  For  comparison,  also  given  in  Table  1  are  the 
mean  avoidance  scores  of  each  smi  mutation,  at  the 
same  concentration  of  odorant  used  to  assess  trans¬ 
heterozygote  olfactory  behavior  (HOM;  Anholt  et  al 
1996).  All  homozygous  smi  mutations  have  reduced 
avoidance  scores  relative  to  the  transheterozygotes. 


TABLE  3 


Analysis  of  variance  of  general  and  specific  combining 
abilities  of  transheterozygous  smi  lines 


Source 

d.f. 

SS 

F 

P 

Sex 

1 

4.983 

11.15 

0.0009 

GCA 

11 

71.189 

14.49 

<0.0001 

SCA 

54 

39.728 

1.65 

0.0025 

GCA  X  Sex 

11 

2.711 

0.55 

0.87 

SCA  X  Sex 

54 

29.782 

1.23 

0.13 

Error 

1188 

530.695 

Therefore,  negative  GCA  effects  reflect  lower  mean  het¬ 
erozygous  avoidance  scores  and  a  more  mutant  hetero¬ 
zygous  phenotype;  conversely,  positive  GCA  effects  re¬ 
flect  higher  than  average  mean  heterozygous  scores  and 
a  more  wild-type  phenotype.  This  variation  in  GCA 
among  the  smi  mutations  is  highly  significant  (P  < 
0.0001,  Table  3),  and  from  this  we  can  infer  that  all  the 
smi  mutations  are  not  completely  recessive. 

The  overall  mean  avoidance  score  of  the  transhetero¬ 
zygous  genotypes,  3.99  ±  0.08  (Table  1),  is  significantly 
higher  than  that  of  the  Sam;  ry506  strain,  3.65  ±  0.08 
(Anholt  et  al  1996).  Technically,  this  could  be  inter¬ 
preted  as  overdominance  for  olfactory  behavior.  How¬ 
ever,  the  selectable  visible  marker  used  in  this  system 
of  P-element  mutagenesis  is  ry+ ,  and  there  is  a  concern 
that  this  marker  has  a  direct  effect  on  fitness  and  other 
quantitative  traits  relative  to  the  ry~  mutant  background 
of  the  control  strain  (Lyman  etal  1996) .  For  this  reason, 
we  cannot  use  these  data  to  estimate  d,  the  value  of 
the  heterozygote  expressed  as  a  deviation  from  the  mean 
mutant  and  control  strain  value  (Falconer  and  Mackay 
1996)  for  each  smi  mutation.  However,  we  can  estimate 
the  average  degree  of  dominance  of  the  smi  mutations 
from  the  slope  of  the  regression  b  of  GCA  on  homozy¬ 
gous  avoidance  score  of  smi  mutations,  where  each  ho¬ 
mozygous  score  is  expressed  as  a  deviation  from  the 
overall  homozygous  mutant  mean.  The  estimate  of  the 
average  degree  of  dominance  k  is  2  (b  —  0.5)  (Mackay 
et  al  1992;  Lyman  et  al  1996),  where  k  ranges  from 
—  1  (completely  recessive  mutations)  through  0  (strict 
additivity)  to  1  (completely  dominant  mutations).  For 
these  data,  b  =  0.197  (0.01  <  P<  0.05)  and  k  =  -0.605. 
On  average,  the  smi  mutations  are  partially  recessive. 

The  SCA  of  a  pair  of  mutations  reflects  the  extent  to 
which  the  mean  avoidance  score  of  the  double  heterozy¬ 
gote,  expressed  as  a  deviation  from  the  mean  of  the 
total  population  of  heterozygous  genotypes,  departs 
from  that  expected  given  the  sum  of  the  GCAs  of  the 
two  mutant  parents.  Typically,  diallel  crosses  are  made 
among  inbred  lines  that  each  vary  at  a  number  of  loci 
affecting  the  measured  trait,  and  significant  SCA  effects 
can  only  be  attributed  to  nonadditive  interactions  in 
general,  including  dominance  and  epistasis  (Falconer 
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TABLE  4 


Estimates  of  specific  combining  abilities  of 
transheterozygous  smi  lines 


Parent  1 

Parent  2 

SCA 

P 

smi21F 

smi26D 

0.354 

0.009 

smi21F 

smi27E 

0.301 

0.026 

smi21F 

smi45E 

-0.474 

0.001 

smi21F 

smi97B 

-0.369 

0.006 

smi27E 

smi26D 

0.284 

0.036 

smi28E 

smi27E 

-0.319 

0.018 

sm,i51A 

smi97B 

0.275 

0.042 

smi61A 

smi60E 

-0.429 

0.002 

smi65A 

smi51A 

-0.296 

0.029 

and  Mackay  1996).  However,  in  this  experimental  de¬ 
sign  the  genetic  background  has  been  standardized, 
and  SCA  interactions  can  only  result  from  epistasis,  i.e., 
variation  in  heterozygous  effects  of  a  smi  mutation  de¬ 
pending  on  the  genetic  background  with  respect  to 
other  smi  mutations. 

We  observed  highly  significant  SCA  effects  (P=  0.0025, 
Table  3)  for  olfactory  avoidance  among  the  transhetero¬ 
zygote  genotypes.  This  observation  is  not  a  scale  effect. 
The  effect  of  SCA  was  also  highly  significant  if  log, 
square  root,  and  square  transformations  are  applied  to 
the  data  (data  not  shown).  This  suggests  that  epistatic 
interactions  among  loci  affecting  olfactory  behavior  are 
very  common,  because  we  have  sampled  only  a  small 
fraction  of  the  total  number  of  possible  genotypes  at 
12  loci,  each  with  two  alleles  (66/1728  =  3.8%).  To 
determine  which  interacting  mutations  contributed  to 
the  overall  variation  in  SCA,  we  determined  for  which 
transheterozygote  lines  SCA  effects  are  significantly  dif¬ 
ferent  from  zero.  The  results  are  given  in  Table  4.  Nine 
transheterozygous  crosses  reveal  statistically  significant 
epistatic  interactions  between  smi  loci.  In  addition,  the 
smi98B/ smi60E  transheterozygote  has  an  SCA  value 
(—0.251)  that  is  nearly  formally  significant  (P  =  0.063). 
In  five  of  the  nine  statistically  significant  cases,  the  differ¬ 
ence  between  the  observed  and  expected  avoidance 
scores  (SCA)  is  negative;  i.e .,  the  avoidance  response  of 
the  double  heterozygote  is  more  mutant  than  would  be 
expected  given  the  average  degrees  of  dominance  of 
both  parents.  In  four  cases,  the  SCA  estimates  were 
positive,  indicating  better  olfactory  responses  of  the  hy¬ 
brid  offspring  than  expected  from  the  average  heterozy¬ 
gous  effects  of  parental  mutations.  It  should  be  noted 
that  all  of  the  transheterozygotes  show  avoidance  scores 
within  wild-type  range,  i.e.,  complementation,  but  it  is 
the  quantitative  analysis  of  the  degree  of  complementa¬ 
tion  that  reveals  epistatic  effects.  The  negative  and  posi¬ 
tive  interactions  are  quantitative  genetic  analogues  of 
mutations  that  enhance  or  suppress,  respectively,  the 
effects  of  other  mutations  affecting  the  same  phenotype. 


smi65A 
© 

smiSlA 
© 

smi97B 
© 

© 

smi45E  -  smi21F 

© 

smi27E 
© 

smi28E 

Figure  1. — Interaction  diagram  of  smi  loci.  The  0  and  © 
symbols  indicate  epistatic  effects  that  suppress  and  enhance 
the  homozygous  mutant  phenotype,  respectively.  Two  loci, 
smi60E  and  smi61A ,  form  an  independent  pair  with  a  positive 
epistatic  effect  (not  shown). 


The  observed  epistatic  effects  are  quite  large;  the 
mean  of  the  absolute  values  of  significant  SCA  effects 
is  0.34.  This  value  is  of  the  same  magnitude  as  the  mean 
of  the  absolute  values  of  significant  GCA  effects  (0.23), 
and  is  one-half  of  the  environmental  standard  deviation. 
However,  it  is  clear  that  these  “large”  quantitative  effects 
are  very  subtle  in  absolute  terms  and  cannot  be  dis¬ 
cerned  without  quantitative  genetic  analysis  of  the  phe¬ 
notypes,  or  in  variable  genetic  backgrounds.  The  magni¬ 
tude  of  the  epistatic  effects  are  not  necessarily  correlated 
with  the  size  of  the  homozygous  mutant  effects,  smi 
loci  with  relatively  small  effects  on  olfactory  behavior  of 
homozygotes,  e.g.,  smi2  IF  and  smi45E  (Anholt  et  ah  1996), 
produce  large  effects  in  double  heterozygous  progeny. 
We  cannot,  however,  determine  to  what  extent  each 
locus  of  an  interacting  pair  of  loci  contributes  to  the 
observed  epistatic  effect.  Furthermore,  we  do  not  know 
to  what  extent  the  P[lArB]  insertion  limits  the  expression 
of  the  gene  it  affects.  We  predict,  therefore,  that  epistatic 
effects  will  be  stronger  in  double  heterozygotes  that  con¬ 
tain  null  mutations,  such  as  deletions,  at  the  smi  loci. 

The  pattern  of  interactions  observed  is  interesting. 
Of  the  12  smi  loci,  10  interact  with  at  least  one  other. 
Epistatic  interactions  between  eight  smi  loci  can  be  rep¬ 
resented  in  a  simple  interaction  diagram  (Figure  1). 
smi60E  and  smi61A  interact,  but  are  independent  of  the 
others.  It  is  possible  that  smi98B  interacts  with  smi60E 
(the  P  value  of  the  SCA  is  on  the  borderline  of  formal 
statistical  significance),  which  would  place  11  of  the 


1890 


G.  M.  Fedorowicz  et  al. 


12  smi  genes  in  two  interacting  groups.  It  is  somewhat 
surprising  that  the  mutation  that  interacts  most  exten¬ 
sively  with  other  smi  mutants,  smi21F, \  has  itself  very  weak 
homozygous  effects.  The  mutant  phenotype  of  this  gene 
is  only  apparent  at  a  low  concentration  of  benzaldehyde 
and  is  strongly  sexually  dimorphic  (only  females  display 
aberrant  olfactory  responses;  males  are  not  significantly 
different  from  wild  type)  (Anholt  et  al.  1996).  Yet  it 
elicits  strong  interactions  in  transheterozygotes  with 
four  of  the  smi  mutations,  and  in  both  sexes. 

These  loci  represent  only  a  small  sample  of  the  genes 
that  affect  olfactory  behavior.  The  frequency  with  which 
smi  lines  were  detected  in  our  previous  P-element  muta¬ 
genesis  screen  indicated  that  ^4%  of  the  Drosophila 
genome  participates  in  shaping  odor-guided  behavior 
(Anholt  et  al  1996),  which  corresponds  to  a  conserva¬ 
tive  estimate  of  about  400  genes.  Most  likely,  the  ensem¬ 
ble  of  genes  illustrated  in  Figure  1  is  integrated  into 
a  more  extensive  network  of  interactions  within  the 
olfactory  subgenome.  Thus,  loci  that  appear  noninter¬ 
active,  i.e.y  smi35A  and  (possibly)  smi98B,  and  loci  that 
interact  independently  of  the  larger  ensemble,  i.e., 
smi60E  and  smi61A ,  may  prove  to  be  part  of  a  wider 
network  of  interacting  genes  once  more  olfactory  genes 
are  identified.  Such  extensive  epistatic  interactions  be¬ 
tween  smi\ oci  indicate  that  they  form  a  complex  network 
of  genes  that  together  shape  odor-guided  behavior. 

In  recent  years,  other  investigators  have  identified 
olfactory  mutants  in  D.  melanogaster,  mostly  with  muta¬ 
tions  located  on  the  X  chromosome  (Vihtelic  et  al 
1993;  Woodard  et  al  1992;  Riesgo-Escovar  etal  1995; 
Lilly  and  Carlson  1989;  Lilly  et  al  1994a,b),  but 
epistatic  interactions  among  them  and  their  effects  on 
phenotypic  variation  have  not  been  assessed.  Our  obser¬ 
vations  suggest  that  olfactory  genes  identified  on  the 
X  chromosome  might  also  interact  within  functional 
genetic  networks  and  these  possible  interactions  could 
also  include  the  smiloci  described  here.  However,  differ¬ 
ent  genetic  backgrounds  may  render  the  detection  of 
such  epistatic  effects  more  difficult  than  detection  of 
epistasis  among  smi  genes  in  a  coisogenic  background. 

Because  each  of  the  smi  genes  used  in  this  study  is 
tagged  by  a  P element,  it  will,  in  principle,  be  possible  in 
future  studies  to  characterize  their  expression  products 
and  to  obtain  an  understanding  of  the  molecular  basis 
for  the  observed  genetic  interactions.  Moreover,  our 
ability  to  use  coisogenic  P[lArB]-insertion  lines  for  the 
characterization  of  networks  of  interacting  genes  in  the 
olfactory  subgenome  will  enable  the  future  identifica¬ 
tion  of  new  olfactory  genes  by  virtue  of  epistatic  interac¬ 
tions  with  known  smi  genes.  Thus,  these  experiments 
pave  the  road  for  the  use  of  quantitative  genetic  analysis 
of  subtle  phenotypes  as  a  tool  for  targeted  gene  dis¬ 
covery. 
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Olfactomedin  is  a  secreted  polymeric  glycoprotein  of  unknown  function,  originally  discovered  at  the  mucociliary 
surface  of  the  amphibian  olfactory  neuroepithelium  and  subsequently  found  throughout  the  mammalian  brain.  As 
a  first  step  toward  elucidating  the  function  of  olfactomedin,  its  phylogenetic  history  was  examined  to  identify 
conserved  structural  motifs.  Such  conserved  motifs  may  have  functional  significance  and  provide  targets  for  future 
mutagenesis  studies  aimed  at  establishing  the  function  of  this  protein.  Previous  studies  revealed  33%  amino  acid 
sequence  identity  between  rat  and  frog  olfactomedins  in  their  carboxyl  terminal  segments.  Further  analysis,  however, 
reveals  more  extensive  homologies  throughout  the  molecule.  Despite  significant  sequence  divergence,  cysteines 
essential  for  homopolymer  formation  such  as  the  CXC  motif  near  the  amino  terminus  are  conserved,  as  is  the 
characteristic  glycosylation  pattern,  suggesting  that  these  posttranslational  modifications  are  essential  for  function. 
Furthermore,  evolutionary  analysis  of  a  region  of  53  amino  acids  of  fish,  frog,  rat,  mouse,  and  human  olfactomedins 
indicates  that  an  ancestral  olfactomedin  gene  arose  before  the  evolution  of  terrestrial  vertebrates  and  evolved  in¬ 
dependently  in  teleost,  amphibian,  and  mammalian  lineages.  Indeed,  a  distant  olfactomedin  homolog  was  identified 
in  Caenorhabditis  elegans.  Although  the  amino  acid  sequence  of  this  invertebrate  protein  is  longer  and  highly 
divergent  compared  with  its  vertebrate  homologs,  the  protein  from  C.  elegans  shows  remarkable  similarities  in 
terms  of  conserved  motifs  and  posttranslational  modification  sites.  Six  universally  conserved  motifs  were  identified, 
and  five  of  these  are  clustered  in  the  carboxyl  terminal  half  of  the  protein.  Sequence  comparisons  indicate  that 
evolution  of  the  N-terminal  half  of  the  molecule  involved  extensive  insertions  and  deletions;  the  C-terminal  segment 
evolved  mostly  through  point  mutations,  at  least  during  vertebrate  evolution.  The  widespread  occurrence  of  olfac¬ 
tomedin  among  vertebrates  and  invertebrates  underscores  the  notion  that  this  protein  has  a  function  of  universal 
importance.  Furthermore,  extensive  modification  of  its  N-terminal  half  and  the  acquisition  of  a  C-terminal  SDEL 
endoplasmic-reticulum-targeting  sequence  may  have  enabled  olfactomedin  to  adopt  new  functions  in  the  mammalian 
central  nervous  system. 


Introduction 

As  genome-sequencing  projects  proceed  at  an  un¬ 
precedented  rate,  increasing  numbers  of  genes  that  en¬ 
code  proteins  of  unknown  function  are  being  discovered. 
These  include  a  host  of  orphan  receptors,  transcription 
factors  and  proteins  of  which  the  functions  are  open  to 
conjecture.  Uncovering  the  functions  of  many  of  these 
novel  gene  products  would  increase  our  understanding 
and  may  change  our  current  views  of  many  cellular  and 
intercellular  processes.  Evolutionary  analysis  is  of  value 
in  this  endeavor,  since  it  can  reveal  conserved  motifs 
throughout  the  phylogenetic  history  of  a  protein,  thereby 
providing  insights  into  its  potential  function  or,  at  the 
very  least,  identifying  domains  as  targets  for  future  mo¬ 
lecular  studies,  such  as  site-directed  mutagenesis  stud¬ 
ies,  to  evaluate  structure-function  relationships.  Here, 
we  report  a  study  on  the  molecular  evolution  of  olfac¬ 
tomedin,  an  extracellular  matrix  protein  of  yet  unknown 
function,  thought  to  play  an  important  role  in  olfaction, 
neurosecretion,  and  neural  development. 

The  discovery  of  olfactomedin  resulted  from  stud¬ 
ies  aimed  at  the  identification  of  novel  proteins  involved 
in  olfaction.  A  library  of  monoclonal  antibodies  was 
raised  against  chemosensory  cilia  from  the  frog  (Rana 
catesbeiana)  to  find  antigens  uniquely  associated  with 
olfactory  membranes  (Anholt,  Petro,  and  Rivers  1990). 
A  group  of  polypeptides  that  met  this  criterion  appeared 
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to  consist  of  monomeric,  dimeric,  and  polymeric  forms 
of  the  same  protein,  which  was  named  “olfactomedin” 
(Snyder  et  al.  1991).  Immunohistochemical  studies 
showed  that  olfactomedin  is  synthesized  in  Bowman’s 
glands  and  sustentacular  cells  of  the  olfactory  epitheli¬ 
um  and  is  secreted  into  the  viscous  lower  mucus  layer 
that  surrounds  the  apical  dendrites  of  olfactory  neurons 
(Snyder  et  al.  1991;  Bal  and  Anholt  1993).  In  the  frog, 
olfactomedin  immunoreactivity  was  not  detected  in  ol¬ 
factory  nerve  membranes,  nonchemosensory  cilia  from 
respiratory  epithelium,  or  membranes  from  brain,  heart, 
kidney,  or  lung  (Anholt,  Petro,  and  Rivers  1990).  Sub¬ 
sequently,  molecular  cloning  of  olfactomedin  showed  its 
amino  acid  sequence  to  have  no  homologies  to  any  other 
known  protein  (Yokoe  and  Anholt  1993). 

Based  on  the  sequence  of  the  mature  polypeptide 
(448  amino  acids),  together  with  previous  biochemical 
and  immunohistochemical  observations,  a  structural 
model  was  proposed  in  which  cysteines  at  positions  160, 
177,  283,  and  290  form  intramolecular  disulfide  bonds 
dividing  the  molecule  into  central  “head”  and  “neck” 
loops  and  two  “legs”  of  similar  length.  The  N-terminal 
leg  contains  two  nearby  cysteines,  each  forming  disul¬ 
fide  bonds  with  an  adjacent  olfactomedin  molecule.  An 
additional  disulfide  bond  is  made  by  cysteine  377  on  the 
C-terminal  leg,  creating  large  polymers  that  constitute 
the  architecture  of  the  extracellular  matrix  of  the  chemo¬ 
sensory  interface  of  the  olfactory  epithelium  (Yokoe  and 
Anholt  1993).  Olfactomedin  contains  six  potential  N- 
linked  glycosylation  sites  scattered  throughout  its  se¬ 
quence,  and  all  of  them  carry  carbohydrate  moieties  (Bal 
and  Anholt  1993). 

Despite  a  wealth  of  structural  information,  the 
function  of  olfactomedin  remains  unknown.  By  analogy 
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to  functions  of  other  extracellular  matrix  proteins  of  the 
nervous  system,  it  was  proposed  that  olfactomedin 
might  play  a  neurotropic  role  during  development  of  the 
apical  olfactory  dendrite  (Yokoe  and  Anholt  1993). 
Adult  olfactory  neurons  undergo  continuous  replace¬ 
ment  from  a  population  of  neurogenic  stem  cells  (Gra- 
ziadei  and  Monti-Graziadei  1978,  1979;  Farbman  1990; 
Calof  et  al.  1996).  Whereas  questions  about  axon  sorting 
during  formation  of  the  primary  olfactory  projection 
have  received  considerable  interest  (Ngai  et  al.  1993; 
Vassar  et  al.  1994;  Mombaerts  et  al.  1996),  little  atten¬ 
tion  has  been  paid  to  neurotropic  signals  that  induce 
differentiation  of  the  apical  dendrite  to  form  a  knob  with 
chemosensory  cilia  that  mediate  odor  recognition  and 
olfactory  transduction  (reviewed  by  Firestein  1991;  An¬ 
holt  1993;  Axel  1995;  Buck  1996).  The  existence  of 
such  signals  is  inferred  from  the  necessity  of  a  nascent 
dendrite  in  vivo  to  sense  its  arrival  at  the  mucociliary 
surface  and  the  failure  of  isolated  olfactory  neurons  in 
long-term  cell  cultures  to  develop  ciliated  dendritic 
knobs  (Calof  and  Chikaraishi  1989;  Ronnett,  Hester,  and 
Snyder  1991;  Calof  et  al.,  1996).  It  has  been  estimated 
that  olfactomedin  makes  up  5%  of  the  total  tissue  pro¬ 
tein  in  membrane  homogenates  from  frog  olfactory  tis¬ 
sue  (Bal  and  Anholt  1993).  As  a  prominent  component 
of  the  extracellular  mucus  matrix,  olfactomedin  would 
be  in  an  ideal  position  to  deliver  such  a  dendritic  dif¬ 
ferentiation  signal. 

The  hypothesis  that  olfactomedin  might  play  a  role 
in  neuronal  growth  and  differentiation  received  support 
from  the  subsequent  discovery  of  an  olfactomedin  ho¬ 
molog  in  the  mammalian  brain.  In  the  rat  central  ner¬ 
vous  system,  Danielson  et  al.  (1994)  discovered  four 
glycoproteins,  splice  variants  encoded  by  a  single  gene, 
that  are  expressed  throughout  the  brain  and  in  neuroen¬ 
docrine  tissues  such  as  the  anterior  pituitary  and  adrenal 
glands.  Based  on  in  situ  hybridization  patterns,  they  con¬ 
cluded  that  these  glycoproteins  are  produced  by  neurons. 
They  further  observed  that  the  longest  of  these  glyco¬ 
proteins  (the  AMZ  protein)  shares  33%  sequence  iden¬ 
tity  at  its  carboxyl  terminal  end  with  the  carboxyl  ter¬ 
minal  half  of  frog  olfactomedin.  Using  the  rat  AMZ  pro¬ 
tein  as  a  probe,  we  identified  an  olfactomedin  homolog 
from  a  mouse  olfactory  cDNA  library.  Examination  of 
the  GenBank  database  revealed  additional  homologs:  an 
open  reading  frame  from  a  cosmid  (Wilson  et  al.  1994) 
of  Caenorhabditis  elegans,  as  well  as  partial  cDNA  se¬ 
quences  from  fish  pituitary  and  human  brain  libraries. 
To  set  the  stage  for  studies  aimed  at  elucidating  the  func¬ 
tion  of  olfactomedin,  sequences  of  all  known  olfacto¬ 
medin  homologs  were  examined  to  identify  conserved 
motifs. 

This  paper  reports  that  sequence  homologies  be¬ 
tween  the  amphibian  protein  and  its  rodent  counterparts 
are  far  more  extensive  than  previously  appreciated  (Dan¬ 
ielson  et  al.  1994)  and  not  only  cover  the  C-terminal 
region,  but  extend  throughout  the  entire  molecule.  Fur¬ 
thermore,  identification  of  a  distant  olfactomedin  rela¬ 
tive  in  C.  elegans  and  sequence  comparisons  among  ver¬ 
tebrate  olfactomedins  indicate  one  highly  conserved  mo¬ 
tif  near  the  amino  terminus  and  five  conserved  domains 


in  the  carboxyl  terminal  half  of  the  protein.  Evolution 
of  the  N-terminal  half  of  olfactomedin  is  characterized 
by  large  insertions  of  amino  acids,  in  sharp  contrast  to 
the  C-terminal  half,  which  shows  greater  sequence  con¬ 
servation,  suggesting  that  this  extracellular  protein  may 
have  acquired  new  functions  during  its  evolution. 

Methods 

The  following  olfactomedin  amino  acid  sequences 
were  analyzed  (GenBank  accession  numbers  appear  in 
parentheses):  the  sequence  of  frog  ( Rana  catesbeiana) 
olfactomedin  obtained  from  an  olfactory  cDNA  clone 
(L13595;  Yokoe  and  Anholt  1993);  the  sequence  of  the 
rat  ( Rattus  norvegicus)  AMZ  protein,  an  olfactomedin 
homolog  identified  from  brain  cDNA  (U03416;  Daniel¬ 
son  et  al.  1994);  the  translated  partial  sequence  of  an 
olfactomedin  homolog  from  halibut  ( Hippoglossus  hip¬ 
po  glossus)  pituitary  (T23140,  unpublished  sequence, 
Zhiyuan  Gong,  National  University  of  Singapore,  1994); 
the  translated  partial  sequence  of  a  human  ( Homo  sap¬ 
iens)  olfactomedin  homolog  from  an  infant  brain  cDNA 
library  (HI 0467;  unpublished  sequence,  Hillier  et  al., 
the  Washington  University  Merck  EST  Project,  1995); 
and  the  sequence  of  an  olfactomedin  homolog  from  C. 
elegans  (Z81499;  coding  sequence  in  cosmid  F11C3, 
Nematode  Sequencing  Project;  Wilson  et  al.  1994). 

A  probe  corresponding  to  the  carboxyl  terminal  re¬ 
gion  of  the  rat  AMZ  protein,  kindly  donated  by  Dr.  Pa- 
tda  Danielson  (The  Scripps  Research  Institute,  La  Jolla, 
Calif.),  was  used  to  screen  a  mouse  olfactory  cDNA 
library  constructed  in  Lambda  ZAPII  (Stratagene,  La 
Jolla,  Calif.).  Standard  molecular  biological  procedures 
were  performed  as  described  by  Sambrook,  Fritsch,  and 
Maniatis  (1989).  A  clone  with  a  full-length  open  reading 
frame  was  isolated  and  sequenced  at  the  DNA  Sequenc¬ 
ing  Facility  of  Iowa  State  University  (Ames,  Iowa). 

Carboxyl  terminal  sequence  homologies  between 
amphibian  and  rodent  olfactomedins  were  identified 
with  the  BLAST  program  (Altschul  et  al.  1990),  and 
their  alignment  was  verified  by  eye.  To  enable  a  com¬ 
parison  between  vertebrate  olfactomedins  and  a  longer, 
distant  homolog  from  C.  elegans ,  sequences  were  scaled 
to  100%  of  their  lengths.  The  relative  positions  of  con¬ 
served  sequence  motifs  and  posttranslational  modifica¬ 
tion  sites  were  then  placed  along  linear  maps  of  the 
scaled  sequences.  Consensus  glycosylation  sites  and 
leader  sequences  were  verified  with  the  PC/Gene  pro¬ 
gram  from  IntelliGenetics,  Inc.  (Mountain  View,  Calif.). 
Aligned  sequences  were  analyzed  with  the  Molecular 
Evolutionary  Genetics  Analysis  program,  version  1.02 
(Kumar,  Tamura,  and  Nei  1993). 

Results 

Figure  1  shows  the  alignment  of  the  complete  ami¬ 
no  acid  sequences  of  frog  olfactomedin  and  its  rat  ho¬ 
molog.  Previously,  Danielson  et  al.  (1994)  observed 
33%  sequence  identity  between  the  carboxyl  terminal 
regions  of  frog  olfactomedin  and  the  rat  AMZ  protein. 
Reexamination  of  sequence  homologies  between  these 
sequences  indicates  that  homology  between  these  pro- 
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Fig.  1. — Alignment  of  the  amino  acid  sequences  of  frog  and  rat  olfactomedin.  Leader  peptides  are  indicated  in  italics  by  the  boxed  areas. 
The  leader  peptide  of  frog  olfactomedin  was  empirically  determined  (Yokoe  and  Anholt,  1993).  The  leader  peptide  of  rat  olfactomedin  was 
assigned  according  to  Von  Heyne’s  (1986)  rules  as  reported  previously  (Danielson  et  al.  1994).  The  carboxyl  terminal  SDEL  endoplasmic 
reticulum  targeting  sequence  in  rat  olfactomedin  is  in  italics  and  underlined.  Identical  amino  acids  are  highlighted  in  bold  print. 


teins  is  far  more  extensive  than  previously  appreciated, 
although  less  obvious  in  the  N-terminal  halves  of  the 
molecules.  With  the  introduction  of  gaps  in  the  sequenc¬ 
es,  54%  overall  similarity  and  24%  amino  acid  identity 
is  found  scattered  throughout  the  entire  molecule,  albeit 
more  prominently  near  the  carboxyl  terminal  half.  More¬ 
over,  cysteines  essential  for  the  formation  of  character¬ 
istic  homopolymers  are  conserved,  as  are  many  of  the 
glycosylation  sites.  This  can  be  illustrated  by  fitting  the 
sequence  of  rat  olfactomedin  to  the  previously  proposed 
structural  model  of  frog  olfactomedin  (Yokoe  and  An¬ 
holt  1993;  fig.  2).  It  is  clear  that  cysteines  of  the  CXC 


motif  near  the  amino  terminus  are  conserved.  These  cys¬ 
teines  form  neighboring  disulfides  with  an  adjacent  ol¬ 
factomedin  molecule  creating  a  stable  homodimer  (Sny¬ 
der  et  al.  1991;  Yokoe  and  Anholt  1993).  A  cysteine  on 
the  opposite  carboxyl  terminal  leg  is  also  conserved  and 
is  responsible  for  linking  the  homodimers  together  as 
polymers  (Snyder  et  al.  1991;  Yokoe  and  Anholt  1993). 
The  glycosylation  pattern  also  appears  to  be  highly  con¬ 
served.  Carbohydrate  moieties  2,  3,  5,  and  6  of  frog 
olfactomedin  (fig.  2A)  are  found  in  the  same  or  similar 
locations  in  the  rodent  homolog,  corresponding  to  I,  II, 
III,  and  V,  respectively  (fig.  2 B).  Moreover,  folding  of 
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Fig.  2. — Diagrammatic  representation  of  frog  olfactomedin  (A)  and  its  rodent  homolog  ( B ).  Amino  acids  are  represented  as  circles.  Closed 
circles  in  B  indicate  amino  acid  identities  with  corresponding  residues  in  panel  A.  Carbohydrate  moieties  are  numbered  with  Arabic  numerals 
in  A  and  with  Roman  numerals  in  B.  Optimal  alignment  was  achieved  by  the  introduction  of  gaps  in  the  sequences,  which  are  illustrated  as 
sequence  deletions  and  insertions  in  B. 


the  head  region  of  olfactomedin  back  on  itself  across 
the  central  axis  of  the  molecule  would  superimpose  car¬ 
bohydrate  moiety  4  in  figure  2 A  on  moiety  IV  in  figure 
2 B.  Similarly,  intertwining  of  the  amino  terminal  and 
carboxyl  terminal  legs  would  make  the  positions  of  car¬ 
bohydrate  moieties  1  (fig.  L4)  and  VI  and/or  VII  (fig. 
IB)  comparable.  Thus,  despite  considerable  sequence  di¬ 
vergence  between  frog  and  rodent  olfactomedins,  post- 
translational  modifications  that  enable  the  formation  of 
disulfide-linked  polymers  with  distinct  glycosylation 
patterns  are  accurately  preserved.  In  contrast,  intramo¬ 
lecular  disulfides  of  frog  olfactomedin  that  generate  the 
head  and  neck  domains  are  absent  from  the  mammalian 
homolog,  suggesting  that  these  disulfides  are  not  essen¬ 
tial,  although  they  may  help  stabilize  the  conformation 
of  amphibian  olfactomedin. 

To  verify  that  olfactomedin  in  the  rodent,  as  in  the 
frog,  is  also  expressed  in  olfactory  tissue,  a  cDNA  clone 
encoding  olfactomedin  was  isolated  from  a  mouse  ol¬ 
factory  cDNA  library.  Sequence  analysis  of  this  clone 
revealed  that  it  is  closely  related  to  the  rat  AMZ  protein, 
showing  99.1%  amino  acid  sequence  identity.  Of  40  nu¬ 
cleotide  substitutions  between  rat  and  mouse,  36  are  in 
silent  third  codon  positions.  Only  four  amino  acids  were 


different  between  rat  and  mouse  olfactomedin,  at  posi¬ 
tions  106,  183,  302,  and  402  (using  the  numbering  for 
the  rat  sequence  shown  in  fig.  1).  At  these  positions,  the 
mouse  sequence  has  threonine,  serine,  methionine,  and 
arginine,  respectively. 

To  extend  our  analysis  of  the  evolution  of  olfac¬ 
tomedin,  partial  sequences  from  fish,  human,  rat,  and 
frog  olfactomedin  were  aligned.  Sequence  information 
about  the  region  corresponding  to  residues  293  and  343 
of  frog  olfactomedin  (Yokoe  and  Anholt  1993)  is  avail¬ 
able  for  all  five  species.  Sequence  alignments  of  this 
region  are  shown  in  figure  3,  and  the  percentages  of 
sequence  identities  between  these  sequences  are  pre¬ 
sented  in  table  1.  Within  this  region,  seven  invariant 
amino  acid  residues  are  evident,  including  hydrophobic 
residues  leu295,  leu306,  and  ala309  and  a  cluster  of  nega¬ 
tively  charged  residues  between  positions  324  and  331 
forming  the  motif  DXDXXXDE,  in  which  X  is  apolar 
(numbering  is  according  to  the  frog  olfactomedin  se¬ 
quence).  The  occurrence  of  this  invariant  cluster  of  neg¬ 
atively  charged  residues  in  all  four  vertebrate  species 
suggests  that  it  may  play  an  important  structural  and/or 
functional  role. 
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Fig.  3. — Amino  acid  sequence  alignments  of  overlapping  segments  from  partial  sequences  of  fish,  frog,  rat,  mouse,  and  human  olfacto- 
medins.  The  sequences  were  aligned  by  eye  to  maximize  amino  acid  identities.  Identical  amino  acids  are  highlighted  in  bold  print.  Percentages 
of  amino  acid  identities  for  these  sequences  are  listed  in  table  1.  Invariant  amino  acid  residues  are  placed  against  a  shaded  background  and 
indicated  by  arrowheads. 


Amino  acid  conservation  between  frog  and  fish  is 
no  greater  than  that  between  frog  and  rat,  showing 
24.5%  and  32.7%  sequence  identities,  respectively  (ta¬ 
ble  1).  In  contrast,  human  and  rat  olfactomedins  are 
closely  related,  showing  98.1%  identity  over  this  region 
and  96.8%  identity  when  compared  over  a  longer  se¬ 
quence  including  an  additional  71  amino  acids  toward 
the  carboxyl  terminus.  When  nucleotide  differences  be¬ 
tween  rat  and  human  olfactomedin  cDNAs  are  analyzed, 
51  base  substitutions  are  found  within  a  sequence  of  373 
nucleotides.  Of  these  substitutions,  82%  occur  in  third 
codon  positions  and  do  not  alter  the  amino  acid  se¬ 
quence. 

Evolutionary  distances  at  the  DNA  level  were  es¬ 
timated  for  the  sequences  in  figure  3  based  on  the  two- 
parameter  maximum-likelihood  model  of  Kimura 
(1980),  removing  gap  sites  only  in  pairwise  compari¬ 
sons.  Poisson-corrected  amino  acid  substitutions  were 
used  to  estimate  evolutionary  distances  at  the  protein 
level.  Rates  of  evolution  calculated  as  nonsynonymous 
amino  acid  substitutions  per  site  per  109  years  are  pre¬ 
sented  in  table  2.  At  the  amino  acid  level,  there  are  no 
statistically  significant  differences  between  evolutionary 
rates  between  teleosts  and  amphibia  or  between  teleosts 
and  mammals.  From  Table  2,  we  calculate  that  the  evo- 

Table  1 

Sequence  Identities  Among  Teleost,  Amphibian,  and 
Mammalian  Olfactomedins 

Fish  Frog  Rat  Mouse  Human 


Fish .  —  24.5  22.6  22.6  22.6 

Frog .  412  —  32.7  32.7  32.7 

Rat .  412  345  —  99.1  98.1 

Mouse....  412  345  10  —  98.1 

Human  ...  412  345  80  80 


Note. — Numbers  in  the  table  in  the  upper  right  matrix  indicate  the  per¬ 
centage  of  amino  acid  sequence  identity  calculated  for  the  partial  sequences 
shown  in  figure  3.  Numbers  in  the  lower  left  matrix  indicate  estimated  times 
since  divergence  (MYA;  Graybeal  1994). 


lutionary  rate  between  fish  and  frog  is  0.004%/Myr, 
compared  to  0.0015%/Myr  for  the  evolutionary  rate  be¬ 
tween  rat  and  human.  This  suggests  that  evolution  of  the 
protein  within  the  mammalian  lineage  has  slowed  during 
the  last  80  Myr.  However,  it  is  not  clear  whether  these 
differences  in  evolutionary  rates  are  significant,  since 
we  are  not  able,  at  present,  to  distinguish  between  or- 
thologous  and  paralogous  relationships  among  olfacto¬ 
medins. 

Our  data  suggest  that  vertebrate  olfactomedins 
evolved  from  an  ancestral  gene  that  predates  the  evo¬ 
lution  of  terrestrial  vertebrates  about  400  MYA  (table 
1).  The  analysis  presented  in  table  1  led  us  to  hypoth¬ 
esize  that  olfactomedin  homologs  may  also  occur  in  in¬ 
vertebrates.  Indeed,  a  search  through  the  GenBank  da¬ 
tabase  revealed  a  coding  sequence  of  C.  elegans  that 
encodes  a  distant  homolog  of  the  vertebrate  olfactome¬ 
dins.  This  gene  encodes  a  protein  of  654  amino  acids, 
30%  longer  than  the  vertebrate  olfactomedins.  To  iden¬ 
tify  conserved  features  between  vertebrate  olfactome¬ 
dins  and  their  nematode  homolog,  amino  acid  sequences 
of  each  species  were  scaled  to  100%,  and  conserved 
motifs,  sulfhydryls,  and  glycosylation  sites  were 


Table  2 

Rates  of  Nonsynonymous  Substitutions  During  the 
Evolution  of  Olfactomedin 


Species 

Compared 

Time 

Since 

Diver¬ 

gence 

(MYA) 

Kimura 

Two-Parameter 

Model 

Poisson-Corrected 
Amino  Acid 
Substitutions 

Fish  and  frog 

, .  412 

0.46  ±  0.08 

1.66  ±  0.29 

Fish  and  rat . 

.  412 

1.07  ±  0.18 

1.78  ±  0.31 

Rat  and  Human  .  . 

80 

0.64  ±  0.17 

0.12  ±  0.12 

Note. — Evolutionary  rates  are  calculated  in  units  of  substitutions/site  per 
109  years.  Data  presented  are  based  on  a  comparison  of  53  codons  and  are 
arithmetic  means  ±  standard  errors.  Nonsynomymous  nucleotide  substitution 
rates  were  estimated  using  the  model  of  Kimura  (1980). 
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Fig.  4. — Relative  alignments  of  the  amino  acid  sequences  of  Caenorhcibditis  elegans,  frog,  rat,  and  mouse  olfactomedins.  To  accommodate 
the  varying  lengths  of  the  amino  acid  sequences,  they  were  standardized  to  a  scale  of  100%  of  their  lengths  for  each  species.  Relative  locations 
of  cysteines  are  shown  as  red  circles,  and  giycosylation  sites  are  shown  as  blue  diamonds.  Markers  1-6  indicate  areas  of  highly  conserved 
sequences  which  are  identified  below  the  diagrams.  Rat  and  mouse  sequences  are  identical  for  all  of  these  conserved  motifs.  Amino  acids  shown 
in  green  are  identical  for  all  four  species.  Amino  acids  shown  in  magenta  are  the  same  among  three  of  the  four  species  examined. 


mapped  to  the  scaled  sequences  (fig.  4).  Six  highly  con¬ 
served  sequence  motifs  were  identified  at  similar  rela¬ 
tive  locations.  One  conserved  motif  occurs  near  the  N- 
terminus,  and  the  other  five  are  clustered  in  the  C-ter- 
minal  half  of  the  protein.  The  DXDXXXDE  motif  de¬ 
scribed  above  is  contained  in  conserved  motif  3  and  is 
also  present  in  C.  elegans .  The  C  elegans  protein  has 
10  potential  N-linked  giycosylation  sites,  evenly  spaced 
throughout  the  molecule  as  is  the  case  for  vertebrate 
olfactomedins.  Giycosylation  sites  near  motifs  2,  5,  and 
6  occur  at  similar  relative  locations  in  mammalian  ol¬ 
factomedins  and  the  protein  from  C.  elegans  (fig.  4). 
Similarly,  a  single  N-terminal  CXC  motif  is  apparent 
for  the  C.  elegans  homolog,  as  is  a  highly  conserved 
sulfhydryl  in  conserved  motif  5,  proposed  to  participate 
in  the  formation  of  homopolymers  (Yokoe  and  Anholt 
1993).  Thus,  despite  considerable  sequence  diversity 
and  difference  in  polypeptide  length,  the  structural  or¬ 
ganization  of  this  invertebrate  olfactomedin  shows  re¬ 
markable  similarities  to  those  of  its  vertebrate  homologs. 


Discussion 

We  investigated  evolutionary  relationships  among 
vertebrate  olfactomedins  by  comparing  the  sequences  of 
fish,  frog,  rat,  mouse,  and  human  olfactomedins.  Despite 
considerable  amino  acid  sequence  divergence  between 
frog  and  rat,  we  identified  several  invariant  motifs,  in¬ 
cluding  the  CXC  motif  on  the  amino  terminal  leg  and  a 
cysteine  on  the  C-terminal  leg  for  the  formation  of  ho¬ 
mopolymers,  and  NX(T/S)  consensus  giycosylation  sites 
that  preserve  a  distinct  pattern  of  giycosylation  (figs.  1 
and  2;  Bal  and  Anholt  1993;  Yokoe  and  Anholt  1993). 
The  evolution  of  olfactomedin  within  the  mammalian 
lineage  appears  to  be  subject  to  a  strong  evolutionary 
pressure  toward  conservation  of  primary  structure.  Since 
we  cannot,  at  present,  discriminate  paralogous  genes 
from  orthologous  genes,  it  is  difficult  to  draw  solid  con¬ 
clusions  regarding  evolutionary  constraints  from  calcu¬ 
lations  of  evolutionary  rates.  Nonetheless,  a  constant 
rate  of  substitutions  would  predict  approximately  14% 
amino  acid  sequence  divergence  between  human  and  ro- 
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dent  olfactomedins.  Instead,  only  2%  divergence  is  ob¬ 
served,  with  the  vast  majority  of  nucleotide  substitutions 
occurring  in  silent  third  codon  positions.  Invariant  struc¬ 
tural  motifs,  such  as  the  six  motifs  identified  in  figure 
4,  the  CXC  sequence,  the  conserved  sulfhydryl  in  the 
carboxyl  terminal  leg,  and  the  glycosylation  pattern,  are 
likely  to  be  functionally  important  and  can  serve  as  foci 
for  future  studies  exploring  structure/function  relation¬ 
ships. 

Evolutionary  analysis  indicated  that  an  ancestral  ol- 
factomedin  gene  arose  before  the  evolution  of  terrestrial 
vertebrates  and  proceeded  to  evolve  independently  in 
teleost,  amphibian,  and  mammalian  lineages.  Indeed,  a 
distant  homolog  of  olfactomedin  was  identified  in  C. 
elegans.  The  overall  amino  acid  sequence  of  this  inver¬ 
tebrate  protein  was  longer  and  highly  divergent  com¬ 
pared  with  its  vertebrate  homologs.  Conserved  motifs, 
however,  could  be  identified  and  mapped  relative  to  the 
length  of  the  polypeptide  (fig.  4).  This  analysis  revealed 
that  the  protein  from  C.  elegans  shows  remarkable  sim¬ 
ilarities  in  terms  of  conserved  motifs  and  posttransla- 
tional  modification  sites.  Most  of  the  conserved  motifs 
occur  in  the  carboxyl  terminal  half  of  the  protein,  in  line 
with  the  strong  homology  observed  in  this  region  be¬ 
tween  frog  and  rat  olfactomedins  (Danielson  et  al. 
1994).  Thus,  it  appears  that  the  N- terminal  half  of  the 
molecule  has  undergone  more  rapid  evolution,  charac¬ 
terized  by  insertions  and  deletions,  than  the  C-terminal 
region,  which  evolved  slower,  mostly  through  point  mu¬ 
tations,  at  least  during  vertebrate  evolution.  As  dis¬ 
cussed  below,  it  is  conceivable  that  evolution  of  the  N- 
terminal  half  of  olfactomedin  may  have  enabled  the  ac¬ 
quisition  of  new  functions.  The  identification  of  an  ol¬ 
factomedin  homolog  in  C.  elegans  predicts  that  a  similar 
protein  may  also  exist  in  Drosophila  melanogaster.  The 
identification  of  olfactomedin  in  model  organisms  that 
can  be  readily  manipulated  genetically  can  accelerate 
investigations  of  its  function. 

In  addition  to  structural  similarities,  amphibian  and 
mammalian  olfactomedins  also  show  differences.  Cys¬ 
teines  at  positions  283  and  290  of  frog  olfactomedin, 
which  form  intramolecular  disulfides  to  create  the  pro¬ 
posed  head  and  neck  regions  (Yokoe  and  Anholt  1993), 
are  missing  in  its  mammalian  homolog.  This  indicates 
either  that  the  proposed  positions  of  intramolecular  di¬ 
sulfides  in  the  frog  molecule  have  to  be  reevaluated  or 
that  these  intramolecular  disulfide  bonds  are  not  a  priori 
essential  for  the  functional  conformation  of  olfactome¬ 
din.  Furthermore,  the  endoplasmic  reticulum  targeting 
sequence  SDEL  at  the  carboxyl  terminus  of  rat  olfac¬ 
tomedin  does  not  occur  in  the  protein  from  frog.  This 
suggests  either  that  secretion  of  amphibian  olfactomedin 
is  controlled  by  different,  as  yet  unidentified,  targeting 
sequences  within  the  molecule,  or  simply  that  the  poly¬ 
peptide  is  secreted  via  a  default  pathway  which  does  not 
need  to  be  specified  by  an  intracellular  trafficking  signal. 

Previous  studies  have  shown  that  mammalian  ol¬ 
factomedin  occurs  throughout  the  central  nervous  sys¬ 
tem  (CNS)  and  in  neuroendocrine  glands,  such  as  the 
anterior  pituitary  and  adrenal  glands  (Danielson  et  al. 
1994).  Sustentacular  cells  of  the  olfactory  neuroepithel¬ 


ium  and  Bowman’s  glands,  which  are  the  source  of  frog 
olfactomedin  (Snyder  et  al.  1991;  Bal  and  Anholt  1993), 
are  also  derived  from  neural  ectoderm.  The  close  asso¬ 
ciation  of  olfactomedin  with  chemosensory  dendritic 
endings  of  olfactory  neurons  (Snyder  et  al.  1991;  Bal 
and  Anholt  1993),  as  mentioned  earlier,  suggests  that 
olfactomedin  exerts  a  neurotrophic  or  neurotropic  effect 
on  these  neurons.  It  is  of  interest  to  note  that,  thus  far, 
olfactomedin  homologs  have  been  found  in  CNS  neu¬ 
rons  only  in  mammalian  species  (Danielson  et  al.  1994). 
The  fish  homolog  was  identified  from  halibut  pituitary 
cDNA,  and  frog  olfactomedin  has  been  found  only  in 
olfactory  tissue  (Anholt,  Petro,  and  Rivers  1990;  Snyder 
et  al.  1991;  Yokoe  and  Anholt  1993).  The  expression 
pattern  of  the  olfactomedin  homolog  in  C.  elegans  re¬ 
mains  to  be  determined. 

Recently,  the  trabecular-meshwork-inducible  glu¬ 
cocorticoid  response  gene  (TIGR)  in  the  eye  was  shown 
to  contain  olfactomedin  homology  domains.  Mutations 
in  these  olfactomedin-related  domains  of  TIGR  correlate 
closely  with  a  large  fraction  of  cases  of  primary  open- 
angle  glaucoma,  a  leading  cause  of  blindness  and  a  dis¬ 
ease  that  strikes  1  out  of  every  100  individuals  over  the 
age  of  40  years  (Adam  et  al.  1997;  Polanski  et  al.  1997; 
Stone  et  al.  1997).  Like  olfactomedin,  this  protein  is 
associated  with  ciliary  structures,  i.e.,  the  ciliary  rootlet 
in  the  rod  photoreceptor  cell  (Kubota  et  al.  1997).  It  is 
therefore  possible  that  the  conserved  motifs  which  we 
identified  (fig.  4)  may  occur  in  families  of  related  pro¬ 
teins.  A  complete  characterization  of  the  many  proteins 
that  contain  olfactomedin  homology  domains  in  the  spe¬ 
cies  examined  in  this  paper  will  be  needed  to  construct 
a  more  accurate  phylogenetic  history  and  to  determine 
unambiguously  whether  the  olfactomedins  we  analyzed 
are  encoded  by  orthologous  or,  possibly,  paralogous 
genes. 

Failure  of  immunochemical  studies  to  detect  olfac¬ 
tomedin  in  frog  brain  tissue  may  be  due  to  tissue-spe¬ 
cific  glycosylation  and  the  fact  that  unique  carbohydrate 
moieties  dominate  the  immunogenicity  of  olfactomedin 
(Anholt,  Petro,  and  Rivers  1990;  Snyder  et  al.  1991). 
Olfactomedin  is  produced  in  vast  quantities  in  frog  ol¬ 
factory  tissue,  representing  approximately  5%  of  total 
tissue  protein  (Snyder  et  al.  1991).  Thus,  it  is  possible, 
although  unlikely,  that  the  failure  to  detect  olfactomedin 
mRNA  in  frog  brain  previously  by  northern  blotting 
(Yokoe  and  Anholt  1993)  may  have  been  due  to  insuf¬ 
ficient  exposure  of  the  blots  because  of  the  great  abun¬ 
dance  and  easy  detection  of  olfactomedin  message  in 
olfactory  tissue.  A  more  plausible  possibility,  however, 
is  that  in  teleosts  and  amphibia,  olfactomedin  might  ful¬ 
fill  a  function  restricted  only  to  neurosecretory  tissues, 
but  that  during  evolution  of  the  mammalian  brain,  ol¬ 
factomedin  acquired  a  broader  function,  utilized  also  by 
the  nervous  system.  The  expression  of  olfactomedin  in 
neurons  would  necessitate  the  acquisition  of  the  SDEL 
intracellular  trafficking  signal,  and  consolidation  of  a 
new  neural  function  would  constrain  its  mutation  rate. 
Preservation  of  the  characteristic  carbohydrate  groups 
and  the  sulfhydryls  that  are  essential  for  polymerization 
indicates  that  these  posttranslational  modifications  are 
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essential  for  functional  integrity.  The  functions  of  the 
six  conserved  motifs  remain  unknown.  One  possibility 
is  that  they  may  be  necessary  for  protein-protein  inter¬ 
actions  either  with  other  proteins  or  with  adjacent  ol¬ 
factomedin  molecules  of  the  homooligomeric  array. 

The  hypothesis  that  olfactomedin  could  perform 
multiple  functions  in  neural  and  neuroendocrine  tissues 
gains  some  support  from  the  existence  of  splice  variants 
(Danielson  et  al.  1994).  In  the  rat  brain,  olfactomedin 
can  be  expressed  as  either  a  full-length  (Z)  or  a  trun¬ 
cated  (Y)  form.  Whereas  each  of  these  forms  share  the 
same  core  region  (M),  they  can  be  coupled  to  either  one 
of  two  different  amino  termini,  mostly  encoding  leader 
peptide  sequences  (A  and  B).  Thus,  differential  splicing 
gives  rise  to  four  splice  variants  (Danielson  et  al.  1994). 
Two  different  leader  peptides  have  also  been  observed 
for  frog  olfactomedin  (Yokoe  and  Anholt  1993).  The 
significance  of  differences  in  leader  peptide  is  not  clear, 
but  may  bear  on  translational  control  or  intracellular 
trafficking.  Interestingly,  in  the  rat  pituitary  gland,  only 
the  A  form  has  been  observed,  whereas  in  the  adrenal 
gland,  only  the  B  form  is  found  (Danielson  et  al.  1994). 
The  truncated  forms  of  olfactomedin  are  especially  in¬ 
triguing  in  that  they  contain  the  CXC  motif  essential  for 
the  formation  of  olfactomedin  dimers  but  lack  several 
of  the  conserved  motifs  in  the  carboxyl  terminal  region 
of  the  molecule,  including  the  sulfhydryl  group  that  me¬ 
diates  the  formation  of  oligomers.  The  coexpression  of 
truncated  and  full-length  forms  of  olfactomedin  in  the 
mammalian  brain  lends  credence  to  the  notion  that  dif¬ 
ferent  functions  could  be  associated  with  the  amino  ter¬ 
minal  and  carboxyl  terminal  segments  of  the  protein. 

First  identified  in  frog  olfactory  mucus,  it  is  now 
clear  that  olfactomedin  occurs  universally  and  is  widely 
expressed  in  the  mammalian  brain.  The  evolutionary 
analysis  presented  in  this  paper  has  delineated  several 
conserved  motifs  as  targets  for  future  functional  studies 
on  this  intriguing  newly  discovered  extracellular  matrix 
protein. 
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Abstract 

Transmembrane  signaling  events  at  the  dendrites  and  axons  of  olfactory  receptor  neurons  mediate  distinct  functions.  Whereas  odorant 
recognition  and  chemosensory  transduction  occur  at  the  dendritic  membranes  of  olfactory  neurons,  signal  propagation,  axon  sorting  and 
target  innervation  are  functions  of  their  axons.  The  roles  of  G  proteins  in  transmembrane  signaling  at  the  dendrites  have  been  studied 
extensively,  but  axonal  G  proteins  have  not  been  investigated  in  detail.  We  used  immunohistochemistry  to  visualize  expression  of  a 
subunits  of  G0  and  Gi2  in  the  mouse  olfactory  system.  G0  is  expressed  ubiquitously  on  axons  of  olfactory  receptor  neurons  throughout  the 
olfactory  neuroepithelium  and  in  virtually  all  glomeruli  in  the  main  olfactory  bulb.  In  contrast,  expression  of  Gi2  is  restricted  to  a 
sub-population  of  olfactory  neurons,  along  the  dorsal  septum  and  the  dorsal  recess  of  the  nasal  cavity,  which  projects  primarily  to  medial 
regions  of  the  olfactory  bulb,  with  the  exception  of  glomeruli  adjacent  to  the  pathway  of  the  vomeronasal  nerve.  In  contrast  to  the 
overlapping  expression  patterns  of  Gc  and  Gi2  in  the  main  olfactory  system,  neurons  expressing  G0  and  those  expressing  Gi2  in  the 
accessory  olfactory  bulb  are  more  clearly  separated,  in  agreement  with  previous  studies.  Vomeronasal  axons  terminating  in  glomeruli  in 
the  rostral  region  of  the  accessory  olfactory  bulb  express  Gi2,  whereas  those  projecting  to  the  caudal  region  express  G0.  Characterization 
of  the  expression  patterns  of  Gi2  and  G0  in  the  olfactory  projection  is  essential  for  future  studies  aimed  at  relating  transmembrane 
signaling  events  to  signal  propagation,  axon  sorting  and  target  innervation.  ©  1999  Elsevier  Science  B.V.  All  rights  reserved. 

Keywords:  G  protein;  Olfaction;  Olfactory  bulb;  Accessory  olfactory  bulb;  Olfactory  neuron;  Chemotopic  projection 


1.  Introduction 

Odor  recognition  is  essential  for  the  survival  and  pro¬ 
creation  of  most  animals.  Chemical  information  transferred 
via  olfactory  neurons  to  the  olfactory  bulb  is  transformed 
into  a  chemotopic  map  [5,10,11,19].  This  map  is  repre¬ 
sented  by  modular  neural  arrays,  glomeruli,  which  in  the 
mouse  number  approximately  1800  [31]  and  each  of  which 
represents  the  convergent  projection  of  neurons  of  similar 
chemosensory  specificity  (see  Ref.  [16]  for  review). 

The  olfactory  system  contains  two  interacting  compo¬ 
nents:  the  main  olfactory  system  is  dedicated  to  general 
odorant  discrimination,  whereas  the  accessory  olfactory 
system  primarily  processes  chemical  cues  that  guide  social 
behaviors  (see  Ref.  [26]  for  review).  In  the  main  olfactory 
system,  odor  recognition  is  mediated  by  receptors,  which 
belong  to  the  superfamily  of  G  protein  coupled  receptors 
and  are  encoded  by  as  many  as  1000  different  genes 
[2,6,22,40-42].  Two  distinct  families  of  G  protein  coupled 
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receptors  have  also  been  identified  as  putative  pheromone 
receptors  in  chemosensory  neurons  of  the  vomeronasal 
organ  (VNO),  the  chemosensory  organ  of  the  accessory 
olfactory  system  [8,15,25,34]. 

The  dendritic  and  axonal  compartments  of  olfactory 
neurons  fulfill  distinct  functions  for  the  acquisition  of 
chemosensory  information.  Whereas  dendritic  specializa¬ 
tions  mediate  odorant  recognition  and  chemosensory  trans¬ 
duction,  the  axonal  compartment  regulates  signal  propaga¬ 
tion,  axon  sorting  and  target  innervation.  The  roles  of  G 
proteins  in  chemosensory  transduction  at  the  dendritic 
compartments  of  chemosensory  neurons  have  been  studied 
extensively,  but  the  functions  of  axonal  G  proteins  have 
not  been  investigated  in  detail. 

An  important  role  for  axonal  G  proteins  in  mammalian 
chemoreception  is  implicated  by  the  differential  distribu¬ 
tion  of  Gi2  and  G0  in  the  accessory  olfactory  system  (see 
Ref.  [18]  for  review).  The  apical  layer  of  the  VNO  ex¬ 
presses  a  distinct  family  of  about  100  putative  pheromone 
receptors  (VN1  receptors  [8]),  whereas  VNO  neurons  in 
the  basal  layer  express  a  different  family  of  receptors 
(VN2  receptors)  that  resemble  metabotropic  glutamate  re- 
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ceptors  [15,25,34].  Previously,  in  situ  hybridization  studies 
in  the  murine  VNO  [4]  and  immunohistochemical  studies 
in  opossum,  rat  and  mouse  [13,14,18]  showed  that  neurons 
in  the  apical  layer  of  the  VNO  express  Gi2  and  project  to 
the  rostral  region  of  the  accessory  olfactory  bulb  (AOB), 
whereas  neurons  in  the  basal  layer  express  G0  and  project 
to  the  caudal  region  of  the  AOB. 

Here  we  confirm  and  extend  these  previous  observa¬ 
tions.  We  show  that  differential  expression  patterns  of  Gi2 
and  G0  are  not  unique  to  the  AOB,  but  also  occur  in  the 
main  olfactory  bulb  (MOB).  Characterization  of  the  ex¬ 
pression  patterns  of  Gi2  and  G0  in  the  main  olfactory 
system  is  a  prerequisite  for  any  future  studies  aimed  at 
delineating  transmembrane  signaling  events  at  the  axons  of 
olfactory  neurons  and  uncovering  their  relationships  to 
signal  propagation  and/or  odor  coding. 


2.  Materials  and  methods 

2.7.  Animals 

CD-I  mice  were  originally  obtained  from  Charles  River 
Laboratories  (Kingston,  NY)  and  maintained  in  the  breed¬ 
ing  colony  of  Dr.  J.G.  Vandenbergh  (Department  of  Zool¬ 
ogy,  North  Carolina  State  University,  Raleigh,  NC).  Ani¬ 
mals  were  housed  in  IACUC  and  USDA  inspected  and 
approved  facilities  and  cared  for  according  to  the  NIH 
Guide  for  Care  and  Use  of  Laboratory  Animals  (1997). 

2.2.  Antibodies 

Antibodies  against  a  subunits  of  G-proteins  were  ob¬ 
tained  from  Calbiochem  (La  Jolla,  CA).  The  antibody 
against  the  a  subunit  of  Gi2  was  raised  against  the  C- 
terminal  decapeptide  KNNLDCGLF,  which  is  also  found 
in  the  a  subunit  of  Gu.  To  distinguish  Gai2  and  Gan  an 
antibody  against  a  Gau-specific  internal  peptide  was  used. 
The  antibody  against  the  a  subunit  of  G0  was  raised 
against  the  C-terminal  peptide  KNNLKECGLY  of  Gai3, 
which  shares  the  sequence  CGLY  with  Ga0  and  recog¬ 
nizes  the  a  subunits  of  both  G0  and  Gi3.  To  distinguish  G0 
and  Gi3  an  antibody  against  a  Ga  i3-specific  internal  pep¬ 
tide  was  used.  Antibodies  reactive  with  Gai2  and  Ga0  did 
not  cross-react.  Specificity  of  the  antisera  was  confirmed 
by  Western  blotting. 

2.3 .  Western  blotting 

Freshly  dissected  olfactory  bulbs  were  homogenized  in 
ice-cold  0.1  M  phosphate  buffered  saline,  pH  7.3  (PBS) 
with  a  Teflon  homogenizer  and  membranes  were  collected 
by  centrifugation,  washed  once,  recentrifuged  and  resus¬ 
pended  in  PBS.  Protein  was  determined  according  to  Lowry 
et  al.  [23]  using  bovine  serum  albumin  as  standard,  and 
membranes  were  subjected  to  SDS-PAGE,  followed  by 


electrophoretic  transfer  onto  a  nitrocellulose  membrane  [1]. 
Membrane  strips  containing  approximately  20  jxg  protein 
were  incubated  with  a  1000-fold  dilution  of  rabbit  anti¬ 
serum  against  Gai2  or  Ga0.  Bound  antibody  was  visual¬ 
ized  with  Amer sham’s  chemiluminescent  ECL  detection 
system  (Amersham,  Arlington  Heights,  IL).  Migration  dis¬ 
tances  were  calibrated  against  Kaleidoscope  prestained 
molecular  weight  markers  (BioRad,  Richmond,  CA). 

2.4.  Tissue  preparation 

Mice  were  given  a  lethal  injection  of  sodium  pento¬ 
barbital  (50  mg/kg,  i.p.)  and  perfused  intracardially  with 


Fig.  1.  Specificity  of  antibodies.  (A):  Western  blot  showing  immunoreac- 
tive  bands  (arrowheads)  with  antibodies  against  Gai2  and  Ga0  in 
olfactory  bulb  membranes.  (B)-(F):  parasagittal  sections  through  mouse 
olfactory  bulbs  stained  with  normal  rabbit  serum  (B)  or  antisera  against  a 
subunits  of  Gu  (C),  Gi2  (D;  crossreactive  with  Gan),  Gi3  (E),  and  G0  (F; 
crossreactive  with  Gai3).  ONF,  olfactory  nerve  fiber  layer;  GL,  glomeru¬ 
lar  layer;  EP,  external  plexiform  layer;  the  arrowhead  indicates  the  mitral 
cell  body  layer.  Note  the  staining  of  glomeruli  in  panels  D  and  F 
(arrows).  The  scale  bar  in  B  represents  500  |Jim  and  applies  to  B-F. 
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PBS,  followed  by  extensive  perfusion  with  10%  buffered 
neutral  formalin.  The  olfactory  bulbs  and/or  nasal  cavities 
were  dissected  and  fixed  overnight  in  10%  buffered  neutral 
formalin.  Decalcification  of  nasal  tissue  was  performed  for 
3  days  at  ambient  temperature  using  the  formic  acid-sodium 
citrate  method  [24]. 

2.5.  Immunohistochemistry 

Formalin-fixed  and  paraffin-embedded,  5  jx  thick  sec¬ 
tions  through  olfactory  bulbs  and  nasal  tissue  were  de- 
paraffinized  in  xylene  and  rehydrated  through  graded  alco¬ 
hols.  The  sections  were  pretreated  with  0.1%  pepsin 
(Sigma,  St.  Louis,  MO)  in  0.01  M  HC1,  pH  2.3,  for  20  min 
to  facilitate  epitope  access.  Following  pepsin  treatment,  the 
sections  were  blocked  for  30  min  with  BEAT  blocking 
solutions  A  and  B  (Zymed  Laboratories,  San  Francisco, 
CA).  They  were  then  incubated  with  a  250-fold  dilution  of 
normal  rabbit  serum  or  antiserum  in  PBS,  0.05%  Triton 
X-100  (Boehringer  Mannheim,  Indianapolis,  IN)  overnight 


at  4°C.  Bound  antibody  was  visualized  using  the  Histo- 
mouse  SP  kit  from  Zymed  Laboratories.  Following  incuba¬ 
tion  with  the  primary  antibody,  sections  were  washed 
extensively  in  PBS,  0.05%  Triton  X-100,  and  incubated  for 
10  min  at  ambient  temperature  with  affinity  purified  bio¬ 
tinylated  goat  anti-rabbit  antibody.  Following  10  min  incu¬ 
bation  with  horseradish  peroxidase-conjugated  strepta- 
vidin,  antibody  complexes  were  visualized  using  3'-amino- 
9'-ethylcarbazole  as  chromogenic  substrate.  This  generates 
a  red  deposit  at  the  site  of  antibody  binding.  Sections  were 
counterstained  with  hematoxylin,  if  desired,  and  viewed 
and  photographed  under  a  Zeiss  Axiophot  microscope. 

3.  Results 

3.1.  Specificity  of  antibodies 

In  pilot  experiments,  parasagittal  sections  through  olfac¬ 
tory  bulbs  revealed  staining  with  antisera  against  a  sub¬ 
units  of  G  proteins  (Fig.  1).  Before  proceeding  with  a 


Fig.  2.  Immunohistochemical  localization  of  Ga0  in  a  coronal  section  through  the  mouse  nasal  cavity.  S  designates  the  septum  and  the  turbinates  are 
indicated  with  roman  numerals.  Note  that  staining  with  antibodies  against  Got  0  is  apparent  in  all  axon  bundles  in  the  submucosa  (arrowheads).  The  boxed 
area  delineates  the  region  in  which  axon  bundles  stain  with  antiserum  against  Ga0,  but  not  with  antiserum  against  Gai2  (see  Fig.  4).  The  scale  bar 
represents  200  p.m. 
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Fig.  3.  Immunohistochemical  localization  of  Ga0  in  horizontal  sections  through  the  mouse  olfactory  bulb.  Panels  A,  D  and  G  show  low  magnification 
views  of  the  entire  olfactory  bulb  from  a  dorsal  level  (A)  to  progressively  more  ventral  regions  of  the  bulb  (D  and  G).  The  scale  bar  in  G  represents  500 
pum  and  applies  also  to  panels  A  and  D.  Panels  B  and  C  show  higher  magnification  views  of  the  medial  and  lateral  region  of  the  olfactory  bulb  shown  in 
panel  A,  respectively.  Panels  E  and  F  show  higher  magnification  views  of  the  medial  and  lateral  region  of  the  olfactory  bulb  shown  in  panel  D, 
respectively.  Panels  H  and  I  show  higher  magnification  views  of  the  medial  and  lateral  region  of  the  olfactory  bulb  shown  in  panel  G,  respectively.  The 
scale  bars  in  panels  H  and  I  correspond  to  40  p.m  and  apply  also  to  panels  B,  C,  E  and  F.  In  panel  D,  ONF  designates  the  olfactory  nerve  fiber  layer,  GL, 
the  glomerular  layer,  EP,  the  external  plexiform  layer,  and  VN,  the  vomeronasal  nerve.  The  arrowhead  indicates  the  mitral  cell  body  layer.  The  boxed  areas 
in  panels  D  and  G  indicate  the  AOB.  Note  the  ubiquitous  staining  of  the  glomeruli  at  all  three  horizontal  levels  examined,  both  medially  and  laterally 
(arrows). 
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detailed  characterization  of  staining  patterns  we  assessed 
the  specificity  of  our  antisera.  Western  blotting  with  olfac¬ 
tory  bulb  homogenates  reveals  a  single  band  with  an 
apparent  molecular  weight  of  41  kDa  with  the  antiserum 
against  Gai2  and  a  single  band  with  an  apparent  molecular 
weight  of  39  kDa  with  the  antiserum  against  Ga0  (Fig. 
1A).  There  is  no  evidence  for  non-specific  staining  when 
olfactory  bulb  sections  are  incubated  with  normal  rabbit 
serum  (Fig.  IB).  To  evaluate  cross-reactivity  of  the  anti¬ 
serum  against  Gai2  with  Gau,  adjacent  sagittal  sections 
were  stained  with  antiserum  specific  for  Gau  (Fig.  1C) 
and  with  antiserum  against  Gai2  (Fig.  ID).  No  staining  is 
observed  with  the  Gau -specific  antiserum,  whereas  stain¬ 
ing  in  glomeruli  is  evident  with  the  antiserum  against 
Gai2.  The  pattern  of  glomerular  staining  in  this  particular 
region  is  not  uniform,  with  some  glomeruli  staining 
strongly,  whereas  others  are  not  stained  or  show  only  faint 
staining  (Fig.  ID).  Staining  within  individual  glomeruli. 


however,  is  uniform.  To  determine  whether  glomerular 
staining  with  the  antibody  against  Ga0  indeed  is  due  to 
Ga0  rather  than  Gai3  with  which  this  antiserum  crossre- 
acts,  we  stained  adjacent  parasagittal  sections  with  an 
antiserum  that  is  specific  for  Gai3  (Fig.  IE)  and  the 
antiserum  against  Ga0  (Fig.  IF).  Antiserum  raised  against 
an  internal  peptide  of  Gai3  shows  staining  in  the  olfactory 
nerve  fiber  layer  and  faint  staining  in  the  external  plexi- 
form  layer  of  the  MOB,  but  does  not  stain  the  glomeruli. 
In  contrast,  antiserum  against  Ga0  stains  both  the  olfac¬ 
tory  nerve  fiber  layer  and  the  glomeruli.  Although  different 
glomeruli  do  not  stain  with  equal  intensity,  staining  within 
individual  glomeruli  again  is  uniform  (Fig.  IF).  Thus, 
staining  detected  in  MOB  glomeruli  with  the  antiserum 
against  Ga0  indeed  is  due  to  the  presence  of  Ga0.  The 
faint  staining  observed  with  the  Gai3-specific  antiserum 
may  reflect  the  presence  of  Gai3  on  glia.  These  results 
instill  confidence  that  our  antisera  against  Gai2  and  Ga0 


Fig.  4.  Immunohistochemical  localization  of  Gai2  in  a  coronal  section  through  the  mouse  nasal  cavity.  The  section  shown  in  this  figure  is  adjacent  to  the 
one  shown  in  Fig.  2.  S  designates  the  septum  and  the  turbinates  are  indicated  with  roman  numerals.  Note  that  staining  with  antibodies  against  Ga  i2  is 
apparent  in  axon  bundles  in  the  submucosa  (arrowheads)  along  the  dorsal  septum  and  around  the  dorsal  recess  of  the  nasal  cavity.  The  boxed  area 
delineates  the  region  in  which  axon  bundles  do  not  stain  with  antiserum  against  Gai2,  but  stain  with  antiserum  against  Ga0  (see  Fig.  2).  The  scale  bar 
represents  200  ^m. 
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specifically  recognize  the  a  subunits  of  Gi2  and  G0, 
respectively,  in  MOB  glomeruli. 

3.2.  Expression  of  Ga0  in  the  main  olfactory  projection 

Immunohistochemical  staining  of  coronal  sections 
through  the  murine  nasal  cavity  with  an  antiserum  against 
Ga0  shows  extensive  staining  of  axon  bundles  in  the 
submucosa  of  the  olfactory  neuroepithelium  (Fig.  2).  Vir¬ 
tually  all  axon  bundles  along  the  septum  and  all  five 
turbinates  (I,  lb,  II,  lib  and  III)  reveal  expression  of  this  G 
protein. 

As  predicted  from  the  extensive  expression  pattern  of 
Ga0  in  the  olfactory  neuroepithelium,  immunohistochemi¬ 
cal  localization  of  this  G  protein  in  the  olfactory  projection 
to  the  MOB  is  widespread.  Fig.  3  shows  immunohisto¬ 
chemical  staining  patterns  for  Ga0  at  three  horizontal 
levels  through  the  olfactory  bulb  with  views  of  the  entire 
bulb  (A,  D  and  G)  and  higher  magnification  views  of  the 
medial  (B,  E  and  H)  and  lateral  (C,  F  and  I)  aspects  of  the 
bulb.  At  all  three  levels  observed,  from  the  most  dorsal  (A, 
B  and  C)  to  the  most  ventral  (G,  H  and  I),  expression  of 
Ga  0  appears  extensive  and  uniform  within  the  nerve  fiber 
layer  and  glomeruli.  Although  staining  for  Ga0  is  particu¬ 
larly  prominent  in  the  glomeruli,  staining  is  observed  in 
many  of  the  MOB  laminae,  especially  the  granular  cell 
layer,  with  exception  of  the  mitral  cell  body  layer.  This  is 
not  surprising  given  the  prominence  of  G0  as  one  of  the 
major  membrane  proteins  in  the  central  nervous  system 
[17].  Thus,  Ga0  appears  to  be  universally  expressed  in 
olfactory  axons  that  project  to  the  MOB.  In  contrast, 
expression  of  Ga0  in  glomeruli  of  the  AOB  (boxed  areas 
in  Fig.  3D  and  G)  is  not  uniform,  but  localized  to  the 
caudal  region  of  the  AOB  (see  below). 

3.3.  Expression  of  Gai2  in  the  main  olfactory  projection 

In  contrast  to  the  ubiquitous  expression  of  Ga0,  im¬ 
munohistochemical  localization  of  Gai2  in  coronal  sec¬ 
tions  through  the  murine  nasal  cavity  shows  a  more  re¬ 
stricted  expression  pattern  (Fig.  4).  The  expression  pattern 
of  Gai2  appears  nested  within  the  center  of  the  expression 
pattern  of  its  counterpart,  Ga0.  Whereas  virtually  all  axon 
bundles  stain  with  antiserum  against  Gai2  in  the  submu¬ 


cosa  along  the  dorsal  septum  and  around  the  dorsal  recess 
of  the  nasal  cavity,  no  expression  of  this  G  protein  is 
observed  in  olfactory  neurons  originating  from  turbinates  I, 
lb,  lib  and  III  and  the  most  ventral  region  along  the 
septum  (Fig.  4).  Thus,  olfactory  neurons  in  the  more  lateral 
regions  of  the  nasal  cavity  express  only  Ga0,  but  not 
Gai2,  as  indicated  by  the  boxed  areas  in  Figs.  2  and  4. 

In  the  MOB,  the  expression  pattern  of  Gai2  contrasts 
sharply  with  the  uniform  pattern  of  expression  of  Ga0. 
Fig.  5  shows  immunohistochemical  staining  patterns  for 
Gai2  at  three  horizontal  levels  through  the  olfactory  bulb 
with  views  of  the  entire  bulb  (A,  D  and  G)  and  higher 
magnification  views  of  the  medial  (B,  E  and  H)  and  lateral 
(C,  F  and  I)  aspects  of  the  bulb.  At  the  most  dorsal 
horizontal  level  examined,  only  a  subset  of  glomeruli 
expresses  Gai2.  Whereas  different  glomeruli  stain  with 
different  intensity,  within  a  single  glomerulus  the  staining 
intensity  appears  uniform.  Many  glomeruli  immunoreac- 
tive  with  antiserum  against  Gai2  are  located  in  the  medial 
region  of  the  MOB,  whereas  in  the  lateral  region  of  the 
bulb  more  unstained  glomeruli  are  interspersed  with  those 
that  express  Gai2  (Fig.  5B  and  C).  Thus,  the  expression  of 
Gai2  appears  biased  toward  the  medial  aspect  of  the 
MOB.  This  preponderance  of  the  expression  pattern  of 
Gai2  in  the  medial  rather  than  the  lateral  region  of  the 
MOB  is  especially  evident  in  more  ventral  layers  of  the 
bulb  (Fig.  5H  and  I).  Interestingly,  virtually  no  expression 
of  Gai2  is  observed  in  glomeruli  that  border  the  region 
through  which  the  vomeronasal  nerve  penetrates  toward 
the  AOB  (Fig.  5D  and  E).  The  vomeronasal  nerve  stains 
here  prominently  and  projects  to  glomeruli  in  the  anterior 
region  of  the  AOB  (Fig.  5D).  Thus,  the  projection  of 
Ga  i2-expressing  olfactory  neurons  from  the  dorsomedial 
region  of  the  nasal  cavity  fans  out  across  the  medial  aspect 
of  the  MOB,  sparing  the  penetration  corridor  of  the 
vomeronasal  nerve  and  extending  sparsely  to  scattered 
glomeruli  in  the  lateral  region  of  the  MOB. 

3.4.  Expression  of  Ga0  and  Gai2  in  the  accessory  olfac¬ 
tory  projection 

Whereas  in  the  MOB  the  expression  of  Gai2  is  more 
restricted  than  the  distribution  of  Gac,  the  reverse  pattern 
is  evident  in  the  AOB,  where  the  expression  pattern  of 


Fig.  5.  Immunohistochemical  localization  of  Gai2  in  horizontal  sections  through  the  mouse  olfactory  bulb.  Panels  A,  D  and  G  show  low  magnification 
views  of  the  entire  olfactory  bulb  from  a  dorsal  level  (A)  to  progressively  more  ventral  regions  of  the  bulb  (D  and  G).  These  panels  represent  sections 
adjacent  to  their  counterparts  shown  in  Fig.  3.  The  scale  bar  in  G  represents  500  |xm  and  applies  also  to  panels  A  and  D.  Panels  B  and  C  show  higher 
magnification  views  of  the  medial  and  lateral  region  of  the  olfactory  bulb  shown  in  panel  A,  respectively.  Panels  E  and  F  show  higher  magnification  views 
of  the  medial  and  lateral  region  of  the  olfactory  bulb  shown  in  panel  D,  respectively.  Panels  H  and  I  show  higher  magnification  views  of  the  medial  and 
lateral  region  of  the  olfactory  bulb  shown  in  panel  G,  respectively.  The  scale  bars  in  panels  H  and  I  correspond  to  40  jxm  and  apply  also  to  panels  B,  C,  E 
and  F.  In  panel  D,  ONF  designates  the  olfactory  nerve  fiber  layer,  GL,  the  glomerular  layer,  EP,  the  external  plexiform  layer,  and  VN,  the  vomeronasal 
nerve.  The  arrowhead  indicates  the  mitral  cell  body  layer.  The  boxed  areas  in  panels  D  and  G  indicate  the  AOB.  Note  the  heterogeneous  staining  of 
glomeruli.  Gairimmunoreactivity  is  mostly  found  in  medial  glomeruli  of  the  olfactory  bulb  (arrows  in  panels  B  and  H),  although  some  lateral  glomeruli  in 
the  dorsal  region  of  the  bulb  also  stain  (arrows  in  panel  C).  Note  the  virtual  absence  of  expression  of  Gai2  in  panels  D,  E  and  F  where  the  vomeronasal 
nerve  (VN)  travels  along  the  medial  aspect  of  the  MOB  toward  the  AOB  (boxed  area). 
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Gai2  appears  more  extensive  than  that  of  Ga0.  This  is 
especially  evident  in  parasagittal  sections  through  the  AOB 
(Fig.  6A  and  B).  The  arrows  in  panels  A  and  B  of  Fig.  6 
indicate  the  level  of  the  horizontal  sections  shown  in 
panels  C  and  D.  These  horizontal  views  are  higher  magni¬ 
fications  of  the  boxed  areas  of  panels  D  in  Figs.  3  and  5, 
respectively.  The  necklace  glomeruli,  a  group  of  special¬ 


ized  glomeruli  in  the  caudal  region  of  the  MOB  bordering 
the  AOB  [37],  appear  to  express  Ga0,  but  not  Gai2. 
Expression  patterns  in  the  AOB  are  similar  to  those  previ¬ 
ously  observed  in  opossum,  mouse  and  rat  [13,14,18]. 
Neuronal  projections  expressing  Gai2  and  Ga0  are  clearly 
demarcated  in  the  murine  AOB,  with  vomeronasal  axons 
expressing  Gai2  projecting  to  glomeruli  in  the  rostral 
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Fig.  6.  Immunohistochemical  localization  of  Gai2  (panels  A  and  C)  and  Ga0  (panels  B  and  D  in  the  mouse  AOB.  Panels  A  and  B  show  parasagittal 
sections  close  to  the  lateral  border  of  the  AOB.  Panels  C  and  D  are  high  magnification  views  of  the  boxed  areas  of  Figs.  3D  and  5D  and  show  horizontal 
views  of  the  projection  of  the  vomeronasal  nerve  (VN)  at  the  levels  indicated  by  the  arrowheads  in  panels  A  and  B.  Note  that  expression  of  Gai2  is  mostly 
confined  to  the  rostral  region  of  the  AOB,  indicated  by  the  arrow  in  panel  C  and  the  arrowhead  in  panel  D,  whereas  expression  of  Ga0  is  observed  in  the 
caudal  region  of  the  AOB.  The  scale  bar  in  panel  B  represents  25  \xm  and  applies  also  to  panel  A.  The  scale  bar  in  panel  D  represents  100  p,m  and  applies 
also  to  panel  C. 


region  of  the  AOB  and  those  expressing  Ga0  being  re¬ 
stricted  to  the  caudal  region  of  the  AOB  (Fig.  6B  and  D). 


4.  Discussion 

4.1.  Differential  expression  of  G  proteins  in  chemosensory 
projections  to  the  MOB  and  AOB 

We  have  used  immunohistochemistry  to  visualize  the 
expression  patterns  of  the  a  subunits  of  G0  and  Gi2  in  the 
mouse  olfactory  system.  G0  is  expressed  ubiquitously  in 
olfactory  receptor  neurons  throughout  the  olfactory  neu¬ 
roepithelium  and,  consequently,  is  visualized  in  virtually 
all  glomeruli  in  the  MOB  (Figs.  2  and  3).  In  contrast,  the 
expression  of  Gi2  is  restricted  to  a  subset  of  olfactory 
neurons,  located  along  the  dorsal  septum  and  around  the 
dorsal  recess  of  the  nasal  cavity  and  projecting  primarily  to 
medial  regions  of  the  MOB,  with  the  exception  of  glomeruli 
in  the  immediate  vicinity  of  the  vomeronasal  nerve  (Figs.  4 
and  5). 

In  contrast  to  the  overlapping  expression  patterns  of  G0 
and  Gi2  in  the  main  olfactory  system,  the  projection  fields 
of  neurons  expressing  G0  and  those  expressing  Gi2  in  the 
AOB  are  more  clearly  demarcated,  in  agreement  with 
previous  studies  (Fig.  6)  (see  Ref.  [13]  for  review).  In  the 
AOB,  vomeronasal  axons  terminating  in  glomeruli  in  the 


rostral  region  of  the  AOB  express  Gi2,  whereas  axons 
projecting  to  the  caudal  region  of  the  AOB  express  G0. 

Our  results  further  indicate  absence  of  Gu  and  Gi3  in 
primary  olfactory  afferents,  although  the  latter  is  expressed 
in  the  olfactory  nerve  fiber  layer  of  the  MOB  (Fig.  1). 
Previously,  we  localized  Gq/n  to  chemosensory  neurons  in 
the  porcine  VNO,  where  it  is  likely  involved  in  regulating 
inositol- 1, 4, 5-trisphosphate-mediated  responses  to  VNO 
stimuli  [45].  We  did,  however,  not  detect  Gq/11  immuno- 
reactivity  in  the  mouse  MOB  or  AOB.  An  antiserum 
against  Gas  showed  only  weak  staining  in  MOB  glomeruli 
and  staining  patterns  observed  with  this  antiserum  were 
inconclusive.  Recently,  homologous  recombinant  mice 
lacking  the  a  subunit  of  Golf  [20]  were  shown  to  be 
anosmic  [3].  However,  they  retained  a  morphologically 
intact  olfactory  neuroepithelium  and  MOB  [3],  suggesting 
that  Golf  does  not  play  a  role  in  the  formation  of  the 
olfactory  projection. 

The  prominent  presence  of  Gi2  and  G0  in  the  MOB  had 
been  noted  previously  [38],  but  differential  expression 
patterns  of  these  G  proteins  in  populations  of  olfactory 
neurons  throughout  the  MOB  have  not  been  characterized 
before.  Since  these  G  proteins  are  expressed  in  different 
populations  of  neurons  in  the  VNO  and  AOB  [4,13,14,18], 
we  decided  to  focus  on  the  characterization  of  the  expres¬ 
sion  patterns  of  these  two  transduction  proteins.  It  should 
be  noted,  however,  that  we  cannot  exclude  that  some 
glomeruli  may  express  different  G  proteins  that  have  es¬ 
caped  detection  in  this  study. 
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4.2 .  What  are  the  functions  of  G  proteins  on  the  axons  of 
olfactory  and  vomeronasal  neurons? 

Segregation  of  lateral  and  medial  olfactory  projections 
has  been  noted  in  studies  with  monoclonal  antibodies 
against  surface  antigens  [12,28,35,36]  and  N-CAM  iso¬ 
forms  [21].  Characterization  of  segregated  populations  of 
neurons  expressing  different  G  proteins  is  especially  sig¬ 
nificant,  since  these  proteins  are  essential  components  of 
signal  transduction  pathways.  Thus,  the  presence  of  two 
distinct  G  proteins  on  afferent  axons  of  olfactory  neurons 
raises  questions  regarding  the  receptors  these  proteins  are 
coupled  to,  the  signals  they  transduce  and  the  effector 
systems  they  activate. 

G  proteins  on  the  axons  of  olfactory  neurons  could  play 
a  role  either  in  modulation  of  signal  propagation  or  in  the 
assembly  of  convergent  axonal  projections  during  the  for¬ 
mation  and  maintenance  of  the  chemotopic  map.  In  the 
main  olfactory  system,  olfactory  neurons  expressing  the 
same  odorant  receptor  converge  on  two  glomeruli  in  the 
MOB,  one  medial  and  one  lateral  [27,32].  Thus,  the  olfac¬ 
tory  projection  forms  two  chemotopic  representations  in 
each  olfactory  bulb,  one  lateral  map  and  one  medial  map. 
It  is  tempting  to  speculate  that,  whereas  G0  may  play  a 
role  in  axonal  targeting  in  all  olfactory  neurons  in  the 
MOB,  expression  of  Gi2  may  direct  a  subpopulation  of 
neurons  to  form  a  dorsomedial  projection  that  may  define, 
at  least  in  part,  the  medial  chemotopic  map.  Involvement 
of  these  G  proteins  in  the  formation  of  chemotopic  projec¬ 
tions  is  especially  intriguing,  since  odorant  receptors  them¬ 
selves  may  play  a  central  role  in  forming  convergent 
glomerular  projections  [27,44],  As  odorant  receptors  are 
heptahelical  receptors,  axonal  signaling  via  these  receptors 
is  expected  to  be  mediated  via  G  proteins.  Furthermore,  in 
the  VNO  neurons  with  VN1  receptors  express  Gi2,  whereas 
neurons  with  VN2  receptors  express  G0  [15,25,34].  Thus, 
it  is  tempting  to  speculate  that  in  the  main  olfactory 
projection  Gi2  and  G0  may  be  coupled  to  odorant  receptors 
and  play  a  role  in  the  formation  of  the  glomerular  map. 
This  hypothesis  implies  that  odorant  receptors  and  putative 
pheromone  receptors  would  change  their  coupling  speci¬ 
ficity  from  Golf  on  olfactory  cilia  [20]  and  Gq/ll  on  VNO 
microvilli  [45],  respectively,  to  Gi2  or  G0  on  olfactory  and 
vomeronasal  axons.  Such  G  protein  switching  appears 
plausible,  since  (32-adrenergic  receptors  can  switch  their  G 
protein  specificity  upon  phosphorylation  by  cyclic  AMP- 
dependent  protein  kinase  [7].  We  are  well  aware,  however, 
that  the  hypothesis  that  axonal  odorant  receptors  may 
couple  to  G0  or  Gi2  is  at  present  highly  speculative. 

Duncan  et  al.  [9]  showed  that  small  injections  of  an 
anterograde  tracer  in  the  frog  ventral  olfactory  epithelium 
labeled  predominantly  the  lateral  portion  of  the  MOB, 
whereas  injections  in  the  dorsal  olfactory  epithelium  gener¬ 
ated  heavier  labeling  in  the  medial  region.  Mori  provided 
evidence  that  neighboring  glomeruli  in  the  MOB  receive 
inputs  from  olfactory  neurons  with  related  odorant  re¬ 


sponse  properties  [30].  A  recent  fura-2  calcium  imaging 
study  in  mouse  showed  that  retrogradely  labeled  olfactory 
neurons  projecting  to  defined  dorsomedial  or  dorsolateral 
glomerular  regions  respond  to  different  classes  of  odorants 
[5].  Dorsolaterally  projecting  neurons  respond  to  struc¬ 
turally  diverse  odorants,  including  carvone,  eugenol,  cin- 
namaldehyde  and  acetophenone,  whereas  dorsomedially 
projecting  neurons  respond  preferentially  to  organic  acids 
[5].  These  observations  are  consistent  with  earlier  electro- 
physiological  studies  by  Mori  et  al.  [29],  which  showed 
that  mitral  and  tufted  cells  in  the  dorsomedial  region  of  the 
rabbit  MOB  respond  preferentially  to  long  chain  car¬ 
boxylic  acids,  but  not  to  their  corresponding  alcohols. 
These  findings  suggest  that  the  chemotopic  glomerular 
map  may  comprise  dorsolateral  and  dorsomedial  projec¬ 
tions  of  neurons  expressing  functionally  distinct  classes  of 
odorant  receptors. 

As  noted  above,  it  is  possible  that  the  expression  pattern 
of  Gi2,  nested  within  the  projection  of  G0-expressing  neu¬ 
rons,  is  important  for  the  construction  of  lateral  and  medial 
chemosensory  projection  fields.  The  widespread  expres¬ 
sion  of  G0  on  axons  of  olfactory  neurons  is  perhaps  not 
surprising,  since  it  has  long  been  known  that  G0  interacts 
with  GAP43  on  axonal  growth  cones  and  may  play  a  role 
in  linking  transduction  of  extracellular  signals  to  growth 
cone  function  [39].  However,  homologous  recombinant 
mice  deficient  in  the  a  subunit  of  G0  did  not  show  gross 
behavioral  or  neurological  disorders  and  their  optic  nerves 
were  structurally  indistinguishable  from  those  of  heterozy¬ 
gous  litter  mates  [43].  Similarly,  homologous  recombinant 
mice  deficient  in  the  a  subunit  of  Gi2  are  viable  due  to 
compensatory  expression  of  related  G  proteins  [33].  The 
functional  integrity  of  the  chemotopic  maps  and  the  ability 
of  these  animals  to  discriminate  odorants  or  respond  to 
pheromones,  however,  has  not  yet  been  investigated.  More 
detailed  studies  on  these  mice  could  in  the  future  yield 
important  information  about  the  function  of  these  axonal  G 
proteins  in  the  olfactory  system. 
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Legends  to  Figures 


Figure  1:  Diagrammatic  representation  of  the  “dipstick”  assay  for  measurements  of  olfactory 
avoidance  responses.  After  a  1 5s  recovery  period,  the  number  of  flies  in  the  compartment  away  from 
the  odor  source  is  counted  at  5s  intervals  and  the  average  of  10  consecutive  measurements  is 
recorded  as  the  avoidance  score.  The  example  shows  atypical  wild-type  avoidance  response.  Details 
of  the  assay  are  described  in  the  text. 

Figure  2:  Variation  for  avoidance  response  to  benzaldehyde  among  isogenic  chromosome  1  ( X ;  left 
panel)  and  chromosome  3  (right  panel)  substitution  lines  of  Drosophila  melanogaster.  The  male  and 
female  avoidance  scores  of  each  line  are  connected.  Adapted  from  Mackay  et  al.,  1996. 

Figure  3:  Interaction  diagram  of  smi  loci.  The  dotted  and  uninterrupted  lines  indicate  epistatic 
effects  that  enhance  and  suppress  the  homozygous  mutant  phenotype,  respectively.  Two  loci,  smi60E 
and  smi61A,  form  an  independent  pair  with  a  positive  epistatic  effect  (not  shown).  Adapted  from 
Fedorowicz  et  al.,  1998. 

Figure  4:  Unrooted  neighbor-joining  trees  of  18  olfactomedin-related  proteins  (A)  and  their 
olfactomedin  homology  domains  (B).  The  bootstrap  values  are  coded  such  that  a  closed  circle  refers 
to  a  value  of 95  - 1 00%  and  an  open  circle  refers  to  values  between  70  and  94%.  GenBank  accession 
numbers  for  these  protein  sequences  are  as  follows: 
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gi|3695 129|gb|AAC62657. 1 ; 
gi|3  695 1 43  |gb|  AAC62664. 1 ; 
gi|4240 1 28|dbj  |B  AA74844. 1 ; 
gi|4 1 85802|gb|  AAD09 191.1; 
gi[22 1 3659|gb|  AAC98700. 1 ; 
gi|4034486|emb|CAAl  0458.1; 
gi|4 1 64059|gb|AAD05324. 1 ; 
gi|58561  l|sp|Q07081; 
gi|1947128|gb|AAB52933.1; 
gi|3875750|emb|CAB04088.1; 
gi|3024228|sp|Q99784; 
gi|2599125|gb|AAB84058.1; 
gi|3024210|sp|Q62609; 
gi|21 59929|gb|AA447264. 1 ; 
gi|4406679|gb|AAD20056. 1 ; 
gi|1349928|gb|  W53028. 1 ; 
gi|3024209|sp|Q99972; 
gi|3218528|dbj|BAA28767.1 


Figure  5:  Tissue-specific  expression  of  olfactomedin-related  gene  products.  Non-crosshybridizing 
radiolabeled  oligonucleotide  probes  were  hybridized  to  aNorthem  blot  containing  samples  of  human 
heart  (H),  brain  (B),  placenta  (PI),  lung  (Lu),  liver  (Li),  skeletal  muscle  (M),  kidney  (K),  pancreas 
(Pa),  spleen  (S),  thyroid  (Th),  prostate  gland  (P),  testis  (Ts),  ovary  (O),  small  intestine  (Si),  colon 
(C),  and  peripheral  blood  leukocyte  (Pb).  Two  blots  were  hybridized  sequentially  to  different  probes 
and  complete  removal  of  the  previous  probe  was  verified  between  hybridizations.  A  faint  band  of 
hybridization  is  detected  in  pancreas  with  the  TIGR/myocilin  probe.  This  band  is  not  an  artifact  due 
to  incomplete  erasion  of  the  hOlfB  probe,  since  hybridization  with  TIGR/myocilin  was  preceded  by 
hybridization  with  hOlfA  to  this  blot. 


Example: 


Time  (s):  15  20  25  30  35  40  45  50  55  60 

#  flies:  3344443545 


Avoidance  score:  3.9 
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Table  1:  Transposon  insertion  sites  and  candidate  smi  genes  * 


/’-element  insertion 
smi21F 

smi26D 

smi27E 

smi28E 


smi35A 

smi45E 

smi51A 

smi60E 


smi61A 


smi65A 


candidate  genes  at  the  cytological  location 


CG5397  (carboxyl  esterase);  CG4523  (cell  adhesion  protein);  Acap  (adenyl  cyclase  - 
associated  protein);  CG4887  and  CG4896  (RNA  binding  proteins);  CG5001  (heat  shock 
protein);  CG5080  (cytoskeletal  protein);  CG5 1 05  (phospholipase  A2  activating  protein);  up 
to  18  unknown  gene  products,  including  a  putative  odorant  binding  protein. 

CG9493  (protein  phosphatase);  CG9499  and  CG9501  (putative  ion  channels);  CG9507  and 
CG9505  (endopeptidases);  CG9500  (structural  protein);  Tig  (Tiggrin;  extracellular  matrix 
protein);  CG9527  (acyl  coenzyme  A  oxidase  homologue);  CG9508  (neprilysin);  Cpr 
(cytochrome  P450  reductase);  CG9490,  CGI  1573  and  CG9488  (protein  kinases); 
CG9491(c AMP-dependent  Rapl  guanine  nucleotide  exchange  factor);  up  to  16  unknown 
gene  products. 

CG4496  (zinc  finger  transcription  factor);  CG4675  (transport  protein);  Wnt4  (Wnt  oncogene 
analog  4;  up  to  12  unknown  gene  products. 

CG7219  (serpin);  CG7221  (putative  dehydrogenase  enzyme);  CG7367  (lipase  homologue); 
CG7392  (calmodulin  binding  protein  homologue);  CG7424  (ribosomal  protein);  CG7466 
(cell  adhesion  protein);  CG7227  (lysosome  membrane  protein  homologue);  CG7356  (y- 
glutamyl  transferase);  Calo  (calmodulin  binding  protein);  poe  (transmembrane 
protein);CG7586  a2-macroglobulin  homologue);  Trf  (RNA  polymerase  II  transcription 
factor);  CG8668  and  CG8673  (putative  galactosyl  transferase);  up  to  17  unknown  gene 
products  (including  gel  and  belt). 

This  is  the  alcohol  dehydrogenase  ( Adh )  region,  which  has  been  extensively  annotated  by 
Ashbumer  et  al.  (1999).  The  P[lArB]  insertion  site  is  near  wb  (laminin)  and  1(2) 3 4 Fa  (dyrk2 
kinase  homologue). 

Wnt2  (Wnt  oncogene  analog  2);  cro  ( croaker ,  a  courtship  impaired  and  slow  mating  mutant); 
CGI 931  (cytoskeletal  protein);  rdgG  {retinal  degeneration  G ;  unknown  gene  product);  up 
to  5  unknown  gene  products,  including  an  olfactory  receptor. 

CG8151  (RNA  polymerase  II  transcription  factor);  CG8422  (G  protein-coupled  receptor); 
CGI 01 04  (endopeptidase);  CGI 73 85  and  CGI  7390  (zinc  finger  transcription  factors); phyl 
{phyllopod;  nuclear  protein);  cpsf  (component  of  the  cleavage  and  adenylation  specificity 
factor  complex);  Asx  {Additional sex  combs;  chromatin  binding  protein);  ttv  (glucuronyl  N- 
acetylglucosaminyl  transferase  homologue);  CG10110  (RNA  binding  protein);  up  to  23 
unknown  gene  products  (including  the  oho51 ,  auk ,  L ,  and  xen  loci). 
gsb  and  gsb-n  {gooseberry;  RNA  polymerase  II  transcription  factor);  uzip  (integral 
membrane  axon  guidance  protein);  gol  {goliath)  and  Tkr  (zinc  finger  transcription  factors); 
CG2803  (troponin  homologue);  BcDNA:GH04753  (glutathione-S-transferase  homologue); 
CG12850  (transcription  factor);  CG2811  and  CG9358  (putative  ligand  carrier  proteins); 
RpL19  (ribosomal  protein);  CGI 0142,  CG9047  and  ESTS:17F2S  (peptidases);  emp 
(epithelial  membrane  protein);  zip  (non-muscle  myosin);  ETH (ecdysis  triggering  hormone); 
NaCP60E  (sodium  channel  protein);  up  to  20  unknown  gene  products. 

CGI 201,  BcDNA:GH04978  and  Pk61C  (protein  kinases);  CG1216  and  Gyk  (glycerol 
kinase);  CGI  1 869  (putative  microtubule-associated  protein);  CGI 3406  (G  protein-coupled 
receptor);  miple2  (midline/pleiotrophin  family  protein);  Lsply  (larval  serum  protein  1  y- 
subunit);  CG1212  (putative  signal  transduction  protein);  CG7051  (dynein-like  motor 
protein);  CG7036  (putative  transcription  factor);  Mtch  (mitochondrial  carrier  protein);  NitFlit 
(nitrilase  and  fragile  histidine  triad  fusion  protein);  CGI 7 142  (cytoskeletal  structural 
protein);  KazJ  (serine  protease  inhibitor);  up  to  16  unknown  gene  products  (including/vv/). 
CGI 0541  (cytoskeletal  structural  protein);  CGI 0546  (ligand  carrier  binding  protein); 
CG 1 7498  (cell  cycle  regulator);  CG5537  (uracil  phosphoribosyl  transferase);  S6k  (ribosomal 
protein  S-p70-protein  kinase);  vn  (neuregulin-like  protein);  Bjl  (chromatin  binding  protein); 
33-13  and  Ets65A  (DNA  binding  proteins);  CG10486,  CG5592,  CG6600  and  CG10226 
(transport  proteins);  CGI 0487  (receptor  guanylate  cyclase);  CGI 0489  (DNA  replication 
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